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Fabrication  Technology  and  Measurement  of  Coupled 

Quantum  Dot  Devices 


by 

Martin  Burkhardt 

This  report  describes  the  fabrication  and  measurement  of  planar  tunneling  devices.  X-ray 
lithography  was  used  to  define  gate  patterns  in  order  to  achieve  lateral  electrostatic  confinement 
in  a  two-dimensional  electron  gas.  Technologies  were  developed  for  the  printing  of  features  with 
linewidth  of  50  nm  and  below,  a  lithograpic  resolution  which  is  necessary  for  the  fabrication  of 
narrow  tunneling  barriers.  Development  of  technologies  such  as  this  can  also  be  used  for  large- 
scale  fabrication  of  silicon  and  GaAs  devices  and  circuits  with  critical  dimensions  of  100  nm  and 
below. 


Quantum  dots  in  which  the  capacitances  to  the  dot  were  minimized,  were  fabricated  using  high 
resolution  lithography.  Decreased  capacitances  to  the  dot  increase  the  charging  energy  of  a 
quantum  dot,  making  it  possible  to  observe  single  electron  effects  at  elevated  temperatures.  The 
conductance  of  a  device,  featuring  eight  electrodes  to  control  size  and  shape  of  a  quantum  dot, 
was  measured  in  a  Heliox  insertion  probe  at  a  temperature  of  300  mK.  Measurements  of  several 
quantum  dot  sizes  were  performed  and  the  results  were  discussed.  The  same  device  was  biased 
to  produce  two  unequal  quantum  dots  in  series.  The  results  are  discussed  and  compared  with 
theoretical  predictions. 
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Chapter  1 


Introduction 


Since  the  invention  of  the  planar  silicon  process  about  three  decades  ago,  device 
dimensions  have  shrunk  dramatically,  enabling  circuit  designers  to  increase  the  com¬ 
plexity  of  integrated  circuits.  This  progress  was  achieved  through  novel  techniques  in 
device  design  and  layout  as  well  as  lithographic  technologies.  Considering  the  critical 
dimensions  of  present  day  circuits  and  the  wavelength  of  light  used  in  printing  these 
circuits,  it  is  clear  that  lithography  will  be  the  single  most  important  factor  in  the 
miniaturization  of  future  circuits.  Figure  1-1  shows  the  shrinking  of  the  critical  di¬ 
mension  (CD)  over  the  last  few  years  for  successive  generations  of  devices,  including 
projections  for  future  generations.  Each  generation  is  labelled  by  its  corresponding 
dynamic  random  access  memory  (DRAM)  generation.  Up  to  now,  DRAM  develop¬ 
ment  was  the  main  driving  force  for  advances  in  lithography,  but  lithography  for  logic 
applications  is  becoming  equally  demanding.  Requirements  on  overlay  control  -  the 
ability  to  match  different  mask  levels  on  a  substrate  -  are  also  getting  more  stringent 
with  every  new  generation.  Control  of  the  linewidth  is  typically  set  at  10  %  of  the 
critical  dimension.  The  ability  of  a  lithographic  technology,  whether  it  be  optical, 
e-beam  or  x-ray  lithography,  to  meet  these  requirements  will  ultimately  decide  which 
one  will  find  its  way  into  the  manufacturing  environment. 

While  it  is  very  likely  that  optical  lithography  will  be  sufficient  for  the  0.25  ^m 
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Figure  1-1:  Critical  dimension  (CD)  for  successive  DRAM  generations.  Critical  di¬ 
mensions  follow  approximately  a  straight  line  in  this  semi-log  plot.  Other  important 
lithography  parameters,  such  as  overlay  tolerance  and  critical  dimension  control  are 
parallel  to  this  line  (Typical  CD-control  is  10%  of  CD). 

technology,  which  is  anticipated  to  be  in  production  in  1997,  the  next  two  generations 
of  0.18  nm  and  0.12  ixm  will  have  to  be  realized  with  extremely  sophisticated  optics, 
using  deep  UV  lasers  at  KrF-  and  ArF-  wavelengths  close  to  the  diffraction  limit.  At 
such  short  wavelengths,  absorption  within  the  lens  can  cause  heating,  changing  the 
index  of  refraction,  which  in  turn  can  move  the  focal  plane.  It  also  becomes  increas¬ 
ingly  difficult  to  design  lenses  with  high  numerical  aperture  in  the  deep  UV  regime. 
Furthermore,  new  resist  systems  have  to  be  developed  for  every  new  wavelength  gener¬ 
ation  (i.e.  G-line,  I-line,  KrF,  etc.).  The  optical  lithography  community  is  constantly 
developing  new  schemes  to  enhance  the  performance  of  tools  for  a  given  wavelength, 
such  as  off-axis  illumination  and  phase- shifting  masks.  Even  if  it  is  possible  to  design 
such  optics  with  the  necessary  field  size  and  absence  of  distortions  and  aberrations, 
the  requirements  for  the  wafer  and  underlying  relief  structure  will  be  extremely  strict, 
requiring  expensive  processing  steps.  One  of  the  main  problems  is  the  reduced  depth 
of  focus  for  shorter  wavelength,  which  will  be  well  below  0.5 //m  and  will  require  ex¬ 
tremely  fiat  wafers  and  polished  surfaces.  Another  problem  for  optical  lithography  is 
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reflection  off  surfaces  and  interfaces  to  be  exposed,  which  will  require  anti-reflective 
coatings.  It  is  questionable  if  a  solution  which  uses  deep-UV  lithography  will  still  be 
economically  feasible  at  CDs  of  0.18 /xm  and  0.12 /zm.  X-ray  lithography  could  be  in 
a  strong  position  as  the  only  available  alternative  high-throughput,  high  resolution 
lithography.  While  problems  concerning  overlay  and  mask  production  have  to  be  ad¬ 
dressed,  x-ray  lithography  does  have  significant  advantages  in  relief  coverage  [1],  and 
absence  of  reflections. 

For  this  reason,  it  is  essential  to  develop  technologies  needed  for  x-ray  lithography 
and  to  use  these  techniques  for  applications  that  make  use  of  its  advantages.  In  this 
area,  some  ground-breaking  work  has  been  done  at  MIT’s  Nanostructures  Laboratory. 
One  of  the  advantages,  superb  resolution,  has  been  used  to  fabricate  quantum-effect 
devices  [2-10],  high  performance  silicon  devices  [11,12],  and  T-gates  for  MESFET 
devices  [13].  X-ray  lithography  has  also  been  used  to  pattern  large  area  coherent 
deep  submicron  patterns  [14-16],  a  task  difficult  to  achieve  with  other  lithography 
tools  except  for  direct  exposure  using  holographic  lithography. 

It  is  also  questionable  if  it  will  be  possible  to  use  old  paradigms  in  the  search 
for  future  electronics.  The  shrinking  of  device  dimensions  seemed  to  be  limitless  for 
a  long  time,  but  in  the  future  it  might  not  be  feasible  to  continue  to  scale  down 
devices.  There  are  classical  issues  to  be  dealt  with,  e.g.  problems  arising  from  the 
power  dissipation  of  the  devices  which  make  it  very  difficult  to  cool  highly-integrated 
circuits  and  more  fundamental  issues,  e.g.  the  quantization  of  charge,  which  makes  the 
movement  of  small  electron  clouds  between  devices  a  statistically  unreliable  process. 
One  of  the  most  problematic  issues  is  the  so-called  wiring  crisis  in  traditional  circuit 
design.  Even  if  one  develops  technologies  to  implement  extremely  high  integration 
densities,  one  still  has  to  connect  all  the  devices,  leading  to  wiring  problems  on 
the  chip.  One  promising  solution  for  these  problems  is  the  use  of  quantum-effect 
devices[17].  The  nonlinear  characteristics  of  quantum-effect  devices  such  as  resonant 
tunneling  devices  and  single-electron  transistors  could  be  used  for  logic  applications. 
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Arrays  of  quantum-effect  devices  might  be  able  to  do  computation  by  nearest-neighbor 
interaction,  thus  offering  a  solution  to  the  wiring  crisis.  For  this  reason,  it  is  essential 
to  investigate  the  area  of  mesoscopic  physics  and  try  to  transfer  as  much  into  the 
engineering  community  as  necessary  to  create  new  paradigms. 

In  this  report,  I  will  describe  the  fabrication  and  measurement  technology  that  is 
necessary  to  build  planar  tunneling  devices  which  manifest  quantum  effects.  I  have 
chosen  an  unusual  format  in  that  the  core  of  the  report  is  contained  in  only  two 
chapters,  one  chapter  containing  x-ray  lithography  and  device  fabrication,  and  the 
other  device  operation  and  measurement.  This  division  reflects  the  split  of  my  work 
between  fabrication  work,  working  on  x-ray  lithography  and  device  fabrication  in  the 
Nanostructures  Laboratory  on  the  one  hand,  and  measuring  the  fabricated  devices 
in  the  measurement  laboratory  on  the  other.  A  further  division  of  each  of  these  chap¬ 
ters  seemed  artificial  and  not  appropriate.  The  fabrication  chapter  contains  a  short 
introduction  to  x-ray  lithography,  then  describes  the  patterning  and  electroplating 
of  x-ray  masks,  and  finally  details  the  device  fabrication.  The  measurement  chapter 
begins  with  an  introduction  to  the  problem,  briefly  reviewing  reduced-dimensional 
systems.  It  then  describes  the  theory  and  measurement  of  a  quantum  dot  device,  fol¬ 
lowed  by  the  theory  and  measurement  of  a  double  quantum  dot  device.  After  these 
core  chapters,  the  work  is  summarized  and  put  into  perspective  with  general  devel¬ 
opments  in  device  research,  and  logical  continuations  of  this  research  are  proposed. 
In  order  not  to  interrupt  the  train  of  thought,  extensive  parts  describing  details  have 
been  transferred  to  appendices.  They  describe  the  fabrication  sequence  for  devices, 
x-ray  mask  patterning,  the  layout  transfer  from  MIT  to  NRL,  diagnostic  patterns  for 
e-beam  mask  writing,  an  x-ray  lithography  run-sheet,  and  electronic  circuits  used  for 
device  measurement,  x-ray  mask  gap  control,  and  electroplating. 


Chapter  2 


X-Ray  Lithography  and  Device 
Fabrication 


2.1  Introduction 

The  x-rays  typically  used  for  lithography  are  soft  x-rays  with  wavelengths  in  the 
1-nm  region.  In  our  laboratory,  we  use  Cu-L  x-rays  with  a  wavelength  of  1.3  nm. 
Refractive  optics  for  use  in  the  x-ray  regime  do  not  exist  since  all  materials  have 
approximately  the  same  index  of  refraction.  This  is  actually  a  significant  advantage 
of  x-ray  lithography  over  optical  lithography  because  of  the  absence  of  reflections  and 
standing  waves  during  exposure.  The  only  significant  difference  between  materials  at 
these  wavelengths  arises  in  absorption.  X-rays  can  be  projected  using  zone  plates  or 
partially  reflective  multilayer  mirrors,  but  the  approach  taken  by  most  groups  is  to 
go  from  the  farfield  into  the  nearfield  regime  and  use  proximity  x-ray  lithography. 

2.1.1  Diffraction  and  Penumbra 

It  is  well  established  that  the  relationship  between  the  maximum  mask-to-substrate 
gap  G  and  minimum  feature  size  w  in  proximity  x-ray  lithography  is  given  by  the 
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relationship 

G  =  a^  ,  (2.1) 

where  A  is  the  x-ray  wavelength.  In  practice,  the  maximum  value  for  the  parameter 
a  is  equal  to  unity  [18-20],  in  some  cases  1.5  or  even  larger  [21,22].  For  A=lnm,  G 
must  be  5  fim  for  w  <  70  nm.  To  achieve  controlled  gaps  below  5  /im  requires  the 
development  of  new  technologies  such  as  extremely  flat  x-ray  masks  and  the  means 
of  detecting  and  removing  dust  particles  [23]  that  would  otherwise  interfere  with 
the  achievement  of  such  gaps.  Experience  in  other  flelds  (e.g.  magnetic  recording) 
indicates  that  eliminating  dust  and  working  with  small  gaps  is  challenging  but  not  a 
fundamental  problem. 

A  large  gap  can  be  detrimental  to  the  printing  of  flne  features  using  proximity 
x-ray  lithography  for  another  reason:  if  the  source  has  a  finite  extent,  there  exists  a 
penumbral  blur,  6,  of  the  image  cast  onto  the  substrate,  which  can  be  calculated  if 
the  gap,  source  extent,  and  source-to-mask  distance  are  known.  With  this  in  mind, 
one  defines  a  figure  of  merit  /?  such  that  [18] 

5  =  l3w  ,  (2.2) 

where  w  is  the  minimum  feature  size  to  be  printed.  For  a  spatially  coherent  source, 
e.g.  a  plane  wave,  P  vanishes.  These  two  figures  of  merit,  a  and  /3,  are  used- to  describe 
the  effects  of  gap  on  exposure  contrast.  It  turns  out  that  /?  and  a  are  coupled,  that 
is,  for  a  given  w,  a  increases  for  increasing  /3  until  {3  ~  0.5.  To  achieve  a  maximum 
value  of  a  (e.g.  1  <  a  <  1.5),  P  should  be  ~  0.5. 

2.1.2  X-ray  Mask  Standard 

The  x-ray  mask  technology  for  the  fabrication  of  fine  patterns  was  developed  by  previ¬ 
ous  graduate  students  of  the  Nanostructures  Laboratory  [1,24,25]  and  is  described  in 
ref.  [26].  The  present-day  mask  standard  and  most  of  the  work  done  in  our  laboratory 
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membrane 


Figure  2-1:  Schematic  of  the  mask  commonly  used  for  soft  x-ray  proximity  lithography 
at  A=1.3nm.  In  this  case,  the  membrane  consists  of  SiN,  the  absorber  of  Au.  Other 
materials  are  possible,  as  indicated.  Courtesy  of  Bill  Chu. 


is  based  on  this  technology,  with  small  modifications. 

The  mask  consists  of  a  31  mm-diameter  SiN^  membrane  supported  on  a  ~  400  fim- 
thick  Si  ring,  anodically  bonded  to  a  pyrex  support  frame  [26].  A  schematic  of  the 
mask  is  shown  in  Fig.  2-1.  Fabrication  can  be  refined  to  yield  masks  in  which  both 
the  membrane  and  the  mesa  are  optically  flat  to  within  less  than  the  wavelength  of 
visible  light  [27].  A  different  process,  yielding  masks  with  the  same  geometry,  has  been 
shown  to  produce  masks  flat  to  within  500  nm  [28].  The  Si  ring  provides  an  elevation 
of  the  membrane  above  the  plane  of  the  pyrex  ring  of  approximately  400  /xm.  With 
this  “mesa”  configuration,  the  only  part  of  the  mask  that  comes  in  close  proximity 
to  a  substrate,  or  another  mask  during  mask  replication,  is  the  narrow  mesa  rim  and 
the  patterned  membrane  itself. 
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A  plating  base,  consisting  of  5nm  NiCr  and  10  nm  gold  is  evaporated  onto  the 
membrane  to  provide  a  base  for  the  electroplating  of  gold.  Care  is  taken  to  cover  the 
edges  of  the  mesa  with  thick  gold  to  ensure  electrical  conduction  from  the  pyrex  ring 
onto  the  membrane.  This  is  accomplished  with  an  additional  evaporation  of  200  nm 
of  gold  outside  the  membrane. 

2.1.3  Absorber  Material 

The  mask  blank,  described  above,  is  patterned  with  an  absorber  of  10  dB  attenuation 
to  produce  an  x-ray  mask.  The  absorber  is  chosen  to  have  high  atomic  weight  to  keep 
the  absorber  thickness  to  a  minimum  and  for  ease  of  patterning.  Commonly  used 
absorbers  include  gold  [29-32]  and  tungsten  [25,33,34],  used  primarily  in  Europe  and 
the  US,  and  Tantalum,  used  mainly  in  Japan.  While  tungsten  and  tantalum  are 
patterned  using  a  subtractive  dry  etching  process,  gold  is  patterned  via  an  additive 
process  in  which  it  is  electroplated  into  the  interstices  of  a  resist  mold.  Work  described 
in  this  report  is  exclusively  done  with  gold  as  an  absorber,  and  the  patterning  and 
electroplating  of  this  absorber  from  the  mask  blank  to  the  patterned  mask  is  described 
in  detail  in  succeeding  sections.  The  mask  patterning  process  is  described  in  detail 
in  Appendix  A. 


2.2  Mask  Flatness  and  Particle  Control 

Ideally,  the  mesa  rim  is  optically  flat  [28].  For  the  older-style  masks  used  in  most  of 
the  experiments  reported  here,  the  mesa  rim  deviates  from  a  perfect  flat  by  about  1-3 
pm.  This  nonflatness  sets  an  initial  minimum  gap  for  experiments  in  which  the  mask 
is  brought  in  proximity  to  a  wafer  or  another  mask.  It  is  not  a  great  concern  when 
the  mask  is  used  to  print  onto  substrates  that  are  small  compared  to  the  diameter  of 
the  membrane.  If  one  copies  one  mask  onto  another,  the  nonflatness  of  the  mesa  rim 
provides  an  initial  gap  of  2  to  6  pm.  For  an  ideal,  extremely  flat  mesa  rim,  studs  can 
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be  placed  on  the  rim  to  establish  an  initial  gap. 

To  check  for  the  presence  of  particles  and  for  flatness  of  an  x-ray  mask,  a  glass 
optical  flat  with  a  ~5  /zm-thick  spacer  made  of  photoresist  is  used,  as  shown  in  Fig.  2- 
2.  Any  patterned  or  unpatterned  clear  field  x-ray  mask  can  be  placed  onto  this  optical 
flat.  The  nonflatness  of  the  mask  can  easily  be  checked  in  green  light  by  observing 
the  interference  fringes  caused  by  the  gap  between  the  mask  and  flat.  The  gap  can  be 
checked  by  counting  the  fringes  that  appear  by  moving  the  point  of  observation  from 
a  90°  angle  to  a  45°  angle.  The  gap  is  then  approximately  the  number  of  fringes  minus 
one  in  microns.  This  gap  measurement  technique  was  described  by  Schattenburg  and 
others  [35].  The  gap  between  the  mask  and  flat  should  be  verified  to  be  approximately 
the  thickness  of  the  photoresist.  The  gap  typically  turns  out  to  be  a  little  larger  than 
the  resist  thickness  because  of  the  nonflatness  of  the  mesa. 

Particles  on  the  mask  or  flat  that  are  larger  than  the  mask-to-flat  gap  can  easily  be 
detected  by  observing  the  fringe  pattern  in  green  light.  Masks  are  frequently  checked 
using  this  technique  to  detect  possible  particles  on  membranes  or  mesas  before  x-ray 
exposures.  Figure  2-2  shows  a  schematic  and  a  photograph  of  a  mask  on  an  optical 
flat,  illuminated  with  green  light. 


2.3  Patterning  of  Masks 

2.3.1  Process  Outline 

A  typical  sequence  for  the  fabrication  of  the  gate  level  of  a  quantum  effect  device  is 
outlined  in  Fig.  2-3.  The  fabrication  process  evolves  from  the  mother  mask  (a),  as 
written  by  electron  beam  lithography,  to  the  daughter  mask  (b),  to  the  device  (c), 
which  in  this  figure  is  a  monitor  sample  consisting  of  Ti/Au  (20/80  nm)  on  silicon. 
The  mother  mask  (e.g.  Fig.  2-3(a))  is  exposed  with  a  50keV  electron  beam  using 
typically  between  220  and  300  nm-thick  PMMA  as  a  resist.  The  mask  is  then  elec¬ 
troplated  with  200  nm  of  gold.  After  the  completion  of  the  mother  mask,  a  daughter 
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Figure  2-2:  Schematic  and  photograph  of  particle  detection  scheme.  The  x-ray  mask 
is  placed  on  the  optical  flat  to  check  for  particles  and  flatness. 

mask  (e.g.  Fig.  2-3(b))  is  made  using  x-ray  lithography  at  a  gap  between  1  and  5/im. 
Our  x-ray  source  is  a  water-cooled  electron  bombardment  target  emitting  Cu-L  x- 
rays.  For  masks  containing  large  features  (>0.1/im),  a  gap  of  about  5/im  may  be 
used.  Masks  containing  very  fine  features  are  replicated  at  a  gap  below  1.5 /xm.  Be¬ 
cause  x-ray  exposure  is  used  instead  of  e-beam  exposure,  the  process  latitude  in  the 
development  (i.e.  development  time  and  developer  temperature)  is  much  larger.  The 
daughter  mask  is  then  used  to  expose  substrates  in  intimate  contact  (Fig.  2-3  (c)). 
This  intimate  contact  is  achieved  by  creating  a  partial  vacuum  between  the  substrate 
and  the  mask,  as  described  in  earlier  papers  [9,36]. 

2.3.2  Mother  Mask 

The  process  of  writing  high  density  patterns  with  an  e-beam  is  rather  time  consum¬ 
ing.  This  makes  e-beam  lithography  unsuitable  for  large-scale  manufacturing.  In  our 
case,  however,  a  single  x-ray  mask  is  written  by  e-beam  lithography.  This  mask  is 
then  replicated  onto  a  daughter  mask,  which  is  then  used  in  x-ray  lithography.  X-ray 
lithography  is  a  parallel  rather  than  a  serial  exposure  process  (as  is  e-beam  lithogra- 
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Figure  2-3:  a:  mother  mask  (200  nm  gold  thickness);  b:  daughter  mask  (200  nm  gold 
thickness);  c:  finished  gate  pattern  (20/80  nm  Ti/Au)  of  a  device  with  a  point  con¬ 
striction  and  a  tunneling  barrier  on  each  side.  The  center-to-center  distance  between 
the  gates  is  200  nm,  the  linewidth  is  about  70  nm. 

phy)  and,  consequently,  can  have  a  higher  throughput.  Our  laboratory  often  writes 
even  large  area  patterns  such  as  pads  and  bounding  boxes  using  e-beam  lithography. 
If  this  is  not  tolerable,  one  can  utilize  a  mix-and  match  scheme  to  pattern  pads  with 
optical  lithography  and  write  only  fine  features  by.  e-beam  [27]. 

After  a  mask  pattern  has  been  designed,  it  is  converted  into  a  data  format  suitable 
for  a  JEOL  JBX-5DII  electron-beam  lithography  system.  The  mother  mask  is  written 
at  the  Naval  Research  Laboratories  (NRL)  in  Washington  D.C.  The  pattern  data  are 
transferred  electronically  to  NRL  by  means  of  the  Internet.  The  specifics  of  the 
transfer  process  and  the  patterning  have  been  described  in  previous  theses  [24,37], 
and  an  updated  description  is  given  in  Appendix  B.  The  e-beam  writing  at  NRL 
uses  an  acceleration  voltage  of  50  keV  and  different  beam  currents  for  different-sized 
features.  For  coarse  features  (~10//m),  a  beam  current  of  3nA  and  an  areal  dose  of 
approximately  400//C/cm^  is  used.  For  submicron  features,  a  beam  current  of  15  to 
50  p A  and  an  areal  dose  of  400  fxC/cm^  is  typically  used.  Critical  freatures  in  lateral 
tunneling  devices  typically  require  as  narrow  a  linewidth  as  possible.  For  the  finest 
linewidths,  those  below  about  100  nm,  a  line  dose  of  as  low  as  InC/cm  is  used.  In 
this  case,  the  electron  beam  sweeps  across  the  length  of  the  line  only  once,  exposing  it 
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in  a  single  pass.  Typically,  the  same  die  is  written  at  several  locations  on  a  membrane 
with  different  areal  and  line  doses  assigned  to  each  die  to  ensure  that  at  least  one  die 
satisfies  the  design  with  acceptable  linewidth  results. 

After  the  e-beam  exposure,  gold  is  plated  up  to  a  thickness  of  200  nm,  which 
provides  a  contrast  of  approximately  10  dB  for  the  Cu-L  line  x-rays  used  in  our  labo¬ 
ratory.  To  plate  to  a  thickness  of  200  nm,  one  should  use  a  resist  thickness  of  at  least 
220  nm.  This  seems  to  make  the  lithographic  definition  of  fine  lines  more  difficult  than 
in  a  process  using  direct  write  e-beam  which  typically  uses  thin  resist  (~50-100nm). 
This  disadvantage  is  at  least  partially  compensated  by  the  fact  that  in  the  case  of  an 
x-ray  mask  one  is  writing  on  a  thin  membrane  rather  than  a  bulk  substrate.  Most 
of  the  electrons  are  transmitted  through  the  membrane,  thus  greatly  reducing  the 
amount  of  back-scattered  electrons,  and  therefore  the  proximity  effect  compared  to 
bulk  substrates  [9].  Because  of  the  near  absence  of  a  proximity  effect,  it  is  possible 
to  achieve  linewidths  below  50  nm  in  high  density  patterns.  Figure  2-4  shows  such 
an  example,  in  this  case  the  line  is  about  30  nm  wide.  Linewidth  as  narrow  as  25  nm 
have  been  achieved  this  way. 

2.3.3  Daughter  Mask 

After  the  mother  mask  has  been  completed,  a  copy,  the  daughter  mask,  is  produced, 
using  x-ray  lithography.  This  replication  serves  two  purposes:  it  reverses  the  polarity 
and  produces  a  dark-field  mask  suitable  for  liftoff  in  device  processing;  and  it  reduces 
the  time  on  the  e-beam  lithography  tool.  If  a  mask  should  break,  a  new  daughter 
mask  can  be  made  within  a  short  amount  of  time. 

The  x-ray  exposure  for  mask  replication  is  done  with  a  microgap  between  mother 
and  daughter  mask.  The  size  of  the  maximum  practical  gap  is  determined  by  feature 
size  using  eq.  2.1,  with  a  set  equal  to  1  or  1.5.  The  masks  are  placed  on  top  of 
another,  with  the  gap  set  by  either  studs  evaporated  onto  the  mesa  or  by  the  natural 
nonfiatness  of  the  mesa.  In  both  cases,  the  gap  is  targeted  between  2  /zm  and  6  /zm. 
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Figure  2-4:  SEM  micrograph  of  a  device  written  directly  with  an  electron  beam  on  a 
mother  mask.  The  line  is  approximately  30  nm  wide. 

For  masks  with  ideally  flat  mesa  rims,  one  can  put  studs  on  the  rim  of  the  mother 
mask  to  establish  an  initial  gap.  The  gap  can  then  be  reduced  to  about  1  /urn  by 
partially  evacuating  the  space  between  the  masks  [38]. 

In  the  following  paragraphs,  the  calculation  of  gap  and  gap  capacitance  is  demon¬ 
strated  and  the  apparatus  for  replicating  x-ray  masks  at  a  controlled  gap  of  ^  1  fiva 
is  described  [38].  Such  a  gap  is  necessary  for  faithful  replication  of  mask  features 
^  40  nm. 

A  schematic  of  the  microgap  mask  replication  apparatus  is  shown  in  Fig.  2-5.  The 
capacitance  between  the  two  masks  is  sensed  using  a  Boonton  72B  capacitance  meter 
and  fed  back  into  a  circuit  that  adjusts  the  partial  pressure  between  the  masks  via 
an  MKS  248A  control  valve  to  keep  the  capacitance  constant.  The  pressure  between 
the  masks  reaches  a  constant  value  within  a  few  seconds. 

In  order  to  calibrate  the  apparatus,  gaps  were  measured  both  capacitively  and 
interferometrically.  The  latter  is  easily  accomplished  using  a  monochromatic  source 
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Figure  2-5;  Schematic  of  microgap  mask  replication  apparatus.  In  this  case,  the 
mother  mask  features  studs,  in  order  to  keep  the  masks  at  a  certain  gap  across  the 
mesa.  This  technique  is  necessary  for  optically  fiat  mesas,  while  a  nonflat  mesa 
provides  for  a  natural  gap  of  about  the  same  size.  A  Boonton  72B  bridge  was  used 
as  capacitance  sensor,  an  MKS  248A  as  control  valve. 
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while  viewing  the  fringes  between  the  two  mask  membranes  [35].  Before  inserting  the 
masks  into  the  apparatus,  any  large  dust  particles  should  be  removed.  This  can  be 
done  using  an  optical  flat,  as  described  earlier. 


Calculation  of  membrane  deflection 


The  gap  can  be  measured  and  monitored  during  the  exposure  setup  with  the  mo¬ 
nochromatic  green  light.  The  gap  can  also  be  calculated  using  elementary  theory  of 
elasticity.  In  the  case  of  two  circular  membranes,  one  can  use  polar  coordinates  to 
analytically  solve  for  the  bow  of  the  membranes.  The  calculation  for  a  similar  problem 
is  done  in  Ref.  [39].  If  the  pressure  difference  between  the  top  and  the  bottom  of  the 
membrane  is  denoted  as  P,  we  have  to  solve  the  following  equation 


Ijd  ^  _  _P 
r  dr  ^  dr  T  ’ 


(2.3) 


where  T  is  the  absolute  magnitude  of  the  stretching  force  per  unit  length  of  the  edge 
of  the  membrane,  C  is  the  deflection,  and  r  is  the  radius  variable  as  indicated  in  the 
inset  of  Fig.  2-6.  The  solution  for  the  deflection  of  the  membrane  is 

P-4) 

where  R  is  the  radius  of  the  membrane.  Note  that  only  the  amplitude,  but  not 
the  shape,  of  the  deflection  depends  on  parameters  like  pressure  and  forces  on  the 
membrane  edges.  The  deflection  can  also  be  measured  for  the  gap  between  the  two 
masks.  If  the  masks  have  identical  physical  properties,  the  deflection  of  a  single 
membrane  should  be  half  the  gap  between  the  masks.  The  gap  as  a  function  of  radius 
variable  r  can  be  found  by  determining  the  location  of  the  interference  fringes.  The 
measured  deflection  in  Fig.  2-6  is  in  good  agreement  with  the  fit  using  eq.  (2.4). 
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Figure  2-6:  Plot  of  normalized  deflection  as  a  function  of  normalized  radius.  Measured 
values  are  indicated  with  circles,  the  line  shows  the  fit  using  eq.  (2.4).  The  inset  shows 
a  schematic  of  the  system,  indicating  the  variables  used  in  the  calculations. 
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Calculation  of  capacitance  vs.  gap 


Given  eq.  (2.4),  one  can  calculate  the  capacitance  as  a  function  of  minimum  gap.  The 
capacitance  C  of  the  membrane  is  defined  as: 


2n  j-R  rdrdO 
Jo  Go  -  2C(r) 


where  c  is  the  dielectric  constant  and  Go  is  the  gap  without  any  force  applied,  which 
is  taken  to  be  uniform  on  the  circumference  of  the  membrane.  The  deflection  is 


independent  of  9,  so 


rR 

C'  =  27re 

Jo  Cjq  — 


UR^-r^) 


We  observe  that  the  maximum  deflection  of  each  membrane  is  related  to  the  minimum 


gap  Gmin  and  Go  as 


Cmax  —  C(^  0) 


Go  -  Gmin 


This  leads  to 


27r€jR^ 

Go  Gmin  •'0 


Gmin 


1  +  & 


Now  we  can  substitute  X  =  Go/ (Go— Gmm)—l+(^/-R)^j  so  that  dX  —  2{r/R)d{r /R). 
The  boundaries  of  the  integral  are  Go! {Go  —  Gmin)  —  1  and  Go/(Go  —  Gmin)i  and  we 


G'  = 


Go  Grain 


Gn 

'  Gq  Gmin 

‘  gfl _ 

Go  —  Gmin 


This  is  the  capacitance  due  to  the  deflected  membrane.  There  is  also  a  capacitance 
due  to  the  mesa  which  does  not  deflect.  The  total  capacitance  is  then 


'KeR?  .  Go  .  ‘^^{R-m  ~  R^) 

C  =  7^ — y; —  + - 7^ - 

^0  ^min  ^min  '^O 


(2.10) 
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Here  Rm  is  the  radius  of  the  mesa,  which  in  this  case  is  3  mm  larger  than  the  radius 
of  the  membrane.  The  capacitance  C  as  a  function  of  minimum  gap  Gmin  is  plotted 
in  Fig.  2-7  for  a  range  of  initial  gaps  Gq  between  2.5  and  6  /rm.  These  values  were 
calculated  using  eq.  (2.10)  and  are  plotted  as  full  lines.  The  minimum  gap  as  a 
function  of  capacitance  was  also  measured  by  determining  the  initial  gap  Go  and 
then  counting  fringes  for  different  values  of  the  capacitance.  Those  values  are  plotted 
as  circles  in  Fig.  2-7.  The  measured  values  are  offset  from  the  theoretically  predicted 
values  by  about  0.8  nF,  which  can  be  explained  by  an  uncertainty  in  the  calculation 
of  CjTiesa- 

Scanning  electron  micrographs  of  a  pattern  on  the  parent  x-ray  mask  made  by 
e-beam  lithography  and  its  replication  onto  a  daughter  mask  using  the  microgap 
mask  replication  apparatus  at  a  gap  of  ~1  ^m  are  shown  in  Fig.  2-8.  This  pattern 
is  used  to  define  the  gates  in  a  Coulomb-blockade  device  that  can  be  biased  to  yield 
two  quantum  dots,  coupled  to  each  other  via  a  narrow  tunneling  barrier,  or  as  one 
quantum  dot  with  three  leads  [40,41].  The  lithographic  barrier  width  is  less  than 
50  nm,  which  is  critical  for  the  functioning  of  the  device,  i.e.  for  efficient  tunneling. 
While  the  tunneling  barrier  -  the  line  dividing  the  two  dots  -  printed  without  loss  of 
resolution,  the  corners  appear  slightly  rounded.  Since  the  gold  is  plated  almost  to 
the  thickness  of  the  resist,  a  top  view  of  a  plated  mask  accentuates  any  rounding  of 
edges  on  the  top  of  the  photoresist. 


2.4  Development 

After  e-beam  exposure,  the  mother  masks  are  developed  in  a  mixture  of  methyl- 
iso-butyl-kytene  (MIBK)  and  iso-propyl-alcohol  (IPA).  Mother  masks  are  generally 
immersion  developed  because  it  is  easier  to  maintain  controlled  and  repeatable  condi¬ 
tions  in  a  beaker  than  by  puddle  development.  The  main  disadvantage  of  immersion 
development  is  that  the  fluid  is  not  constantly  being  replenished  and  particles  floating 
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Figure  2-7:  Plot  of  capacitance  versus  minimum  gap.  The  full  lines  correspond  to 
theoretical  capacitances  for  different  initial  gaps  between  3  and  6/im.  The  measured 
points  are  determined  by  measuring  the  initial  gap  and  then  counting  interference 
fringes.  The  dashed  line  is  the  capacitance  for  a  4  //m  initial  gap  with  an  additional 
offset  accounting  for  uncertainty  in  Cmesa-  The  dotted  line  gives  the  capacitance  of  a 
parallel  plate  capacitor  for  comparison.  The  images  on  the  bottom  are  photographs 
of  the  setup  at  different  values  of  minimum  gap. 
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Figure  2-8:  Scanning  electron  micrographs  of  parent  and  daughter  masks  of  a  pattern 
containing  40  nm  features.  The  gold  absorber  is  plated  to  a  thickness  of  200  nm. 
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on  the  surface  of  the  developer  can  attach  themselves  to  the  sample.  The  develop¬ 
ment  is  usually  done  in  a  2:1  IPA:MIBK  solution  at  a  temperature  of  21  °C  for  90 
seconds.  Upon  mixing  of  the  two  chemicals,  the  temperature  drops  to  about  16  °C. 
The  temperature  has  to  be  increased  through  heating  and  stabilized  at  the  desired 
point. 

Daughter  masks,  which  are  x-ray  exposed,  have  a  larger  exposure  latitude  than 
e-beam  exposed  mother  masks  and  are  therefore  generally  puddle  developed  with  a 
mixture  of  3:2  IPA:MIBK.  Puddle  development  means  that  a  squirt  bottle  is  used  to 
constantly  supply  fluid  to  a  puddle  on  the  sample,  thereby  constantly  replenishing  the 
fluid.  This  method  has  an  advantage  in  that  reduces  particle  contamination  problems. 

All  masks  are  subjected  to  a  de-scum  process  between  development  and  plating. 
PMMA  is  believed  to  leave  a  very  thin  (mono)  layer  of  scum  on  the  developed  surface 
which  would  cause  poor  electroplating  if  it  were  not  removed.  The  parameters  are 
given  in  Appendix  A. 

2.5  Electroplating 

After  development,  the  masks  are  electroplated  with  gold  absorbers  to  a  thickness  of 
200  nm,  which  corresponds  to  10  dB  attenuation  for  the  Cu-L  x-ray  sources  that  are 
used  in  our  laboratory  (A  =  1.32  nm). 

SEL-REX  BDT510,  a  commercially  available  gold  plating  bath,  is  used  for  elec¬ 
troplating  the  absorbers  on  the  x-ray  masks.  Although  the  plating  solution  and 
its  brightening  agent  are  proprietary,  some  information  can  be  gained  from  Material 
Safety  Data  Sheets.  The  gold  is  contained  in  the  solution  in  the  form  of  gold  sulphite. 
This  is  a  toxic  substance,  but  it  is  still  safer  to  work  with  than  many  competing  gold 
plating  solutions,  which  are  typically  cyanide-based.  If  the  bath  is  not  contaminated 
by  foreign  substances,  the  purity  of  the  deposited  gold  reaches  99.9%,  according  to 


the  manufacturer. 
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The  quality  of  the  plating  must  be  monitored  continuously,  and  quality  indicators 
such  as  plating  roughness  and  film  stress  must  be  correlated  to  parameters  such  as 
current  density,  reference  electrode  voltage,  bath  temperature,  brightener  concentra¬ 
tion,  and  agitation.  The  brightener  is  a  proprietary  additive  to  the  plating  bath  [42] , 
which  induces  grain  nucleation  and  therefore  tends  to  make  the  gold  fine-grained  and 
“bright”.  The  brightener  works  by  bonding  to  the  gold  grain,  terminating  its  growth, 
thus  inducing  nucleation  of  new  grains.  If  the  grains  are  small  compared  to  the  wave¬ 
length  of  visible  light,  there  is  no  diffuse  scattering  on  the  surface  of  the  gold,  making 
it  appear  “bright”.  Since  the  chemical  composition  of  the  brightener  is  not  known,  it 
is  important  to  check  the  characteristics  of  the  plating  bath  periodically. 

One  important  characteristic  for  a  plating  bath  is  the  current- voltage  (IV)  plot. 
The  IV-characteristics  of  BDT510  change  as  a  function  of  brightener  concentration,  as 
indicated  in  Fig.  2-9.  For  reference,  typical  plating  currents  used  in  our  laboratory  are 
between  10  and  20  mA  for  the  plating  geometry  used  in  these  experiments.  When  a 
new  bath  is  used,  such  an  IV-measurement  can  be  done  after  the  addition  of  brightener 
and  periodically  thereafter  to  keep  track  of  the  electrical  changes  that  occur  during 
the  lifetime  of  the  bath.  A  typical  change  of  IV-characteristic  with  time  for  a  bath 
used  in  our  laboratory  is  depicted  in  Fig.  2-10.  The  addition  of  brightener  significantly 
reduces  the  voltages  and  decreases  the  slope  for  the  currents  of  interest  in  the  linear 
portion  of  the  IV-curve.  Over  time,  the  voltage  increases  for  identical  current  density, 
and  is  larger  than  in  the  bath  without  brightener  after  a  couple  of  weeks.  It  is  obvious 
that  the  bath  changes  its  chemistry  during  its  lifetime,  possibly  due  to  the  oxidation 
or  evaporation  of  brightener  and  water.  After  addition  of  more  brightener,  the  plating 
voltage  is  again  reduced,  as  shown  in  Fig.  2-10.  The  characteristics  indicated  here 
are  for  a  bath  which  was  not  hermetically  sealed  when  not  in  use  and  thus  allowed 
for  significant  evaporation  and  surface  reactions  such  as  oxidation.  The  bath  is  now 
sealed  after  each  use  to  reduce  such  rapid  deterioration  of  bath  quality,  which  has 
caused  its  average  lifetime  to  apporoximately  double. 
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Figure  2-9:  Current-voltage-characteristic  of  a  plating  bath  as  a  function  of  bright- 
ener  concentration.  The  plot  shows  the  voltage  of  the  cathode  with  respect  to  a 
saturated  Calomel  electrode.  The  recommended  brightener  concentration  is  10  ml 
per  gallon.  The  cathode  area  in  this  experiment  was  about  25  cm^,  so  that  typical 
current  densities  used  during  plating  correspond  to  between  10  and  20  mA  in  this 
plot. 


2.5.1  The  Plating  Apparatus 

In  order  to  plate  an  x-ray  mask,  a  plating  apparatus  has  to  be  designed  and  built. 
Plating  in  our  laboratory  is  typically  done  in  one  large  bath  which  is  used  for  several 
months  and  which  is  strictly  monitored  during  its  lifetime  [29,24].  Although  efforts 
were  made  to  reduce  fluctuations  in  plating  quality  by  monitoring  the  pH  value,  con¬ 
ductivity,  stress  of  the  plated  film,  and  temperature,  plating  was  at  times  unreliable. 
Problems  were  usually  solved  by  changing  the  bath,  which  consists  of  about  5  liters 
of  plating  solution.  In  order  to  reduce  costs  and  to  have  more  flexibility  with  experi¬ 
ments,  a  new  fixture  was  designed  with  the  explicit  intent  of  reducing  the  amount  of 
plating  fluid  used.  A  schematic  of  the  fixture  is  shown  in  Fig.  2-11.  A  small  beaker. 
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Figure  2-10:  Current-voltage  characteristic  of  a  plating  bath  as  a  function  of  time. 
Also  shown  is  a  curve  for  a  bath  with  no  brightener.  The  addition  of  brightener 
reduces  voltages,  whereas  they  typically  increase  during  the  life  of  a  bath.  Typical 
currents,  corresponding  to  current  densities  used  in  the  fabrication  of  x-ray  masks, 
are  10  to  20  mA. 


containing  approximately  250  ml  of  plating  solution,  is  used.  The  plating  solution  is 
heated  on  a  hotplate  to  a  temperature  of  33  °C.  The  mask  is  then  placed,  membrane- 
side-up,  into  the  beaker,  followed  by  a  teflon  ring,  which  serves  as  a  spacer  between 
mask  and  platinum  anode.  Gold  coated  magnet  wire  on  the  bottom  of  the  teflon  ring 
ensures  electrical  contact  to  the  mask,  which  is  the  cathode  in  the  plating  setup.  The 
insulator  of  the  magnet  wire  is  removed  and  the  copper  is  gold  plated  for  contact  to 
the  mask.  The  insulator  of  the  magnet  wire  is  effective  in  the  solution.  A  platinum 
covered  mesh,  the  anode,  is  then  placed  on  top  of  the  teflon  ring. 

One  of  the  drawbacks  of  this  setup  is  that  the  solution  is  not  constantly  filtered,  as 
is  the  case  in  previous  setups.  In  my  experience,  this  has  not  led  to  mask  defects  due 
to  the  plating  process.  In  most  cases,  particles  get  introduced  during  fabrication  and 
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Figure  2-11:  Schematic  of  the  fixture  used  for  plating  of  x-ray  masks.  A  pyrex  beaker 
is  used  to  contain  the  plating  solution,  limiting  the  amount  of  plating  solution  to 
about  250  ml. 


handling  of  the  mask  blank  rather  than  subsequent  processing  and  blanks  can  easily 
be  screened  for  particles  before  use.  The  major  advantage  of  the  plating  setup  is  that 
the  plating  solution  can  be  changed  much  more  frequently,  such  as  once  a  month  or 
more,  without  excessive  costs.  This  setup  was  assembled  entirely  from  components 
available  in  the  laboratory.  A  more  sophisticated  setup,  incorporating  these  ideas, 
could  easily  be  built. 

A  power  supply,  which  allows  either  galvanostatic  or  potentiostatic  control,  was 
designed  to  be  used  with  the  experimental  plating  setup.  In  most  experiments,  the 
power  supply  was  used  as  a  galvanostat  in  which  the  voltage  is  regulated  to  impose  a 
constant  current.  The  circuit  diagram  of  this  supply  and  a  short  explanation  of  the 
circuit  is  given  in  Appendix  C. 

If  this  power  supply  is  used  for  plating,  the  value  of  the  area  to  be  plated  is  set 
at  the  plating  box.  This  enables  the  electronic  circuit  to  estimate  a  thickness  of  the 
plated  gold  assuming  that  the  plating  efficiency  is  unity.  The  current  is  integrated 
over  time  and  a  voltage  proportional  to  the  plating  current  and  the  plated  thickness 
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is  constantly  monitored  using  voltmeters.  It  is  also  a  good  idea  to  connect  a  chart 
recorder  to  the  plating  setup,  recording  simultaneously  the  plating  current  and  the 
voltage  of  the  saturated  Calomel  electrode  (SCE).  In  this  way,  it  is  easy  to  detect 
problems,  even  after  the  plating  run  is  done.  The  SCE  voltage  should  be  similar 
for  comparable  samples  and  equal  cathode  current  densities.  Any  deviations  suggest 
possible  problems  with  gold  quality  and  bath  deterioration. 

2.5.2  Determination  of  Plated  Thickness 

It  is  very  important  to  determine  the  thickness  of  the  plated  gold  with  great  accuracy. 
One  generally  tries  to  keep  the  resist  thickness  to  a  minimum  in  order  to  avoid  prob¬ 
lems  associated  with  large  aspect  ratio  resist  lines,  e.g.  falling-over  lines,  adhesion 
problems,  or  line  broadening  in  the  e-beam  lithography  due  to  forward  scattering  of 
electrons  in  the  thick  resist.  If  one  wants  to  plate  absorbers  to  a  thickness  of  200  nm, 
the  resist  thickness  used  for  the  e-beam  or  x-ray  exposure  is  about  220  nm,  leaving 
little  room  for  overplating.  It  is  therefore  crucial  to  end  the  plating  within  a  few  tens 
of  nanometers.  If  the  exact  area  of  the  plated  surface  and  the  plating  efficiency  is 
known,  it  is  possible  to  plate  to  within  these  tolerances  using  proper  timing.  The 
plating  efficiency  indicates  how  efficiently  the  ionic  current  in  the  solution  transfers 
metal  onto  the  sample. 

The  most  straightforward  way  to  measure  the  plated  thickness  is  to  use  a  pro- 
filometer,  a  stylus  measurement  system  such  as  an  Alpha-Step.  This  has  been  done  on 
membranes  and  gives  very  consistent  results.  The  main  disadvantage  of  this  method 
is  the  fact  that  the  force  exerted  by  the  stylus  is  not  very  well  controlled  and  can  in 
some  instances  destroy  the  membrane.  To  avoid  such  a  scenario,  a  resist  step  on  the 
mesa  of  the  x-ray  mask  can  be  measured  instead.  The  mesa,  however,  is  sometimes 
far  removed  from  the  actual  region  of  interest  and  there  are  instances  in  which  the 
plating  is  nonuniform  across  the  x-ray  mask.  This  is  especially  problematic  if  fine  fea¬ 
tures  plate  slower  than  large  features,  which  was  occasionally  observed  in  old  plating 
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bath’s.  In  these  cases,  the  voltages  observed  during  plating  tended  to  be  higher  than 
normal.  In  such  situations,  a  profilometer  measurement  on  the  mesa  might  indicate 
sufficient  absorber  thickness,  while  in  the  regions  of  interest,  towards  the  center  of 
the  membrane,  the  absorbers  are  too  thin. 

A  preferred  method  to  measure  the  thickness  of  gold  plated  onto  an  x-ray  mask 
is  to  use  a  microscope  and  measure  the  transmission  of  visible  light  through  the 
plated  region,  as  indicated  in  Fig.  2-12.  A  substage  condenser  lens  is  used  to  image 
a  rectangular  aperture  onto  the  membrane  at  the  position  of  the  absorber.  Gold  is 
partially  transmissive  in  the  visible  region,  with  a  local  maximum  in  the  green  part 
of  the  spectrum,  and  it  is  possible  to  relate  the  intensity  of  the  transmitted  green 
light  to  a  given  gold  absorber  thickness.  The  substage  condenser  lens  is  movable, 
so  that  the  absorber,  which  in  most  cases  consists  of  rectangles,  can  be  imaged  by 
adjusting  the  aperture  to  fit  the  desired  field  of  view.  The  transmission  of  the  green 
light  through  the  sample  can  then  be  compared  to  the  transmission  of  the  green 
light  through  a  known  thickness  of  gold.  A  prerequisite  is  a  similar  transmission 
through  the  SiN  membrane  on  sample  and  reference  membrane.  An  advantage  of  this 
method  is  that  the  lens  chosen  to  observe  the  thickness  through  transmission  can  be 
a  high  magnification,  high  numerical  aperture  lens,  enabling  the  observer  to  recognize 
variations  of  absorber  thickness  on  a  lateral  scale  of  a  few  microns.  Using  a  150  W 
Xenon  light  source,  it  is  possible  to  visually  detect  transmitted  green  light  for  gold 
thicknesses  up  to  about  220  nm  for  a  lOx,  0.2NA  lens,  and  up  to  just  below  200  nm 
for  a  lOOx,  0.9NA  lens. 

It  turns  out  that  the  sensitivity  of  the  human  eye,  especially  at  low  intensities, 
makes  thickness  determination  within  a  hundred  Angstroms  possible.  It  should  be 
noted,  however,  that  this  method  is  not  quite  safe  for  the  human  eye,  as  a  mask  break¬ 
age  would  result  in  a  high  light  intensity  incident  on  the  observer’s  eye.  A  drawback 
of  the  current  system  is  the  fact  that  comparison  is  possible  only  after  exchanging  the 
sample  mask  with  the  mask  containing  the  standard  gold  thicknesses.  This  could  be 
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Figure  2-12:  Schematic  of  method  used  to  determine  plated  gold  thickness.  A  rect¬ 
angular  aperture,  illuminated  with  bright  light,  is  focused  onto  the  backside  of  a  gold 
absorber.  The  intensity  of  the  transmitted  light  is  dependent  on  the  gold  thickness.  A 
high  magnification,  high  NA  lens  can  be  used  for  observation,  permitting  high  lateral 
resolution  in  determining  the  absorber  thickness. 
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improved  easily  by  introducing  glass  slides,  evaporated  with  standard  thicknesses  of 
1000,  500,  200  and  100  Angstroms  for  example,  plus  plating  base,  into  the  illumina¬ 
tion  path  of  the  substage  condenser.  In  this  way,  transparent  areas  of  the  mask  would 
be  dimmed  and  one  would  have  to  match  up  this  intensity  with  the  intensity  of  the 
absorber  region  without  the  additional  filter  present.  This  method  would  eliminate 
the  uncertainty  regarding  SiN  thickness  of  sample  and  reference  membrane  and  the 
accuracy  of  the  thickness  reading  could  be  further  increased. 


2.6  Device  Fabrication 

To  make  quantum  effect  devices,  we  start  off  with  substrates,  grown  using  molecular 
beam  epitaxy  (MBE).  The  layer  sequence  is  similar  to  layer  structures  used  to  build 
High  Electron  Mobility  Transistors  (HEMTs),  with  the  difference  that  the  two  di¬ 
mensional  electron  gas  (2DEG)  is  closer  to  the  surface.  For  a  comprehensive  review 
on  such  structures,  see  ref.  [43].  Since  constrictions  of  the  electron  path  are  achieved 
through  electrostatic  confinement,  the  electron  gas  needs  to  be  close  to  the  gate  elec¬ 
trode.  A  typical  layer  structure  used  in  the  experiments  is  depicted  in  Fig.  2-13.  The 
material  is  grown  by  Prof.  Michael  Melloch  of  Purdue  University.  Starting  with  a 
semi-insulating  wafer  of  GaAs,  a  1  //m-thick  layer  of  undoped  GaAs  is  grown.  On 
top  of  that,  two  layers  of  Alo.3Gao.7As  are  grown,  a  20nm-thick  undoped  layer  and  a 
30  nm-thick  silicon-doped  layer.  The  latter  provides  the  charge  carriers  that  accumu¬ 
late  in  the  2DEG  in  the  GaAs  close  to  the  AlGaAs  interface.  The  layer  structure  is 
terminated  by  a  5  nm  layer  of  n-doped  GaAs.  The  undoped  AlGaAs  layer  provides  a 
spatial  separation  between  the  ionized  impurities  and  the  charge  carriers,  preventing 
extensive  scattering  at  low  temperatures.  Typically,  the  thicker  the  spacer  layer,  the 
larger  the  mobility  and  the  smaller  the  charge  density  in  the  well.  Since  we  need  good 
control  over  the  potential  in  the  2DEG  and  a  short  distance  to  the  surface,  the  layer 
structures  used  in  our  experiments  represent  a  trade-off  between  high  mobility,  high 
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Figure  2-13:  Layer  structure  used  as  substrate  for  the  fabrication  of  quantum  effect 
devices. 

sheet  density,  and  thin  layer  structure.  Typically,  the  layers  used  in  our  experiments 
have  a  mobility  of  about  250,000  cm^/Vs  and  a  carrier  density  of  about  5x10^^  cm 

An  outline  of  the  fabrication  process  is  shown  in  Fig.  2-14.  Starting  with  the  layer 
structure,  photolithography  is  performed  using  an  optical  mask  which  defines  the  mesa 
structure.  The  active  and  the  ohmic-contact  regions  of  the  devices  are  defined  and 
device  isolation  is  achieved.  I  use  an  etch  consisting  of  1000:20:6  H20:NH40H:H202 
for  about  20  s,  which  results  in  a  mesa  of  about  60  nm. 

To  put  down  the  ohmic  contacts,  photolithography  is  performed  using  the  ohmic 
contact  mask,  followed  by  a  short  de-scum  in  a  UV-ozone  cleaning  station  and  the 
removal  of  oxide  from  the  sample  surface  using  a  dip  in  buffered  ammonium  hydroxide 
solution.  Then  the  ohmic  contact  material  is  deposited  by  electron  beam  evaporation. 
The  layer  sequence  for  this  evaporation  is  Ni(2):Au(4):Ge(22):Au(44)Ni(10):Au(30), 
where  thicknesses,  indicated  in  parenthesis,  are  in  nanometers.  The  contacts  are 
subsequently  annealed  in  a  rapid  thermal  annealing  system  for  15  sec  at  430°C. 
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The  last  step  involves  the  deposition  of  the  gates  which  define  the  geometry  of 
the  potential  barriers  of  the  quantum  effect  device.  PMMA  (950K,  4%)  is  spun  onto 
the  substrate  to  a  thickness  of  about  250  nm  and  baked  for  an  hour  at  180  °C.  The 
sample  is  then  aligned  and  exposed  following  the  procedure  developed  by  Ghanbari 
and  Chu  [9,36].  The  exposure  is  done  in  our  laboratory  using  x-ray  lithography.  After 
development,  about  100  nm  of  metal  are  deposited  and  lifted  off.  For  electron  beam 
evaporators,  10  nm  titanium,  followed  by  90  nm  gold,  are  evaporated.  For  thermal 
evaporators,  100  nm  of  aluminum  can  be  used.  There  is  some  indication  that  use  of 
electron  beam  evaporators  cause  damage  in  2DEGs  [44].  Metals  as  thick  as  120  nm  are 
routinely  lifted  off  using  a  resist  thickness  of  about  250  nm.  The  complete  fabrication 
sequence  is  listed  in  more  detail  in  Appendix  D. 

2.6.1  Advantages  of  X-ray  Lithography  in  Fabrication 

There  are  only  a  few  lithography  technologies  that  can  achieve  the  resolution  sufficient 
for  patterning  the  gate  structures  of  quantum  effect  devices.  Optical  lithography  has 
made  great  progress  in  achieving  small  feature  sizes  but  it  is  still  far  from  the  necessary 
linewidth.  Electron-beam  lithography  is  the  tool  most  often  used  in  the  fabrication 
of  quantum-effect  devices.  There  is,  however,  some  indication  that  under  certain 
conditions  damage  is  inflicted  on  devices  during  exposure  with  an  electron-beam  [45]. 
In  the  case  of  ion-beam  lithography,  the  problem  of  damage  is  even  more  severe,  as 
indicated  by  Ismail  [4].  The  only  lithography  technique  that  can  directly  compete 
with  e-beam  lithography  is  x-ray  lithography  which  has  been  shown  to  be  capable 
of  resolving  features  down  to  lO’s  of  nanometers  [1],  and  has  further  been  shown  to 
cause  no  damage  in  modulation  doped  structures  [45]. 

X-ray  lithography  has  some  unique  advantages  over  e-beam  lithography.  Among 
them  is  the  fact  that  x-rays  can  penetrate  very  thick  resist  on  nearly  arbitrary  to¬ 
pography  without  scattering  [1].  Typically,  the  resist  thicknesses  exposed  in  our 
laboratory  are  300  nm,  but  this  is  not  a  fundamental  limit.  The  ability  to  expose 
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large-aspect-ratio  resist  structures  can  be  useful  if  a  thick  film  of  metal  is  lifted  off.  It 
is  not  necessary  to  add  an  additional  gate  metalization  to  ensure  good  step  coverage 
or  contact  pad  thickness  because  x-ray  lithography  provides  the  opportunity  to  lift 
off  such  thick  layers.  Consequently,  one  has  to  worry  less  about  static  charges  that 
could  destroy  thin  gates  or  about  step  coverage  if  one  wants  to  lead  thin  wires  onto 
a  mesa.  This  last  point  is  a  problem  in  making  planar  resonant  tunneling  devices  by 
e-beam  lithography  [46].  As  can  be  seen  in  Fig.  2-15,  step  coverage  of  100  nm  lines 
can  easily  be  achieved  for  a  mesa  height  of  about  120  nm. 


■H  K 

200  nm 


Figure  2-15:  Step  coverage  for  a  two  barrier  planar  resonant  tunneling  device.  The 
gate  fingers  are  separated  by  200  nm  and  the  metal  thickness  is  30  nm.  The  metal 
thickness  could  be  made  thicker.  Films  of  more  than  100  nm  are  frequently  used  in 
the  fabrication  of  fine  features.  This  improves  the  likelihood  of  good  continuity  across 
the  mesa  step  and  reduces  the  sensitivity  of  the  device  to  static  charges. 

It  was  already  mentioned  that  one  can  commit  extensive  time  and  effort  into  writ¬ 
ing  an  x-ray  mask  using  e-beam  lithography  because  subsequent  processing  steps  are 
parallel.  This  advantage  has  been  used  in  the  design  and  fabrication  of  optoelectronic 
devices,  such  as  channel-dropping  filters  [16].  A  “segmented  fiducial  grid”,  a  unique 


44 


CHAPTER  2.  X-RAY  LITHOGRAPHY  AND  DEVICE  EABRICATION 


tool  for  obtaining  long  range  coherence  with  an  e-beam  lithography  tool  [47],  was 
used  to  write  gratings  onto  an  x-ray  mask.  Optoelectronic  devices,  produced  using 
this  technique  have  significant  advantages  concerning  long  range  coherence  compared 
to  devices  made  using  other  techniques. 

Damage  during  processing 

In  order  to  see  quantum  effects,  it  is  crucial  not  only  to  use  high  quality  substrates 
but  also  to  maintain  the  high  quality  during  processing.  One  of  the  best  indicators 
for  substrate  quality  is  the  low-temperature  mobility.  The  higher  the  mobility  in 
a  particular  sample,  the  higher  the  temperature  at  which  quantum  effects  may  be 
observed.  It  is  clear  that  mobility  degradation  by  exposure  to  high  energy  particles 
during  device  fabrication  must  be  avoided. 

Damage  caused  by  electron  beam  evaporation  of  gate  metals  on  HEMT  struc¬ 
tures  [44]  and  on  bulk  materials  [48-50]  has  been  widely  reported.  The  most  impor¬ 
tant  sources  of  radiation  inside  an  electron  beam  evaporator  are  electrons  and  x-rays. 
Damage  caused  by  x-ray  emission  during  evaporation  has  been  previously  suggested 
as  a  possible  mechanism  in  silicon  systems  [51].  In  other  experiments  on  GaAs  [48], 
it  was  shown  that  through  proper  shielding  techniques,  electron  irradiation  could  be 
minimized.  Any  damage  was  then  attributed  to  stray  electrons  back-scattered  from 
the  evaporation  target.  Degradation  of  2DEG  mobility  by  e-beam  lithography  in  the 
voltage  range  of  5-12.5  keV  had  previously  been  observed  [52].  However,  at  higher 
energies  (up  to  20  keV),  the  authors  did  not  observe  any  degradation. 

We  observed  significant  mobility  degradation  in  modulation-doped  AlGaAs/GaAs 
Hall  bars  after  irradiation  by  50  keV  electrons.  Unirradiated  samples  on  the  same 
header  showed  no  degradation.  Similarly,  no  degradation  was  observed  when  devices 
were  exposed  to  Cu-L  x-rays  in  a  helium  atmosphere  [45].  These  results  are  consistent 
with  reports  of  mobility  degradation  due  to  metallization  of  GaAs  in  some  e-beam 
evaporators. 
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We  did  not  study  the  effect  of  e-beam  evaporation  on  device  quality  in  the  fab¬ 
rication  of  our  devices.  It  would  be  difficult  to  replace  the  e-beam  evaporator  with 
a  thermal  evaporator  for  the  ohmic  contact  evaporation  because  of  the  multilayer 
structure  used  in  the  ohmic  contacts.  In  addition,  the  damage  during  this  evapora¬ 
tion  should  not  cause  significant  problems  because  of  the  annealing  step  that  follows 
the  evaporation.  It  would  be  of  interest,  however,  to  evaporate  the  gate  metal  with 
a  thermal  and  an  e-beam  evaporator  and  compare  device  characteristics.  Instead  of 
low  temperature  mobility,  one  could  also  measure  ballistic  transport  through  quan¬ 
tum  point  contacts,  which  could  give  an  indication  of  scattering  centers  introduced 
during  processing. 
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Chapter  3 


Device  Operation  and 
Measurements 


3.1  Introduction 


Reduced  dimensional  systems  are  of  great  importance  in  physics  and  engineering. 
Most  of  mesoscopic  physics  is  based  on  such  systems.  The  physics  community  has 
ventured  into  the  two-,  one-  and  zero-dimensional  regimes  and  explored  exciting  new 
realms  of  physics.  The  engineering  community,  however,  has  not  yet  moved  below  two- 
dimensional  systems  in  large-scale  applications,  of  which  the  most  important  one  is 
the  metal  oxide  semiconductor  field  effect  transistor  (MOSFET).  But  with  increasing 
integration  density  and  reduced  device  dimensions,  there  is  an  end  in  sight  to  this 
route.  Increased  device  density  and  clock  rate  on  a  chip  already  lead  to  problems 
concerning  heat  dissipation  on  a  chip.  The  cooling  of  modern  CPUs  is  one  of  the 
most  important  issues  to  be  addressed  in  the  design  of  a  circuit  board.  It  is  likely 
that  this  problem  will  set  a  fundamental  limit  on  both  the  functionality  and  the  speed 
of  computers.  Novel  devices  might  offer  a  solution  to  these  problems,  necessitating 
continued  research  in  areas  such  as  quantum-effect  devices. 
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3.1.1  Reduced  Dimensional  Systems 

Two  dimensions 

In  a  two-dimensional  electron  gas  (2DEG),  electrons  are  constrained  to  move  in  only 
two  directions  by  a  potential  well  which  prevents  them  from  moving  in  the  third. 
A  2DEG  can  be  created  at  a  silicon-silicon  dioxide  interface,  as  is  the  case  for  the 
MOSFET,  or  at  the  interface  between  materials  of  different  bandgap,  as  is  the  case 
for  modulation-doped  field-effect  transistors  (MODFETs).  By  far  the  most  popular 
MODFET  is  based  on  the  GaAs-AlGaAs  material  system.  GaAs  and  AlGaAs  are 
closely  matched  in  their  lattice  constant,  and  one  material  can  easily  be  grown  on 
the  other  without  creating  large  dislocation  densities.  The  MODFET  structure  is  the 
basis  for  the  devices  that  were  used  in  the  measurements. 

The  electron  density  in  the  two  dimensional  electron  gas  can  be  varied  by  applying 
a  voltage  to  a  gate  located  on  top  of  the  layer  structure.  The  energy  E  of  the 
conduction  electrons  as  a  function  of  the  two  dimensional  wave  vector  variable  q  in 
a  2DEG  is 

fe2^2 

^(0  =  ,  (3.1) 

where  m*  is  the  effective  mass  of  the  electron  which  is  smaller  than  the  electron  mass, 
rrie,  due  to  interaction  with  the  lattice.  For  GaAs,  we  have  m*  ~  O.OGTmg.  The  wave 
vector  variable,  q,  lies  in  the  x-y  plane  if  z  is  the  coordinate  of  confinement,  and  its 
magnitude  is  |q|  =  In  such  a  system,  the  two-dimensional  density  of  states, 

n,  in  the  absence  of  a  magnetic  field  is 

777* 

mn{E)  =  ^  .  (3.2) 


This  means  that  the  density  of  states  is  independent  of  energy  up  to  the  energy  of  the 
first  excited  state  in  the  ^-direction.  The  Fermi  energy  of  a  2DEG  can  be  calculated 


3.1.  INTRODUCTION 


49 


easily  if  the  sheet  electron  density  ns  is  known  as 

^ 

n2D  2m* 


The  Fermi  wavevector,  /cf,  is  then 


(3.3) 


kp  =  'j2'Kns  ■  (3-4) 

At  low  temperatures,  phonon  scattering  is  negligible,  and  electrons  scatter  only  due 
to  ionized  impurities.  An  important  figure  of  merit  for  a  given  layer  structure  is  the 
elastic  scattering  length,  the  distance  that  an  electron  can  travel  on  average  without 
scattering  elastically  with  an  impurity.  To  estimate  the  elastic  scattering  length,  4) 
of  electrons  in  the  2DEG  one  uses  the  Drude  formula  (with  the  lifetime  r  =  le/vp) 
and  relates  the  elastic  scattering  length  to  the  Fermi  wavevector  and  the  mobility  of 
the  substrate, 

^  ^  ^  .  (3.5) 

e  e 

For  the  layer  structure  depicted  in  Fig.  2-13,  a  mobility  of  about  250,000  cm^/Vs  and 
an  electron  density  of  about  5  x  10^^  cm'^  results  in  an  elastic  scattering  length  of 
about  3/im  and  a  Fermi  wavelength  of  about  35  nm.  Layer  structures  with  mobilities 
above  lO^'cm^/Vs  have  been  reported  [53],  leading  to  elastic  scattering  lengths  in  the 
tens  of  microns.  The  structure  used  in  the  experiments  reported  here  represents  a 
trade-off  between  high  mobility  and  control  over  the  2DEG  due  to  its  close  proximity 
to  the  surface. 


One  dimension 

In  a  one-dimensional  conductor,  an  additional  potential  barrier  confines  the  electrons 
to  move  in  only  one  direction.  For  example,  if  the  movement  in  the  a;-direction  is 
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confined  by  a  potential  V{x),  the  energy  is  given  as 

^2  2 

EjiQy)  =  Ej{k^)  +  ^  .  (3.6) 

The  wavevector  Qx  =  kx  is  quantized  and  only  Qy  remains  as  a  continuous  vari¬ 
able.  If  the  potential  well  is  parabolic,  which  is  a  reasonably  good  approximation 
(V’(a:)  =  m*UQx'^/2),  the  energy  levels  are  equidistant  with  Ej  =  {j  —  ^)hujo.  The 
index  j  is  called  the  subband  index  and  indicates  how  many  conducting  channels 
are  open  in  the  confined  region.  If  only  states  in  the  band  j  =  1  are  occupied,  one 
refers  to  one-dimensional  conduction,  and  for  j  >  1  to  quasi-one  dimensional  conduc¬ 
tion.  One-dimensional  conduction  in  devices  has  been  dealt  with  extensively  in  the 
literature  [54-57].  For  a  review,  see  refs.  [58,59]. 

A  one-dimensional  conductor  can  be  fabricated  using  a  2DEG  by  confining  it  with 
two  metal  gates  on  top  of  the  layer  structure,  leaving  a  narrow  gap  of  less  than  0.5  //m 
between  them.  Such  a  device  is  called  a  quantum  point  contact  or  split  gate  device. 
By  applying  a  negative  voltage  to  the  gates,  the  2DEG  is  depleted  underneath  and 
only  a  narrow  channel  between  the  gates  is  left  open.  The  optimum  distance  between 
the  gates  to  observe  subband  quantization  depends  on  the  layer  structure  and  depth  of 
the  2DEG  below  the  surface  [60].  The  length  over  which  one-dimensional  conduction 
can  be  observed  depends  on  the  elastic  mean  free  path  of  the  electron,  which  in  turn 
depends  on  the  layer  structure  and  processing  steps.  To  observe  ballistic  transport 
in  a  quantum  point  contact,  the  length  of  the  one-dimensional  region  has  to  be  small 
compared  to  the  elastic  mean  free  path. 

The  conductance  of  a  ballistic  wire  at  T=0  can  be  calculated  in  the  following  way: 
a  voltage  drop  V,  applied  across  the  constriction,  causes  a  difference  in  the  chemical 
potential  /x  across  the  quantum  point  contact. 


eV  =  /xi  -  /X2  =  A/x  . 


(3.7) 


3.1.  INTRODUCTION 


51 


This  difference  in  chemical  potential  leads  to  a  current  due  to  the  uncompensated 
electron  states  on  one  side  of  the  quantum  point  contact.  The  current  I  can  be 
calculated  as 

/  =  <=!;  .  (3.8) 

j,k=l  ^ 

Here,  N  is  the  number  of  occupied  subbands,  ni£)j{E)  is  the  one-dimensional  den¬ 
sity  of  states  in  the  jth  subband,  Vj{E)  is  the  electron  velocity,  and  Tj^k{E)  is  the 
transmission  coefficient  for  an  electron  from  the  kth  into  the  jth  subband.  This 
formula  can  be  simplified  dramatically  since  the  one-dimensional  density  of  states 
riiD  j  =  2{dEj  jdq)~'^/'K  is  inversely  proportional  to  the  velocity  of  the  electron  Vj  (E)  = 
1/h  dEj/dq.  For  small  applied  voltages,  the  conductance  G  reduces  to 

^  E  •  (3-9) 

“  j,k=l 

This  is  the  Landauer  Formula  for  the  two  terminal  resistance  of  a  quantum  wire  [GI¬ 
GS].  In  the  case  of  ballistic  wires  where  no  scattering  takes  place  Tj^k  becomes  a  unity 
matrix,  and  the  conductance  is 

G  =  ^-N  .  (3.10) 

h 

The  conductance  in  a  ballistic  wire  depends  only  on  the  number  of  occupied  subbands 
and  not  on  the  geometry,  e.g.  length,  as  long  as  the  wire  remains  ballistic.  The 
conductance  is  quantized  and  the  step  height  is 

^  =  77.48//5  =  (1290Gf2)-^  .  (3.11) 

h 

If  one  biases  the  quantum  point  contact  such  that  the  chemical  potential  lies  below 
the  first  subband,  no  conduction  is  possible,  except  for  tunneling.  If  two  of  these 
barriers  are  biased  to  produce  an  island  of  electrons  separated  from  the  leads,  a  zero¬ 
dimensional  electron  gas  results,  as  will  be  discussed  in  the  following  section.  This 
arrangement  is  also  frequently  called  a  quantum  dot.  Coulomb  blockade  device,  or 
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single  electron  transistor. 

Zero  dimensions 

In  1988,  John  Scott-Thomas,  then  a  graduate  student  at  MIT,  fabricated  devices 
featuring  long,  one-dimensional  conductors  on  silicon  using  x-ray  lithography  [64], 
depicted  in  Fig.  3-1.  These  devices  were  not  expected  to  exhibit  ballistic  conduction 


Vq  (Volts) 


Figure  3-1:  Devices  in  which  the  discovery  of  Coulomb  blockade  in  semiconductor 
devices  was  made.  The  bottom  gate  was  used  to  define  a  long,  diffusive  wire  in  the 
silicon,  while  the  top  gate  was  used  to  sweep  the  chemical  potential  within  the  wire. 
The  graph  on  the  bottom  shows  the  conductance  of  the  wire.  After  ref.  [65]. 

characteristics,  since  the  mobility  in  the  silicon  was  not  very  high.  He  was  studying  the 
transport  properties  of  long  wires  in  the  diffusive  regime,  where  the  elastic  scattering 
length  is  longer  than  the  width  but  much  shorter  than  the  length  of  the  wire. 

He  measured  the  conductance  of  these  wires  in  a  dilution  refrigerator  and  expected 
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to  see  evidence  of  conductance  peaks  associated  with  the  population  of  subbands  in 
the  wire.  In  some  of  the  wires,  he  instead  observed  conductance  oscillations  which 
were  periodic  as  a  function  of  the  top  gate  voltage,  which  was  used  to  sweep  the 
chemical  potential  within  the  wire  [65].  According  to  our  current  understanding, 
in  each  of  the  wires  that  showed  these  oscillations  there  were  two  impurities  that 
produced  tunneling  barriers  within  the  wires  by  pinching  off  the  first  subband  [66]. 
In  effect,  an  electron  island  was  produced  and  Scott-Thomas  observed  the  filling  of 
the  island  with  electrons,  with  each  conductance  peak  indicating  the  addition  of  one 
electron  to  the  island.  This  effect  had  been  predicted  and  explained  by  Glazman  and 
Shekhter  [67],  who  applied  knowledge  of  the  Coulomb  blockade  in  superconducting 
tunneling  junctions  [68]  to  zero-dimensional  semiconductor  structures. 

An  intuitive  picture  for  the  explanation  of  conductance  peaks  as  a  function  of 
gate  voltage  is  presented  in  Fig.  3-2.  The  schematic  shows  a  conducting,  quasi-one- 
dimensional  channel  that  is  interrupted  twice,  isolating  a  dot  which  can  be  accessed 
through  tunnel  barriers  to  two  leads.  The  conductance  of  each  tunneling  barrier  has 
to  be  adjusted  below  the  conductance  quantum  of  2e^/h.  If  both  tunneling  barriers 
are  biased  simultaneously  at  these  values,  periodic  conductance  oscillations  can  be 
observed  by  sweeping  the  chemical  potential  within  the  dot  using  the  gate.  The  point 
of  Fig.  3-2  is  that  the  charge  in  the  dot,  rather  than  being  a  continuum,  is  quantized 
to  an  integer  number  of  electrons.  Thus  one  has  a  staircase  with  step  length  e/Cg, 
rather  than  the  classical  line  with  slope  Cg  in  the  Q-V  plot  of  Fig.  3-2.  At  each  step, 
electrons  can  enter  and  leave  the  dot,  resulting  in  current  peaks. 

Using  artificial  structures  that  were  designed  to  have  adjustable  tunneling  barriers, 
the  interpretation  of  the  observed  conductance  oscillations  as  single  electron  tunnel¬ 
ing,  proposed  by  van  Houten  [66]  and  Glazman  [70],  was  soon  confirmed  [71-73].  One 
of  the  first  such  measurement,  made  by  U.  Meirav,  is  shown  in  Fig.  3-3.  An  inverted 
heterostructure  was  used  to  sweep  the  chemical  potential  in  quantum  dots  of  various 
sizes.  The  conductance  oscillations  as  a  function  of  gate  voltage  depend  on  the  spatial 


INTRODUCTION 


55 


3.1. 


Figure  3-3:  Schematic  of  device  and  conductance  curves  for  quantum  dots  of  differ¬ 
ent  spatial  extent.  The  device  was  made  on  an  inverted  heterostructure,  allowing 
modulation  of  the  chemical  potential  from  the  backside  and  control  of  the  tunneling 
barriers  using  the  metal  gates.  The  period  of  conductance  oscillations  depends  on  the 
distance  between  tunneling  barriers.  After  ref.  [71]. 


extent  of  the  quantum  dot.  There  are  several  interesting  review  articles  [74-77]  and 
theses  [69,73,78]  on  this  subject. 

To  observe  the  Coulomb  blockade  effect  and  the  associated  conductance  oscilla¬ 
tions,  the  thermal  energy  should  be  small  compared  to  the  charging  energy, 

(3.12) 

where  Cs  is  the  sum  of  all  capacitances  to  the  dot,  ks  is  Boltzmann’s  constant  and  T 
the  absolute  temperature.  Due  to  the  capacitances  involved  in  semiconductor  struc¬ 
tures  and  the  limitations  imposed  by  eq.  (3.12),  most  measurements  of  this  Coulomb 
blockade  effect  are  done  at  cryogenic  temperatures,  often  in  the  mK-range.  There  is, 
however,  an  interest  to  apply  these  new  effects  in  engineering.  The  phenomenon  that 
only  an  integer  number  of  electrons  can  tunnel  into  the  quantum  dot  has  been  used  in 
so-called  turnstile  devices  [79]  to  produce  quantized  current  plateaus  I  =  ef,  where 
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f  is  the  frequency  by  which  the  tunneling  barriers  are  modulated.  This  device  could 
be  used  as  a  current  standard.  Other  authors  have  proposed  logic  families  [68,80]  and 
storage  circuits  [81]  using  the  Coulomb  blockade  effect.  Practical  applications  have 
not  been  implemented  yet  since  low  temperatures  are  required  to  observe  the  desired 
characteristics.  Researchers  are  trying  to  increase  the  temperatures  at  which  Coulomb 
blockade  can  be  observed  by  reducing  capacitances  of  the  quantum  dot  to  its  environ¬ 
ment  [81,82].  One  of  the  most  exciting  recent  developments  is  the  apparent  discovery 
of  Coulomb  blockade  in  short-channel  silicon  on  insulator  (SOI)  MOSFETs  [83].  The 
authors  used  the  pattern-dependent  oxidation  of  silicon  to  produce  the  tunneling  bar¬ 
riers.  One  of  the  advantages  of  SOI  technology  is  the  reduction  of  capacitances,  which 
also  happens  to  be  the  requirement  for  the  operation  of  a  high  temperature  quantum 
dot.  In  these  devices,  the  total  capacitance  of  the  dot  was  about  2aF,  and  strong 
Coulomb  blockade  oscillations  as  a  function  of  gate  voltage  were  observed  at  77  K  and 
conductance  modulations  remained  even  at  room  temperature. 

This  chapter  continues  with  a  derivation  of  the  conductance  condition  for  the 
single  quantum  dot  in  Sec.  3.2.  The  derivation  is  kept  general  so  it  can  be  used 
for  comparison  during  the  derivation  for  the  double  dot  in  Sec.  3.5,  which  is  done 
in  close  analogy.  A  short  description  of  the  measurement  setup  and  the  probe  used 
for  low  temperature  measurements  are  given  in  Sec.  3.3.  Section  3.4  introduces  the 
device  that  was  measured  in  the  experiments  and  discusses  the  results  of  the  single 
quantum  dot  measurements.  After  discussing  the  physics  of  the  double  dot  in  Sec.  3.5, 
measurements  are  shown  and  discussed  in  Sec.  3.6. 


3.2  The  Quantum  Dot 

Conduction  through  the  quantum  dot  can  occur  only  if  the  probability  of  having  N 
electrons  in  the  dot  is  just  as  high  as  having  A’— 1  electrons  in  the  dot.  The  probability 
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of  having  N  electrons  in  the  dot  is  given  by  the  Grand  Canonical  Ensemble  [76,84] 


P{N)  =  exp(— 


keT  ’ 


(3.13) 


where  E  is  the  grand  partition  function  (a  constant  in  this  case),  fj,  is  the  chemical 
potential,  F  is  the  free  energy,  ks  is  the  Boltzmann  constant  and  T  is  the  temperature. 
At  very  low  temperature,  the  probability  P{N)  is  negligible  for  all  but  one  value  of 
N,  the  one  that  minimizes  {F{N)  —  N/i).  Conductance  is  possible  only  if  P{N)  and 
P{N  —  1)  are  both  nonzero  for  some  N.  This  is  the  necessary  condition  for  tunneling 
to  take  place.  If  this  is  the  case,  a  current  flows  through  the  dot,  and  the  number  of 
electrons  alternates  between  N  and  N  —  1.  This  implies  that 


F{N)-F{N-l)  =  n{N)  .  (3.14) 

At  zero  temperature,  the  free  energy  equals  the  ground-state  energy  of  the  dot,  which 
is  the  sum  of  the  charging  energy  U  and  a  contribution  due  to  the  energy  level 
quantization  within  the  dot, 

F{N)  =  U{N)  +  ^Ep  .  (3.15) 

p=i 

For  the  moment,  we  neglect  the  energy  level  quantization  in  the  dot,  which  is  not 
a  valid  assumption  if  we  are  dealing  with  extremely  small  quantum  dots  of  lateral 
dimensions  of  the  order  of  the  Fermi  wavelength.  This  case  was  treated  in  detail 
by  Beenakker  [85]  and  will  be  left  aside  for  the  moment.  Its  implications  will  be 
discussed  in  Sec.  3.2.1.  Using  eqs.  (3.14)  and  (3.15),  the  chemical  potential  is  deflned 
as  the  energy  cost  of  adding  the  Nth  electron  to  a  dot  with  {N  —  1)  electrons. 


ti{N)  =  U{N)-UiN  -1)  . 


(3.16) 
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For  a  change  in  charging  energy  6U  as  a  function  of  charge  in  the  dot  Qd  and  the 
electrodes  Qi,  we  can  write 


6U  =  Vd5Qd  +  'ZVi5Qi  ,  (3.17) 

i 

where  Vi  are  the  potentials  of  the  electrodes,  and  Vd  is  the  potential  in  the  dot. 
Unfortunately,  we  have  control  only  over  the  voltages  Vi  and  not  over  the  charges 
Qi.  It  is  desirable  to  interchange  the  dependence  of  the  quantities  in  the  conjugate 
pair.  To  accomplish  this,  we  use  the  Legendre  transformation  to  define  a  co-charging 
energy  U  so  that  U  =  U  —  Y,i  ViQi  and  we  get 


SU  =  VdSQd-^QiSVi  . 

i 

It  is  a  straightforward  exercise  to  transform  U  to  U  since  all  the  essential  information 
is  preserved  during  the  Legendre  transformation.  For  now  we  will  rename  U  as  the 
charging  energy  U.  To  calculate  the  magnitude  of  the  charging  energy,  we  integrate 
over  all  charges  in  the  dot  with  all  electrode  voltages  constant  {5Vi  =  0), 

fQd 

U(Qi)  =  /  V((3i,  m)6Q'j  .  (3.18) 

Jo 

The  total  charge  in  the  dot,  Qd,  is  related  to  the  potentials  on  the  electrodes  Vi  by 
the  capacitance  matrix  elements  c* 


Qd  =  J2^'iiVi  .  (3.19) 

i 

The  capacitance  matrix  elements  Cdi  are  called  self  capacitances  for  d  =  i  and  mutual 
capacitances  for  d  ^  i.  The  connection  to  circuit  parameter  capacitances  will  be 
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treated  later  on.  The  potential  within  the  dot  is  then  given  as 


Vd  =  —[Qd-J2caiVi\  = 

^dd  i^d 


Qd  ('di  .y 

Cdd  j  Cdd 


(3.20) 


The  self  capacitance  Cdd,  according  to  electrostatic  theory,  turns  out  to  be  identical  to 
the  sum  off  all  capacitances  to  the  dot,  introduced  as  Cs  in  eq.  3.12.  This  will  be  dis¬ 
cussed  later.  The  total  charging  energy  can  be  calculated  by  adding  up  contributions 
of  all  charges  in  the  dot  using  eq.  (3.18), 

U{Qi)  =  -  Y.  —VMQ'i  ■  (3-21) 

Jo  Cdd  j  ('dd  2^dd  i  ^dd 

If  we  sweep  only  one  electrode  potential,  Vg,  the  charging  energy  U  increases  linearly 
with  this  potential.  One  can  plot  the  energy,  as  is  done  in  Fig.  3-4,  for  several  values 
of  N  and  recognize  that  the  charging  energy  of  the  dot  can  be  determined  easily  under 
the  condition  that  the  charging  energy  is  minimized  for  each  value  of  gate  voltage. 
We  can  pick  out  the  lines  that  represent  N  =  3  and  =  4  electrons,  as  shown  in 


Figure  3-4:  Charging  Energy  of  the  dot  as  a  function  of  gate  potential  for  Cdg  =  10  aF. 
The  system  is  in  the  state  that  minimizes  the  charging  energy.  After  ref.  [75]. 
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Fig.  3-5,  and  see  that  at  the  point  of  intersection  it  is  energetically  just  as  favorable 
to  have  3  or  4  electrons  in  the  dot.  At  this  value  of  the  gate  potential,  we  expect  to 


Figure  3-5:  Close-up  of  Fig.  3-4  for  =  4  electrons. 

see  a  conductance  peak,  since  the  dot  can  oscillate  between  3  and  4  electrons,  and 
current  can  flow.  If  the  gate  potential  is  between  points  of  intersection,  the  current 
is  suppressed  by  the  Coulomb  blockade. 

Using  eqs.  (3.16)  and  (3.21),  we  can  calculate  the  chemical  potential  as 


KN) 


U{-eN)-U{-e{N  -1)) 
-2  1^ 

2' 


^di 


— (JV-rj  +  EeV; 

Cdd  2  .  Cdd 


(3.22) 


This  is  the  chemical  potential  as  a  function  of  electrode  voltages  for  a  single  dot.  The 
chemical  potential  is  swept  as  a  function  of  gate  voltage  Vg  by  a  factor  ecdg/cdd..  The 
ratio  of  gate  capacitance  to  self  capacitance  is  called  the  capacitive  lever-arm.  It  is 
a  measure  of  the  amount  by  which  the  chemical  potential  in  the  dot  is  swept  for  a 
change  of  the  potential  of  a  gate. 

To  calculate  the  linear-response  conductance  G  of  the  quantum  dot,  we  use  G  = 
I /V  for  the  limit  of  U  — >  0.  In  this  case,  source  and  drain  are  at  zero  potential,  and 
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fi  =  0  is  the  condition  for  current  flow.  A  typical  measurement  of  the  quantum  dot 
is  made  by  sweeping  only  one  of  the  electrodes,  i  =  g,  with  all  other  electrodes  kept 
constant.  In  this  case,  the  gate  voltages  at  which  conductance  peaks  occur  are 


^dg  i:^g  ^  Cdg 


(3.23) 


This  means  that  the  conductance  maxima,  defined  as  an  increment  of  N  by  one,  are 
spaced  by 


(3.24) 


So  far,  we  have  only  dealt  with  capacitances.  To  make  the  result  more  intuitive, 
we  transform  the  capacitances  into  capacitors. 


Circuit  Parameter  Representation:  Let  each  electrode  j  be  coupled  to  a  charge 
Qi  by  a  capacitor  Cij,  and  let  all  charges  be  imaged  within  the  system.  The  charge 
Qi  on  conductor  i  is  known  when  all  the  voltages  in  the  system  Vj  and  capacitance 
matrix  elements  Cij  are  known  [86]: 


Qi  =  '£c,iVi  (3.25) 

j 

Now  consider  the  charge  on  the  conductor  i  when  all  voltages  except  the  voltage  on 
conductor  j  are  zero.  The  charge  on  conductor  i  is  then: 

Qi  =  CijVj  (3.26) 

A  positive  voltage  on  conductor  j  induces  a  negative  charge  on  conductor  i.  This 
means  that  for  the  mutual  capacitances,  Cij  <  0.  However,  in  the  case  of  self¬ 
capacitances,  i  =  j,'we  have  Cji  >  0  since  a  positive  voltage  induces  a  positive  charge 
on  the  same  conductor.  Electrostatic  theory  leads  us  to  the  following  relationships  be¬ 
tween  the  capacitance  matrix  elements  and  the  circuit  parameter  capacitances,  which 
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are  defined  to  be  positive  numbers: 


Cdi  =  -Cdi  for  i^d  (3.27) 

Odd  =  +  =  C'e  • 

i 

Using  eq.  (3.22),  the  chemical  potential  as  a  function  of  electrode  voltages  {V^}  is 


(3.28) 


and,  using  eq.  (3.23),  the  gate  voltages  at  which  the  conductance  peaks  occur  are 


Vg  y  Vi  Cdi 

0  f^dg  i:^g  0  Cdg 


(3.29) 


The  gate  voltage  difference  between  adjacent  conductance  peaks  is,  similar  to 
eq.  (3.24), 


AK 


(3.30) 


3.2.1  Energy  Level  Separation  within  the  Dot 

In  a  more  rigorous  treatment,  one  has  to  take  the  energy  levels  within  the  quantum 
dot  into  account,  which  arise  due  to  the  small  spatial  extent  of  the  dot  [85].  Using 
eq.  (3.15),  the  conductance  condition  becomes 

U{N)-U{N -T)  =  yL{N)-EN  .  (3.31) 

The  chemical  potential  in  the  dot  is  then 

^(N)  =  Em  +  U(N)  -  U{N  -1)  =  Ek  +  (N-  i)  a 


(3.32) 
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and  the  gate  voltage  for  a  conduction  peak  is 

Vg  _  N  —  ^  ^  Vj  Cdi  Cs  En 
c  Cdg  6  Cdg  Cdg  6^ 

The  distance  between  adjacent  conductance  peaks  is 


(3.33) 


AV, 


e  ^  CsAE 
Odg  ^dg  ^ 


(3.34) 


The  spin  degeneracy  of  the  levels  is  lifted  by  the  charging  energy.  The  distance 
between  adjacent  peaks  alternates  between  /Cdg  and  /Cdg  +  {CT./Cdg){^E/&)- 
If  the  energy  difference  AE  is  large  compared  to  the  thermal  energy,  only  one 
energy  level  within  the  well  contributes  to  conduction  and  one  approaches  the  limit 
of  resonant  tunneling  through  a  well.  If  AJE^  is  larger  than  the  incremental  charging 
energy,  resonant  tunneling  is  the  predominant  effect.  Resonant  tunneling  occurs  if 
the  Fermi  level  in  the  leads  matches  the  energy  level  of  a  bound  state  within  the  well. 


3.2.2  Thermal  Spread 

The  behavior  of  the  conductance  peaks  with  temperature  was  investigated  by  Bee- 
nakker  [85].  Two  limiting  cases  can  be  given  for  the  conductance  as  a  function  of 
chemical  potential,  depending  on  the  size  of  the  energy  level  separation  NE  compared 
to  the  thermal  energy  ksT.  First,  we  consider  the  lineshape  in  the  resonant  tunneling 
regime,  ksT  AF'  e^/C-^,.  The  thermal  energy  is  smaller  than  the  energy  level 
separation,  as  can  be  the  case  for  very  low  temperatures  and  small  quantum  dots. 
The  conductance  is 


G  = 


■pl  pr 
2  P  P 

F'  +  F’' 


f'it^p) 


r],r; 

4/cbTF',  +  F; 


cosh 


_Jh_\ 

2kBT’  ’ 


(3.35) 


where  /'  is  the  derivative  of  the  Fermi-Dirac  function.  F'  and  F^  are  the  tunneling 
rates  for  the  left  and  right  barriers,  respectively.  For  a  sweep  of  gate  voltage,  we  fit 
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the  conductance  using  the  formula 

G(K,)=EG'«cosh-^(5^5i-j^)  ,  (3.36) 

where  K-i  are  the  position  of  conductance  resonances  with  peak  conductance  Gri-  Note 
that  the  capacitive  lever-arm  a  =  Cdg/C^  causes  a  spreading  of  the  peak,  causing 
distinguishable  peaks  to  disappear  for  large  Cs. 

For  the  classical  limit,  when  AE  -C  ksT  e^fC-z,  a  continuum  of  energy  levels 
contributes  to  the  conductance  and  the  lineshape  is  given  by  [85,87] 


n/ksT 

’^“sinh(/x/itBT) 

r'F’' 

2A£;r'-hr'' 


-  Gmax  cosh" 


-2( - 1± - ) 

^2AkBT’ 


(3.37) 


For  a  sweep  of  the  gate  voltage,  the  conductance  can  be  fit  using 

G(V,)  =  S  .  (3.38) 

Although  the  line  shapes  in  the  classical  and  resonant  tunneling  regimes  are  different, 
they  are  not  distinguishable  if  temperature  is  used  as  a  fitting  parameter. 

3.3  Low  Temperature  Measurements 

A  schematic  of  the  refrigerator  used  to  measure  devices  at  cryogenic  temperatures  is 
shown  in  Fig.  3-6.  Using  this  probe,  it  is  possible  to  cool  a  sample  to  300  mK  within 
a  few  hours  and  to  keep  it  at  base  temperature  for  more  than  8  hours.  Following  is  a 
short  explanation  of  how  the  probe  works.  More  extensive  operating  instructions  for 
this  particular  refrigerator  are  given  in  ref.  [78]. 

The  probe  is  immersed  in  liquid  ^He,  keeping  the  outside  of  the  inner  vacuum  can 
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Figure  3-6:  Schematic  of  the  Oxford  Instruments  Heliox  insertion  probe  by  Oxford 
Instruments.  Temperatures  of  300  mK  are  routinely  achieved  with  a  hold  time  of 
more  than  8  hours. 
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at  a  constant  temperature  of  4.2  K.  The  drinking  straws  are  immersed,  as  indicated 
in  Fig.  3-6.  Cooling  is  achieved  by  reducing  the  vapor  pressure  of  liquid  ^He  using  a 
sorb  pump  which  adsorbs  gas  if  cooled  below  a  critical  temperature  of  about  12  K.  To 
initiate  a  cooling  cycle,  liquid  “^He  is  admitted  into  the  1 K  pot  and  pumped  on  using 
an  external  rotary  pump.  This  reduces  the  temperature  of  the  '^He  to  about  1.2  K. 
The  sorb  pump  is  then  heated  above  its  critical  temperature,  while  the  1 K  pot  and 
the  tube  that  connects  sorb  and  0.3  K  pot  is  kept  cold  using  the  sorb  drinking  straw. 
The  ^He  that  is  released  during  heating  of  the  sorb  is  cooled  to  1.2  K  while  passing 
through  the  1 K  pot,  and  condenses  in  the  0.3  K  pot.  The  sorb  heater  is  then  turned 
off  and  the  temperature  of  the  0.3  K  pot,  which  is  thermally  shorted  to  the  sample 
block,  reaches  the  base  temperature  of  about  300  mK  within  a  few  minutes.  The 
sample  is  glued  to  a  44-pin  chip  carrier  using  silver  paint  and  bonded  with  aluminum 
wire.  The  back  of  the  chip  carrier  is  thermally  sunk  to  the  sample  block.  The  probe  is 
wired  with  manganin  leads  which  are  thermally  shorted  to  the  0.3  K  pot.  Manganin 
is  a  material  of  low  thermal  but  high  electrical  conductivity. 

The  conductance  of  the  device  is  measured  by  applying  a  DC  or  low  frequency  AC 
voltage  to  an  ohmic  contact  (drain)  and  measuring  the  current  through  another  ohmic 
contact  (source).  Because  the  currents  are  very  small,  it  is  important  to  minimize  the 
noise  level.  This  is  also  important  because  any  noise  on  the  wires  can  heat  the  sample 
and  deteriorate  device  performance.  For  this  reason,  all  electrical  lines  that  are  used 
for  device  and  temperature  measurements  must  be  filtered.  In  order  to  measure  the 
high-impedance  devices,  electronics  had  to  be  designed  and  built.  The  design  of  the 
filter  box  is  described  in  Appendix  H.  For  AC  measurements,  a  PAR  5210  model 
lock-in  amplifier  was  used.  Constant  gate  voltages  were  set  using  battery-operated 
voltages  sources,  and  the  variable  gate  voltage  was  set  using  an  HP  3325B  function 
generator.  A  computer-aided  data  acquisition  system  which  controls  the  lock-in  by 
an  IEEE-488  bus  was  used  for  the  measurements. 
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3.4  Measurement  of  Quantum  Dot 

A  device  was  built,  consisting  of  four  consecutive  quantum  point  contacts  which  can 
be  controlled  independently.  An  SEM  photograph  of  a  monitor  sample  with  identical 
geometry  is  shown  in  Fig.  3-7.  A  quantum  dot  is  produced  by  biasing  two  quantum 


Figure  3-7:  SEM  photograph  of  the  completed  device.  The  fingers  of  the  quantum 
point  contacts  have  a  width  of  approximately  70  nm.  A  quantum  dot  can  be  produced 
by  biasing  two  quantum  point  contacts  below  the  the  first  subband,  thus  creating  a 
tunnel  barrier  for  electrons  in  the  dot. 

point  contacts  below  the  first  subband.  This  allows  electron  transport  into  and  out 
of  the  dot  only  via  tunneling.  The  remaining  gates  can  be  used  to  sweep  the  chemical 
potential  within  the  dot.  The  size  of  the  quantum  dot  can  be  adjusted  from  less 
than  200  nm  to  approximately  600  nm  by  pinching  off  the  different  pairs  of  quantum 
point  contacts.  A  schematic  showing  how  the  different  dot  sizes  are  achieved  is  shown 
in  Fig.  3-8.  The  finished  device  has  minimum  linewidths  below  70  nm.  The  small 
linewidth  causes  the  capacitances  of  the  gates  to  the  dot  to  be  small.  According 
to  eq.  (3.12),  the  temperatures  at  which  observation  of  conductance  oscillations  is 
possible  should  be  higher  than  for  coarser  gate  structures. 
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Figure  3-8:  Schematic  of  biasing  conditions  for  different  dot  sizes.  Dots  of  various  sizes 
can  be  made  by  pinching  off  two  pairs  of  the  four  consecutive  quantum  point  contacts. 
Each  pair  of  quantum  point  contacts  has  an  opening  of  300  nm,  the  distance  from 
one  pair  to  the  next  is  300  nm.  The  maximum  possible  length  of  the  dot  is  600  nm. 

A  schematic  of  a  typical  measurement  setup  is  shown  in  Fig.  3-9.  An  AC  voltage,  in 
this  case  100  /xV,  is  imposed  on  the  source  of  the  device,  and  the  current  is  measured 
at  the  drain.  To  bias  one  quantum  dot,  two  of  the  four  pairs  of  quantum  point 
contacts  are  biased  in  the  tunneling  regime.  The  remaining  four  gates  can  be  used 
either  to  sweep  the  chemical  potential  in  the  dot  or  kept  constant.  Since  all  gates 
can  be  controlled  independently,  a  large  variety  of  experiments  can  be  conducted. 
In  addition  to  several  sizes  of  quantum  dot,  as  described  above,  coupling  of  small 
quantum  dots  and  transitions  from  large  to  small  dots  can  be  investigated. 

The  conductance  oscillations  are  strongest  when  both  barriers  are  tuned  to  be 
identical  in  transmission.  To  fine-tune  the  barriers,  the  following  method  can  be 
used.  One  of  the  barriers  should  be  incrementally  made  stronger.  If  the  threshold 
of  the  device  shifts  positive,  the  stronger  barrier  has  been  made  stronger.  If  the 
threshold  stays  the  same,  which  is  the  desired  case,  the  weaker  barrier  has  been  made 
stronger,  and  this  generally  results  in  stronger  conductance  oscillations.  An  example 
of  a  conductance  measurement  versus  gate  voltage  for  the  biasing  conditions  indicated 
in  Fig.  3-9  is  shown  in  Fig.  3-10. 

It  should  be  noted  that  due  to  the  small  spatial  extent  of  the  device  it  is  not 
trivial  to  bias  a  dot.  Large  Coulomb  blockade  devices  can  be  biased  by  determining 
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Figure  3-9:  Schematic  of  the  measurement  setup  used  to  measure  conductances. 


Vg[V] 


Figure  3-10:  Plot  of  conductance  versus  gate  voltage  for  a  quantum  dot,  having  a 
length  of  600  nm.  The  two  outermost  quantum  point  contacts  define  the  dot.  For 
comparison  of  scales,  the  value  of  quantized  conduction  is  2e'^/h  ~  77/xS. 
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the  voltage  of  the  individual  quantum  point  contact  at  which  only  tunneling  takes 
place,  and  then  using  this  voltage  to  bias  a  quantum  dot.  The  close  proximity  between 
the  gates  in  this  structure  causes  a  strong  influence  of  the  change  of  one  tunneling 
barrier  on  an  adjacent  barrier.  For  this  reason,  biasing  has  to  be  done  by  iteration  for 
very  small  quantum  dots.  Another  consequence  of  this  coupling  of  tunneling  barriers 
to  nearby  electrodes  is  that  an  electrode  that  is  used  to  sweep  the  chemical  potential 
within  the  dot  also  alters  the  tunneling  barriers  and  shape  and  size  of  the  dot.  This 
means  that  conductance  oscillations  can  be  seen  over  smaller  ranges  of  gate  voltage 
than  for  larger  devices. 

3.4.1  Gate  Capacitances 

For  the  largest  possible  dot,  with  600  nm  distance  between  the  center  of  the  tunnel¬ 
ing  barriers,  four  gates  access  the  middle  of  the  dot.  Each  and  any  combination  of 
these  four  gates  can  be  swept  and  reveals  a  capacitance  for  the  particular  gates  to 
the  quantum  dot.  Figure  3-11  shows  a  plot  of  the  conductance  for  the  the  dot  with 
one,  two,  three,  and  four  gates  simultaneously  swept.  For  two  gates,  adjacent  gates 
were  swept  simultaneously.  This  experiment  proves  that  the  gates  are  independent  of 
each  other  and  that  each  gate  finger  contributes  approximately  the  same  capacitance 
of  ^ibout  16  aF  to  the  dot-to-gate  capacitance.  All  four  center  gates,  swept  inde¬ 
pendently,  show  an  almost  identical  periodicity  of  conductance  oscillations,  which 
indicates  close  matching  of  their  geometries. 

3.4.2  Measurement  of  a  Variable  Size  Quantum  Dot  Device 

Different-sized  dots,  indicated  in  Fig.  3-8,  can  be  measured  by  biasing  different  quan¬ 
tum  point  contacts  for  each  measurement.  The  conductances  for  600,  400,  and  200  nm 
distance  between  the  defining  quantum  point  contacts  is  plotted  in  Fig.  3-12  as  a  func¬ 
tion  of  gate  voltage.  In  case  of  the  200  nm  dot,  a  gate  outside  the  tunneling  barriers 
was  swept.  The  actual  size  of  the  dots  is  smaller  than  the  distance  between  quantum 
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Figure  3-11:  Conductance  of  a  600  nm  dot  with  one,  two,  three  and  all  four  inner 
gates  swept. 

point  contacts  by  approximately  half  the  length  of  both  tunneling  barriers. 

Considering  that  the  capacitance  is  related  to  the  period  of  conductance  oscilla¬ 
tions,  as  given  in  eq.  (3.30),  we  can  assign  a  dot-to-gate  capacitance  for  the  400  nm 
dot  of  ~8  aF  and  for  the  600  nm  dot  of  ~16  aF.  For  the  smallest  dot,  we  observe  that 
the  spacing  between  the  maxima  is  not  constant,  possibly  indicating  a  suflSciently 
large  energy  level  separation  according  to  eq.  (3.34).  The  capacitance  is  then  deter¬ 
mined  by  the  smallest  spacing  between  conductance  peaks  and  can  be  estimated  to 
be  ~2  aF. 

Fit  with  Fermi  Function 

A  fit  for  the  conductance  of  the  400  and  200  nm  dot  with  the  lineshape  given  in 
eq.  (3.36)  is  shown  in  Fig.  3-13.  It  should  be  noted  that  the  fit  is  only  reasonably 
good  for  the  first  few  conductance  peaks.  For  the  remaining  peaks,  the  off-resonance 
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Figure  3-12:  Conductance  of  quantum  dots  with  separation  of  600,  400,  and  200  nm 
between  tunneling  barriers.  In  all  cases  only  one  of  the  gates,  indicated  as  the  black 
gate  in  the  inserts,  is  swept.  In  the  case  of  the  200  nm  dot  a  gate  beyond  the  tunneling 
barrier  was  swept.  The  curves  are  offset  by  10  fiS  for  clarity. 


Figure  3-13;  Fit  of  measured  conductances  for  the  400  and  200  nm  dot  with  eq.  (3.36). 
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200  nm  (a) 

400  nm  (6) 

C,i  [aF] 

2.0 

8.1 

Ci^[aF] 

310 

400 

Table  3.1:  Capacitances  for  200  nm  dot  (a)  and  400 nm  dot  (ft).  The  sum-capacitances 
are  extracted  from  the  fit  and  the  gate  capacitances  from  the  periodicity  of  the  con¬ 
ductance  peaks. 

conductance  becomes  large.  The  fit  with  the  Fermi  function  enables  us  to  extract 
an  approximate  value  for  the  capacitive  lever-arms  of  the  dots.  If  the  fit  for  the 
resonant  tunneling  regime  is  used,  the  capacitive  lever-arm  for  the  200  nm  dot,  o,  is 
~  1.3  X  10-2,  and  for  the  400  nm  dot,  6,  is  at,  ~  2.0  x  lO'^.  It  must  be  noted 
that  the  lever-arms  for  the  system  are  smaller  than  expected.  If  the  remaining  gates 
have  a  capacitance  of  the  same  order  of  magnitude,  the  total  gate  capacitance  for 
the  400  nm  dot  would  still  be  below  100  aF.  Even  if  one  considers  the  capacitances 
to  the  leads  and  to  the  substrate,  whose  conducting  back-gate  is  separated  from  the 
2DEG  by  the  thickness  of  the  wafer,  it  is  difficult  to  imagine  the  capacitance  to  reach 
400  aF.  It  is  possible  that  the  experiments  were  conducted  at  higher  temperature  than 
was  measured  at  the  sample  block  [88].  However,  the  ratio  between  the  lever-arms, 
important  for  discussion  of  the  two  dots  measured  in  series,  should  be  approximately 
correct.  The  values  for  the  capacitances  for  the  two  dots  are  given  in  table  3.1.  The 
capacitance  given  for  dot  b  is  extracted  from  the  period  of  the  first  few  conductance 
maxima.  For  higher  gate  voltages,  the  dot  becomes  larger  in  size,  resulting  in  a  larger 
capacitance  of  about  9  aF  and  higher  periodicity. 


3.5  The  Double  Dot 

The  device  shown  in  Fig.  3-7  can  be  biased  to  produce  two  quantum  dots  in  series  by 
using  three  pairs  of  quantum  point  contacts  to  produce  three  consecutive  tunneling 
barriers.  A  schematic  of  the  system  is  shown  in  Fig.  3-14.  The  two  dots  were  already 
measured  individually  in  Sec.  3.4  and  are  denoted  as  o  for  the  small,  and  b  for  the 
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Figure  3-14;  Schematic  of  the  double  dot  structure.  The  center-to-center  distance 
between  tunneling  barriers  is  200  nm  and  400  nm  for  the  two  dots.  The  gate  indicated 
in  dark  is  swept. 


large  dot. 

Systems  containing  multiple  dots  have  been  investigated  theoretically  for  a  couple 
of  years  [89-93],  but  experimental  studies  are  just  emerging.  Some  experiments  focus 
on  identical  or  slightly  detuned  quantum  dots  [94],  while  others  focus  on  dots  with 
identical  physical  layout  which  are  strongly  detuned  through  gate  bias  [95].  The 
experiments  described  in  this  work  use  dots  of  unequal  size  and  significantly  unequal 
capacitances.  If  one  stretches  the  physical  analog  of  one  quantum  dot  being  an  atom 
with  its  confined  states  [77],  the  analog  of  a  double  dot  is  a  molecule.  The  coupling  of 
the  two  atoms  of  this  molecule  can  be  modified  by  simply  changing  the  gate  voltage 
of  the  tunneling  barrier  separating  the  dots.  We  expect  the  physical  characteristics  of 
this  double  dot  to  be  different  from  the  characteristics  of  the  isolated  dots.  It  would 
be  of  great  interest  to  map  the  energy  level  spectrum  of  such  structures. 

The  derivation  for  the  conductance  condition  for  a  double  dot  is  done  in  close 
analogy  to  the  derivation  for  the  single  dot  of  Sec.  3.2.  Equations  (3.17)  to  (3.23) 
for  the  single  dot  correspond  to  eqs.  (3.39)  to  (3.50)  for  the  double  dot.  It  is  hoped 
that  this  close  analogy  will  avoid  confusion  over  the  complexity  of  the  formulas.  The 
meaning  of  the  results  will  be  further  explained  by  considering  special  cases.  In  this 
derivation,  we  again  neglect  the  energy  level  separation  caused  by  the  finite  extent  of 
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the  dots,  keeping  in  mind  that  this  can  be  accounted  for  as  in  Sec.  3.2.1.  This  deriva¬ 
tion  considers  the  potentials  and  capacitances  of  both  dots  to  an  arbitrary  number 
of  electrodes.  This  treatment  is  necessary  for  a  quantitative  analysis  of  measurement 
results. 

Consider  a  double  dot,  with  the  two  dots  denoted  by  indices  a  and  b,  surrounded 
by  an  arbitrary  number  of  electrodes  i.  The  change  in  charging  energy  in  the  system 
is 

5U  =  VaSQa  +  VbSQb  +  Vi6Qi  .  (3.39) 

i 

The  first  two  terms  represent  a  change  in  energy  of  the  two  dots.  The  charges  Qi  on 
the  electrodes  are  not  controlled  variables.  To  get  an  equation  for  the  internal  energy 
as  a  function  of  the  controlled  variables  Qa,  Qbi  ^tnd  Vi,  we  use  again  the  Legendre 
Transformation  tj  =  U  —  'Tli  ViQi  to  get 

SU  =  VaSQa +  Vt5Qt-Y.Qi6Vi  .  (3.40) 

i 

We  again  rename  U  as  the  charging  energy  U .  The  experiment  is  done  with  fixed 
{Vi},  so  that 

/*Q  rQb 

I7|(W=0)  =  I  ■  VM,  Q,  =  0,  {V5})5Q1  + 1  H(Q.,  Qi  •  (3.41) 

To  find  V;  and  Vb  as  a  function  of  the  variables  Qa  and  Qb,  we  use  again  the  capaci¬ 
tance  matrix  elements  Camj 

Qa  =  E(^-rnV,n  .  (3-42) 

m 

Solving  for  the  potential  Va  we  get 

V^  =  [Qa  CamVm]  j 

m^a 


(3.43) 
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and  the  dot  potentials  14  and  Vj,  become 


14  = 

— 0. 

1 

1 

I 

1 

1 

1  f>i. 

^aa 

^aa 

i  ^aa 

—Qb 

-^Va-y—\ 

Cbb 

Cbb 

i  Cfcb 

(3.44) 


Solving  for  14  and  14  by  elimination,  we  get 


K 

14 


^  -Qa-  ^  ^ 


1 


^ab  ^ 


V-— - 

Lr  n  . 


^  Can 


KC^b 


Qb  — 


1  1 


K 


^ab  ^ 


Qa 


14  Cit  ^ 

i  «  Cbb  A 


(3.45) 


where  the  parameters  k,  ai,  and  A  are  functions  of  the  capacitance  matrix  elements 
and  are  defined  as 


K  = 


Pr'  = 


1  - 


''ab 


^aa  ^bb 

^bb  ^ai  ^ai  ^ab 

^ab  ^ai  _  ^  ^ai  ^bb 

^aa  ^bi  ^bi  ^ab 


(3.46) 


Using  eq.  (3.41),  we  integrate  over  the  charge  in  both  dots  and  get  for  the  charging 
energy: 


Qt  I  Qb  QaQb  1 _ ^  ^  I  Cfcj  1  ^ 

2c^fl/C  ^Cbbbb  (^ab  ^  j  ^clq.  ^bb  A^ 


(3.47) 


The  chemical  potential  can  be  calculated,  using  eq.  (3.16),  as 


PaiNa,  Nb)  =  U{-eNa,  -eNb)  -  t/(-e(7V„  -  1),  -eNb)  , 


(3.48) 
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and  likewise  for  Nb): 


M.N,)  = 

A/Ufla  ^  UqJ)  n>  ^  L^aa 

^  ^a6  ^  I 


(3.49) 


Cbb/^iK 


The  conductance  is  measured  for  negligible  drain  voltage,  so  source  and  drain  are 
assumed  to  be  at  zero  potential.  Then  =  /^h  =  0  give  the  conditions  for  current 
flow.  If  only  one  of  the  chemical  potentials  in  the  dots  lines  up  with  the  chemical 
potential  in  the  wires,  tunneling  into  and  out  of  the  dot  can  still  happen,  but  current 
flow  is  suppressed.  If  we  define  one  of  the  electrodes  i  as  the  gate  g,  with  voltage  Vg, 
the  conditions  for  current  flow  are 


Ik 

e 


- ^ag - (1  -  Kjag  +  2^ - ^ 

Cag  ^ag  i:^g  ^  ^ag 

(1  -  k)/J,  + 

Cbg  Cab  Cbg  i^g  ^  ^bg  Pi 


(3.50) 


Circuit  Parameter  Representation:  Using  the  formalism  developed  in  Sec.  3.2, 
the  mutual  and  self-capacitances  are  related  to  circuit  elements  by  the  relations 


Cai  =  -Cai  for  z  o 
Cbi  =  -Cbi  for  i^b 

Caa  =  Oab  +  ^  Cat  =  Cab  +  Ca 

i 

Cbb  =  Cab  +  ^2  ^bi  =  Cab  +  Cb  • 


(3.51) 


With  this,  the  parameters  k,  Oj,  /?i,  and  a  newly  defined  C  become 
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{Cab  +  Ca){Cab  +  Cb) 

{Cgb  +  Cb)Cai 

{Cab  Cb)Cai  +  CabCbi 
{Cgb  +  Ca)Cbi 
{Cab  +  Ca)Cbi  +  CabCai 
CabCa  +  CabCb  +  CgCb  • 


(3.52) 


The  charging  energy  is 


Ql{CabVCb)VQl{Cab  +  Ca)  +  2QaQbCab  ^ 

2C2 


(3.53) 


^  {Cab  +  Cb)Cai  +  CabCbi  ,  ^  {Cgb  +  Ca)Cbi  +  CgbCgi  ^ 

- +  0b - - )  ■ 


The  chemical  potential  is  then 


f^a{Nat  Nb') 


+  1,  ,  e^C, 


-  5)  + 

{Cab  +  Cij)Cai  +  CabCbi 


(3,54) 


e^Cab  +  Cg),^^  T  ,  e^C, 


C2  ^  ^  2’  C"^ 

r  {Cgb  +  Ca)Cbi  +  CgbCgi 


and  the  conditions  for  conduction  3.50  become 


NbCgb 

Cag{Cab  +  Cb)  +  CgbCbg 


Ng-^  Cgg{Cgb  +  Cb)  NbCgb 

Cgg  {Cgb  Cb)Cag  +  CgbCbg  C ag{Cab  +  Cb)  +  CgbCl 

^  Vi  {Cgb  +  Cb)Cai  +  CgbCbi 
e  Cgg{Cab  +  Cb)  +  CgbCbg 

Nb  —  \  Cbg{Cab  +  Cg) _ I _ NgCgb _ 

Cbg  {Cgb  +  Ca)Cbg  +  CgbCgg  Cbg{Cab  +  Cg)  +  C gbC , 
^Vj  {Cgb  +  Cg)Cbi  +  CgbCgi 
i^g  e  Cbg{Cab  +  Cg)  +  CgbCgg 


(3.55) 
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These  are  the  general  expressions  for  the  charging  energy  and  conductance  conditions 
in  a  double  dot  system  with  an  arbitrary  number  of  electrodes.  In  the  following 
paragraphs  we  discuss  some  special  cases  which  can  be  deduced  from  these  equations. 


Special  Case:  Single  Dot  Limit  (Cab  ->  oo)  In  the  limit  of  very  strong  coupling 
we  should  get  the  solution  for  one  single  dot  with  Na  +  Nt  electrons  and  capacitance 
Cag+Cbg.  The  expressions  for  charging  energy,  chemical  potential  and  voltage  become: 


Vg  _  Ng  +  Nb  —  ^  ^  Vj  Cqj  +  Cbj  5g\ 

e  Cag  +  Cbg  e  Cag  +  Cbg 

It  can  easily  be  verified  that  the  last  two  expressions  are  identical  to  eqs.  (3.28) 
and  (3.29),  with  Na  +  Nb  =  N  and  Cgi  +  Cu  =  Q. 


Special  Case:  Independent  Dot  Limit  (Coj,  0)  In  the  limit  of  no  coupling 
between  the  dots.  Cab  ->  0,  the  expressions  for  the  charging  energy,  the  chemical 
potential  and  the  voltage  conditions  for  current  transport  become: 

=  ^  +  +  (3-59) 


Cu 


(3.60) 


/AT 


1 
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Note  that  in  this  case  Cq  =  Cas,  because  the  coupling  capacitance  Cab  is  neglected. 
There  are  two  conditions  that  have  to  be  fulfilled  independently  and  there  is  no 
coupling  between  the  two  dots.  The  potential  of  each  dot  is  independent  of  the 
potential  in  the  adjacent  dot. 

Special  Case:  Ruzin  et  al.  If  a  capacitance  C  connects  the  dots  to  each  other 
and  to  the  leads,  as  indicated  in  figure  3-15,  we  should  get  the  same  result  as  Ruzin, 
Glazman  et  al.  [89-91].  We  are  considering  this  as  a  special  case  to  the  formalism 
developed  above.  Appendix  J  derives  the  same  formulas  using  Ruzin’s  approach. 
Since  we  want  to  compare  the  results  to  the  formulas  given  by  Ruzin,  we  adopt  the 
same  notation  (a  — >  1  and  6  — >  2).  In  this  case  we  have 

Cab  =  C 

Ca  =  C  +  Ci 

Cb  =  C  +  C2 

C^  =  3C^  +  2CxC  +  2C2C^CiC2  . 


(3.62) 
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For  the  charging  energy  we  get 


Ql{2C  +  C2)  +  Ql{2C  +  Ci)  +  2QiQ2C  ^ 

^  "  2C2 

_ 2CCi-\-CiC2  +  CC2  ,  ^  2CC2  +  CiC2  +  CCi^ 

VgiQi - ^ - +  Q2 - )  • 


(3.63) 


Equation  (3.63)  should  be  identical  to  equations  (8)-(10)  in  I.  M.  Kuzin’s  paper  [89]^ 
The  two  conditions  for  current  flow  result  in  the  chemical  potentials 


e^{Ni  —  5)  (2(7  +  Cl)  e^N2C  {2CCi  +  CC2  +  C1C2) 

=  -^2  (72 


e\N2  -  1){2C  +  C2)  e^iViC  ,,  {2CC2  +  CCi  +  C1C2) 
- - +  -g5--V  gi 


(3.64) 


and  for  the  voltages,  assuming  the  chemical  potentials  ni  and  fX2  are  zero, 

2C^C2  C 

e  “  ^  ^  2^2CCi  +  (7(72  +  (7iC2  ^2(7Ci  +  (7(^2  +  (71(72 

Yl  _  fiv  S  2C  +  (7i _  _ C - 

e  “  2^2(7(72 +  (7Ci  +  (7i (72  ^2C(72  +  (7Ci  +  Ci(72 


(3.65) 


A  comparison  with  Appendix  J  shows  that  the  two  derivations  have  identical  results. 


Special  Case:  Capacitance  Between  Dots  In  this  case  we  neglect  all  capaci¬ 
tances  except  the  capacitances  from  the  gate  to  the  dots  and  the  coupling  capacitance, 
as  shown  in  Fig.  3-16. 

1  There  axe  two  errors  in  I.  M.  Ruzin’s  paper,  the  first  one  is  a  sign  error  in  eq.  (6),  the  second 
an  indexing  error  in  eq.  (9). 
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Gate 


Figure  3-16:  Easiest  case  for  coupling  of  two  quantum  dots. 


For  the  capacitances  in  this  system  we  get 


Ca 

Cb 

& 


=  c, 


ag 


(3.66) 


=  Cb, 

“  CdbCdg  CdbCbg  "h  CdgCbg 


For  the  charging  energy  we  get 

TT  _  Qli^ab  +  Cbg)  +  QbiCdb  +  Cdg)  +  2QdQbCdb  I  7/  I  ^  \ 
C  —  n-unr'j-r'r'\  *g\Oa  +  Ob) 

^[^ab^ag  •  ^ab^bg  >  ^ag^bg ) 


(3.67) 


TMs  is  identical  to  eq.  (J.6),  which  was  derived  using  Ruzin’s  formalism. 

For  the  chemical  potential  in  the  dots  and  the  conditions  for  conduction  we  get 


fJ'ai^Nd,  Nb) 
f^b{,Ndi  Nb') 


e^{Nd  -  |)(Ca6  +  Cbg)  +  e'^NbCdb 
CdbCdg  A  CdbCbg  “|-  CdgCbg 

e^Nb  -  l){Cdb  +  Cdg)  +  e'^NdCdb 

CdbCdg  “I"  CdbCbg  A  CdgCbg 


-eVd 


eV„ 


(3.68) 


V 


9  _ 


v„ 


{CdbACbg){Nd-^^)ACdbNb 
CdbCdg  “h  CdbCbg  “1“  CdgCbg 
{CdbACdg){Nb-l)ACdbNd 
CdbCdg  ”1“  CdbCbg  "1“  CdgCbg 


(3.69) 
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Figure  3-17:  Energy  ladders  for  a  double  dot  system.  The  left  ladder  shows  the 
smaller  dot,  the  right  the  larger  dot.  Note  that  the  two  ladders  move  with  diflPerent 
speed  when  the  gate  voltage  is  changed,  according  to  the  difference  in  capacitive  lever 
arms. 

3.5.1  Position  of  Conductance  Peaks 

A  double  dot  system  can  only  conduct  current  if  the  chemical  potentials  of  both 
dots  are  aligned  with  the  chemical  potentials  in  the  leads.  A  simple  picture  of  the 
problem  is  shown  in  Fig.  3-17.  In  the  case  of  no  coupling  between  the  quantum  dots 
we  have  two  rigid  energy  scales  which  move  with  respect  to  each  other,  but  whose 
movement  is  only  dependent  on  the  gate  electrode  rather  than  the  charge  on  the 
adjacent  dot.  Conduction  through  the  double  dot  can  be  sustained  if  the  chemical 
potentials,  plotted  here,  are  within  about  ksT  of  each  other  and  of  the  leads. 

Because  the  levels  move  with  respect  to  each  other,  it  is  not  easy  to  extract 
information  on  the  location  of  conductance  maxima  from  Fig.  3-17.  Instead,  one  can 
plot  a  diagram  of  charging  energy  as  a  function  of  gate  voltage,  as  done  in  Fig.  3-4 
for  the  single  dot.  We  use  eq.  (3.53)  to  plot  the  charging  energy  as  a  function  of 
gate  voltage  for  the  double  dot  system.  The  charging  energy  has  to  be  plotted  with 
both  Na  and  Nb  as  parameters.  For  identical  dots  and  no  coupling  capacitance,  the 
states  of  lowest  energy  have  identical  numbers  of  electrons  in  each  dots,  as  shown  in 
Fig.  3-18.  As  can  be  seen  in  Fig.  3-19,  an  unequal  number  of  electrons  in  the  dots  is 
only  possible  at  the  charge  degeneracy  points.  The  energy  levels  for  (2,3)  and  (3,2) 
are  degenerate  for  any  gate  voltage.  At  the  intersection  of  (2,2)  and  (3,3),  current 
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Figure  3-18:  Charging  energy  as  a  function  of  gate  voltage  for  two  identical  dots, 
with  Cag  =  Cbg  =  10  aF. 


Figure  3-19:  Close-up  of  Fig.  3-18.  The  energy  levels  that  contribute  to  current  flow 
are  labeled. 
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Figure  3-20:  Close-up  of  charging  energy  as  a  function  of  gate  voltage  for  two  identical 
dots,  in  case  of  coupling  between  the  dots. 

can  flow  through  the  system,  while  the  system  goes  through  the  two  possible  cycles: 

(2.2)  (3,2)  ^  (2,3)  (2,2),  or  (3,3)  (3,2)  ^  (2,3)  ^  (3,3).  We  have  assumed  that 

the  stable  states  are  (2,2)  and  (3,3)  and  that  electrons  tunnel  from  the  left  to  the 
right.  If  a  coupling  capacitance,  in  this  case  of  1  aF,  is  introduced  into  the  problem, 
the  energy  levels  (2,3)  and  (3,2)  are  still  degenerate,  but  shifted  to  a  slightly  lower 
energy,  as  shown  in  Fig.  3-20.  This  causes  an  intermediate  state,  (2,3)  and  (3,2), 
to  be  most  favorable  for  a  small  range  of  gate  voltages.  For  the  case  of  coupling, 
the  conductance  peak  splits  into  two,  the  first  one  with  the  cycle  involving  the  state 

(2.2) ,  and  the  second  peak  involving  the  cycle  with  the  state  (3,3).  This  splitting  was 
recently  observed  experimentally  by  Waugh  et  al.  [94].  If  there  is  a  detuning  between 
the  dot-to-gate  capacitances,  the  degeneracy  of  the  (2,3)  and  (3,2)  level  is  lifted,  as 
shown  in  fig  3-21.  For  conductance  to  occur  in  this  case,  the  electron  would  have  to 
be  thermally-excited  because  the  state  (3,2)  is  necessary  for  conduction  through  the 
double  dot  system. 

If  the  gate  capacitances  strongly  deviate  from  each  other,  which  is  the  case  for 
unequal  dots,  the  number  of  electrons  added  for  a  sweep  of  the  gate  voltage  is  different 
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Figure  3-21:  Close-up  of  charging  energy  as  a  function  of  gate  voltage  for  two  identical 
dots,  in  case  of  coupling  and  5%  detuning  between  the  dots. 

for  each  dot.  For  the  case  Cag  =  10  aF,  Cbg  =  5Cag  and  Qs  =  lOCip,  the  charging 
energy  as  a  function  of  gate  voltage  is  plotted  in  Fig.  3-22.  As  expected,  for  one 
electron  added  to  the  smaller  dot,  five  electrons  have  to  be  added  to  the  larger  dot.  A 
close-up  of  the  transition  from  the  state  {Na,  Nf,)  =  (1,7)  to  (2, 8)  is  shown  in  Fig.  3- 
23.  Because  we  chose  the  ratio  of  the  capacitances  to  be  an  integer,  the  states  with 
(1,7)  and  (2,8)  happen  to  be  degenerate.  The  conduction  mechanism  is  identical  to 
the  one  discussed  for  identical  dots.  For  slightly  different  choices  of  capacitances,  one 
of  the  two  states  would  be  lower  in  energy  and  the  other  state  could  be  reached  only 
through  thermal  excitation.  Both  states,  (1,8)  and  (2,7),  are  necessary  for  conduction. 

Using  the  values  of  capacitance  measured  for  the  single  dots  given  in  Table  3.1, 
we  can  plot  the  charging  energy  as  a  function  of  gate  voltage,  shown  in  Fig.  3-24. 
The  addition  of  an  electron  to  the  smaller  dot  does  not  coincide  with  the  addition 
to  the  larger  dot.  In  this  case,  we  expect  the  conductance  to  be  negligible  and  only 
thermally-activated.  A  close-up  of  the  transition  from  (1,5)— »(1,6)  —>(2,6)— >-(2,7) 
is  shown  in  Fig.  3-25.  At  the  charge  degeneracy  (1,6), (2,6)  current  can  only  flow 
if  the  intermediate  states  (1,5)  and  (2,5)  are  also  accessible.  At  this  gate  voltage. 
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Figure  3-22:  Charging  energy  as  a  function  of  gate  voltage  for  two  unequal  dots,  with 
Obg  —  5Cag,  CbT,  —  SCaE,  and  Cag  —  10  aF. 


Figure  3-23:  Close-up  of  charging  energy  as  a  function  of  gate  voltage  for  two  unequal 
dots. 


U(N^,N,,)  [e^llC\ 
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Figure  3-24:  Charging  energy  as  a  function  of  gate  voltage  for  two  unequal  dots  using 
capacitances  measured  for  separate  dots  given  in  Sec.  3.4. 
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Figure  3-25:  Close-up  of  charging  energy  as  a  function  of  gate  voltage  for  two  unequal 
dots  using  capacitances  measured  for  separate  dots  given  in  Sec.  3.4. 
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these  states  happen  to  be  degenerate.  They  are  not  lowest  in  energy,  but  still  lower 
than  (1,7)  and  (2,7),  which  could  also  provide  a  conduction  mechanism.  If  the  initial 
state  were  (1,6),  current  could  flow  through  the  cycle  (l,6)->(2,6)-)-(2,5)^(l,6)  or 

(1.6) -)-(l,5)-)'(2,5)-)-(l,6).  It  should  be  noted,  however,  that  conduction  is  more 
likely  at  the  degeneracy  point  of  (1,5)  and  (1,6)  since  the  energy  levels  (2,5)  and 

(2.6)  are  closer  at  this  point.  The  thermal  energy,  using  a  temperature  of  300  mK,  is 
about  0.1e^/2C.  Considering  the  scale  of  Fig.  3-25,  conduction  through  the  double 
dot  system  is  unlikely  at  this  temperature.  The  energy  diflierence  from  the  charge 
degeneracy  point  (1,6)  and  (2,6)  to  the  levels  (1,5)  and  (2,5)  is  approximately  5  A:bT. 
One  would  have  to  look  into  events  where  the  addition  of  an  electron  to  the  smaller 
dot  is  closer  in  gate  voltage  to  the  addition  of  an  electron  to  the  larger  dot.  Ideally, 
this  would  be  a  point  of  multiple  degeneracy,  such  as  the  one  shown  in  Fig.  3-23. 
If  the  ratio  between  the  gate  capacitances  is  a  rational  number,  such  an  event  can 
occur  regularly  as  a  function  of  gate  voltage.  If  the  ratio  is  not  a  rational  number,  we 
expect  to  observe  a  stochastic  Coulomb  blockade  effect,  as  predicted  by  Ruzin  [89] 
and  observed  by  Kemerink  [95]. 

Another  representation  of  the  conductance  condition  emerges  if  one  plots  the 
conditions  on  the  chemical  potential  for  the  independent  dots,  described  in  eq.  (3.60). 
We  neglect  all  electrodes  except  one  gate,  and  get  for  the  chemical  potential 

-  h  -  .  (3.70) 

If  we  consider  only  one  dot,  as  shown  in  Fig.  3-26,  we  can  plot  the  (negative)  chemical 
potential  versus  gate  voltage.  If  electrons  can  tunnel  into  and  out  of  the  dot,  the 
chemical  potential  as  a  function  of  gate  potential  is  a  sawtooth  function.  For  every 
discontinuity,  Na  is  increased  by  one.  If  we  plot  instead  —iJiaP^'^{Na—I/T)/CaT,  versus 
gate  voltage,  where  Na  increments  by  one  at  each  charge  degeneracy,  we  get  a  straight 
line  with  slope  &Cag/Caz  a  function  of  gate  voltage.  The  factor  CaglCaz  is  again 
the  capacitive  lever-arm.  This  straight  line  can  be  collapsed  onto  the  y-axis  without 
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loss  of  information.  The  gate  voltage  is  parallel  to  the  axis,  moving  with  a  speed 
of  —UVgCagICaZi  and  starting  offset  from  the  origin  by  e^/2Cai;.  The  conductance 
peaks  at  every  charge  degeneracy  point,  when  /To  =  0,  and  is  measurable  a  few  ksT 
around  these  points. 

For  two  dots,  this  graph  can  be  extended  to  include  a  second  dot  on  the  other 
axis,  as  shown  in  Fig.  3-27.  The  x-  and  y-axes  represent  a  measure  of  the  chemical 
potential  for  dots  a  and  b,  —iJ.a  +  e^/Ca^{Na  — 1/2)  and  — /Xft+e^/C'6s(-^fc  — 1/2),  where 
Na  and  iVf,  are  the  number  of  electrons  in  the  dots  in  equilibrium.  These  numbers  are 
changing  every  time  that  an  increase  in  gate  voltage  causes  a  charge  degeneracy  to 
happen,  represented  through  vertical  and  horizontal  lines.  The  gate  voltage  changes 
the  chemical  potential  in  the  dots,  indicated  by  the  line  in  the  plot.  The  capacitances 
between  the  gate  and  the  dots  cause  the  chemical  potentials  in  the  two  dots  to  sweep 
with  different  speed.  The  slope  of  the  line  is  given  by  ratio  of  the  capacitive  lever-arms 
of  the  gate  to  the  dots.  Conductance  peaks  appear  when  the  line  of  the  gate  voltage 
intersects  the  lines  of  charge  degeneracy  for  both  dots  simultaneously.  The  circles 
around  these  resonance  points  represent  the  thermal  spread  of  a  few  ksT  around  each 
charge  degeneracy  point,  as  given  by  eq.  (3.38).  Figure  3-27  is  a  picture  of  a  system 
with  Cbg  ^  4:0 ag  and  Cjs  >  Cax-  This  represents  approximately  the  conditions  for 
the  charging  energy  versus  gate  voltage,  plotted  in  Fig.  3-24.  The  ratio  between  the 
gate  capacitances  causes  about  four  electrons  to  be  added  to  the  larger  dot  before  one 
electron  can  be  added  to  the  smaller  dot.  A  few  of  the  electron  numbers  in  the  dots 
are  indicated  in  the  figure.  If  there  were  no  coupling  from  dot  b  to  the  gate  we  would 
obtain  the  situation  depicted  in  Fig.  3-26  if  conductance  through  the  inactive  dot  b 
is  possible,  or  no  conduction  at  all,  if  dot  b  is  blockaded.  A  close-up  of  the  chemical 
potentials  at  the  gate  voltage  where  the  second  electron  is  added  to  the  smaller  dot 
is  shown  in  Fig.  3-28.  In  this  close-up,  the  gate  voltage  range  is  chosen  similar  to 
Fig.  3-25,  and  the  three  charge-degeneracy  points  are  indicated  by  the  black  dots. 
The  double-dot  can  only  conduct  if  the  charge  degeneracy  point  is  close  to  a  common 
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Figure  3-27:  Plot  representing  change  in  chemical  potentials  in  both  dots.  Horizontal 
lines  represent  charge  degeneracy  in  dot  b,  vertical  lines  in  dot  a.  The  slope  of  the 
line  representing  the  change  of  chemical  potential  due  to  gate  voltage  change  depends 
on  the  ratio  of  the  capacitive  leverarms  of  the  dots.  Current  can  flow  if  four  adjacent 
charge  states  are  accessible  within  a  few  ksT. 
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Figure  3-28:  Close-up  of  Fig.  3-27  for  a  similar  gate  voltage  range  as  in  Fig.  3-25. 
Current  can  only  flow  if  all  states  adjacent  to  the  circles  are  accessible.  The  filled 
dots  represent  the  charge-degeneracy  points. 

charge  degeneracy  point  of  both  dots,  indicated  as  circles.  At  the  (1,5)-(1,6)  charge 
degeneracy  point,  the  system  can  conduct  only  if  the  location  of  the  degeneracy  is 
within  a  few  ksT  of  the  circle,  indicating  charge  degeneracy  of  (1,5),  (1,6),  (2,5),  and 
(2,6). 


3.6  Measurement  of  a  Double  Dot 

The  device  shown  in  Fig.  3-7  was  biased  in  a  configuration  shown  in  Fig.  3-14.  The 
conductances,  measured  as  a  function  of  gate  voltage  along  with  conductances  of 
the  individual  dots  for  comparison,  are  shown  in  Fig.  3-29.  A  modulation  of  two 
frequencies  can  be  observed  for  the  double  dot  case.  The  two  frequencies  appear 
similar  to  the  frequencies  observed  for  the  200  and  400  nm  quantum  dots.  An  analysis 
of  the  spacing  between  conductance  peaks  reveals  that  the  capacitances  are  larger 
than  for  the  original  dots.  The  capacitances  for  the  individual  dots  and  the  coupled 
dots,  extracted  from  this  plot,  are  shown  in  Table  3.2.  The  first  tall  peak  for  the 
double  dot  system  can  be  fit  using  eq.  (3.38),  similar  to  the  fits  for  the  individual 
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Figure  3-29:  Conductance  as  a  function  of  gate  voltage  for  separate  dots  of  200,  and 
400  nm  length  and  a  200  and  400  nm  dot  in  series.  For  comparison,  the  conductance 
for  the  600  nm  long  dot  is  plotted  as  well. 


200  nm  (a) 

400  nm  (b) 

individual  dots  [aF] 

2.0 

8.9 

coupled  dots  [aF] 

2.6 

11.8 

Table  3.2:  Capacitances  for  200  nm  dot  (a)  and  400  nm  dot  (b)  for  the  dots  measured 
individually  and  as  coupled  dots. 
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Figure  3-30:  Fit  of  the  first  conductance  peak  in  the  double  dot  case  with  the  Fermi 
function. 


Figure  3-31:  Nomenclature  for  the  electrodes  in  the  experiment. 

dots  in  Fig.  3-13.  The  result  is  shown  in  Fig.  3-30.  The  conductance  peak  fits  the 
same  temperature  and  lever-arm  as  the  400  nm  dot  alone. 

3.6.1  Threshold  Voltage  Shift 

The  threshold  voltage  for  the  onset  of  conduction  is  significantly  higher  for  the  double 
dot  system  than  for  individual  dots.  To  discuss  this,  we  adopt  the  nomenclature  shown 
in  Fig.  3-31  for  the  electrodes.  For  the  explanation  of  the  threshold  shift,  we  consider 
the  conductance  peaks  of  the  smaller  dot,  which  are  easier  to  identify  due  to  their 
sparsity.  The  gate  electrode  5  was  swept  in  all  experiments.  The  capacitance  Cab  sets 
the  distance  between  conductance  peaks.  If  the  potential  on  any  of  the  electrodes  i 
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Figure  3-32:  Shift  of  the  peak  associated  with  the  smaller  dot. 
is  changed  by  AVi,  each  conductance  peak  is  shifted  by 

AV,  =  ,  (3.71) 

^ag 

according  to  ec[.  (3.61).  The  shift  of  a  known  peak  reveals  the  capacitance  of  the  elec¬ 
trode  to  the  dot.  Figure  3-32  shows  the  first  strong  conductance  peak  for  three  dif¬ 
ferent  potentials  of  electrodes  1  and  2,  changed  by  10  mV  increments.  Using  eq.  3.71, 
the  capacitances  Cai  and  Ca2  are  approximately  3.4  aF.  Using  the  same  method,  we 
can  extract  capacitances  for  electrodes  7  and  8  to  be  Ca?  —  C'as  —  1-1  aF.  These  are 
the  electrodes  that  have  to  be  changed  if  one  wants  to  bias  the  large  dot  in  addition  to 
the  already  existing  small  dot.  To  create  a  third  tunneling  barrier,  gates  7  and  8  have 
to  be  swept  negative,  changing  the  gate  voltage  at  which  a  fixed  number  of  electrons 
will  be  in  the  dot.  This  is  the  reason  for  the  threshold  voltage  shift  in  Fig.  3-29. 
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Chapter  4 


Conclusion 


4.1  Summary 

This  report  has  explored  two  critical  aspects  of  future  electronics.  It  has  described 
fabrication  technology  for  nanometer-scale  semiconductor  devices  and  measurements 
of  novel  quantum  dot  devices  that  were  made  using  this  technology. 

In  order  to  make  devices  which  use  lateral  tunneling  structures,  lithography  has 
to  be  pushed  into  the  sub-lOOnm  region.  Lithographic  definition  of  tunneling  barriers 
has  to  be  on  the  order  of  a  few  tens  of  nanometers  and  control  over  these  dimensions 
has  to  be  even  more  stringent.  At  these  dimensions,  the  lithographic  task  to  build 
lateral  tunneling  devices  is  more  demanding  with  regards  to  resolution  than  other 
applications  of  x-ray  lithography,  such  as  MOS  gate  levels  or  optoelectronic  devices. 
These  requirements  leave  x-ray  and  e-beam  lithography  as  the  only  possible  tools 
for  the  definition  of  the  gate  level.  E-beam  lithography  is  a  process  featuring  high 
resolution  and  high  flexibility  in  pattern  formation.  However,  it  is  not  a  technique 
which  can  be  used  in  large-scale  fabrication  processes  because  of  its  serial  nature. 
Other  advantages  of  x-ray  lithography,  directed  especially  towards  the  engineering 
community,  include  superior  process  latitude,  printing  of  high  aspect  ratio  resist  lines, 
and  the  ability  to  write  large  area  patterns  on  x-ray  masks.  E-beam  lithography  can 
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not  easily  be  used  to  write  large  area  coherent  structures  because  of  the  stitching  error 
at  field  boundaries.  This  work  used  e-beam  lithography  only  to  lay  down  a  pattern  on 
a  parent  x-ray  mask,  which  was  then  replicated  once  onto  a  daughter  x-ray  mask  and 
then  printed  onto  a  substrate  using  x-ray  lithography.  In  the  course  of  this  work,  new 
technologies  were  developed  for  the  replication  and  processing  of  x-ray  masks.  Apart 
from  the  diffusive  wires,  made  by  John  Scott-Thomas,  which  initiated  the  discovery 
of  Coulomb  blockade  in  semiconductor  structures,  these  experiments  marked  the  first 
time  that  Coulomb  blockade  devices  were  fabricated  using  x-ray  lithography. 

A  unique  device  consisting  of  four  adjacent  quantum  point  contacts  was  built  and 
measured  at  cryogenic  temperatures.  In  the  course  of  this  work,  the  measurement 
electronics  and  data  acquisition  system  were  designed  and  built.  The  device  was 
measured  at  a  temperature  of  300  mK  in  a  commercially-available  probe.  The  design 
of  the  device  allowed  biasing  of  quantum  dots  of  various  sizes,  as  well  as  biasing  of  a 
double  dot  structure.  Using  two  adjacent  quantum  point  contacts  to  bias  the  device, 
the  center- to-center  distance  of  tunneling  barriers  was  200  nm  for  the  smallest  dot 
and  600  nm  for  the  largest.  Conductance  oscillations  as  a  function  of  gate  voltage 
were  seen  for  quantum  dots  of  all  sizes  and  the  capacitances  of  these  dots  were  ex¬ 
tracted.  The  gate-to-dot  capacitance  of  the  smallest  quantum  dot  was  of  the  order 
of  only  2  aF.  The  device  was  also  biased  to  produce  two  consecutive  quantum  dots 
of  400  and  200  nm  length,  as  measured  from  the  center  of  the  quantum  point  con¬ 
tacts  defining  the  dot.  A  beating  pattern  of  two  different  frequencies  was  observed 
for  the  conductance  oscillations  as  a  function  of  gate  voltage.  For  zero  temperature, 
quasi-random  conduction  resonances  would  be  expected.  Because  of  the  relatively 
high  temperature  during  the  measurements,  the  conductances  peaks  are  broadened, 
leading  to  the  observed  beating  pattern. 
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4.2  Future  Work 

This  report  has  shown  that  it  is  possible  to  make  quantum-effect  devices  using  tech¬ 
nology  that  is  suitable  for  mass  fabrication.  It  has  investigated  the  transition  from 
single  to  coupled  quantum  dots.  A  logical  continuation  is  the  investigation  of  more 
than  two  quantum  dots  in  close  proximity,  leading  eventually  to  large  arrays.  This 
effort  could  merge  previous  experiments,  which  employed  grid  gates  to  produce  peri¬ 
odic  potentials  within  a  2DEG  [7,8,14,96],  with  more  recent  work  on  coupled  quantum 
dots[40, 78,94].  So  far,  most  of  the  work  on  quantum  dot  devices  have  employed  only 
one  high-resoution  lithography  step,  the  gate  level.  With  increasing  device  density,  it 
will  be  necessary  to  integrate  other  lithography  steps,  most  likely  the  mesa  etch,  and 
possibly  the  ohmic  contact  level,  into  the  fabrication  process.  For  this  to  work,  high 
resolution  alignment  methods  have  to  be  used.  All  of  these  technologies  have  been 
developed  in  the  NanoStructures  Laboratory,  but  so  far  have  not  been  integrated  into 
one  whole  process.  Although  the  gate  level  will  most  likely  require  PMMA  as  resist  to 
ensure  high  resolution,  mesa  and  ohmic  contacts  could  employ  chemically-amplified 
resists  instead,  reducing  exposure  times  and  making  the  process  more  compatible 
with  manufacturing  environments.  Another  promising  area  is  to  try  to  diminish  dot 
capacitances  by  using  new  material  systems  such  as  SOI  which  are  specifically  de¬ 
signed  to  reduce  parasitic  capacitances  [83].  Coulomb  blockade  could  potentially  be 
observed  at  liquid  nitrogen  temperatures,  making  the  issue  of  single-electronics  much 
more  plausible  to  the  engineering  community.  This  effort  could  be  based  in  part  on 
the  expertise  of  the  SOI  research  groups  at  MIT. 
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CHAPTER  4.  CONCLUSION 


Appendix  A 


Mask  Fabrication  Process 


The  process  of  patterning,  developing,  and  electroplating  an  x-ray  mask  is  at  the 
moment  done  by  graduate  students.  Ultimately,  as  the  process  becomes  more  common 
and  ceases  to  be  an  area  of  research  interest,  this  task  should  be  transferred  to  a 
technician.  This  process  outline  describes  the  steps  for  the  fabrication  of  an  x-ray 
mask  by  a  technician  under  supervision  of  a  graduate  student. 

An  approximate  time  schedule  is  shown  in  Fig.  El-l.  Depicted  is  a  schedule  for 
the  fabrication  of  an  x-ray  mask  which  could  be  finished  within  about  one  working 
day.  Note  that  the  times  noted  are  estimates  in  case  no  major  problems  appear 
during  the  fabrication  sequence.  It  is  possible  to  cut  down  on  time,  mostly  during 
the  pump-down  for  the  RIE  and  SEM.  It  should  be  noted,  however,  that  because  of 
the  parallelism  of  the  processes  this  might  not  decrease  the  total  time  by  too  much  if 
only  one  person  is  involved  in  the  process.  Characterization  of  the  plating  bath  and 
mixing  of  the  developer  for  example  has  to  be  done  during  the  first  RIE  pump-down, 
no  matter  how  long  that  pump-down  takes. 

•  Turn  on  the  pump  for  the  plating  system  and  the  temperature  controller  for 
the  water  bath.  The  plating  system  needs  a  while  to  stabilize  to  the  desired 
bath  temperature  (33  °C).  The  set  point  for  the  water  bath  might  have  to  be 
adjusted  to  reach  the  desired  bath  temperature.  It  is  generally  around  29  °C. 
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Mask  Plating  Schedule 


RIE 

Plating 

Develop. 

Inspect. 

Superv. 


hours 


pump  down  etch  vent  &pump  down  etch  &  vent 


turn  on  characterize  bath 


characterize  bath  plate 


clean  up 


mix  dev.  develope 


strip  resist 


frontload 


Inspection 

in  optical  microscope 


plate  &  inspect 


inspect 


inspect 


Stud  evaporation  for  mother  mask:  ~2  hrs. 

Corrugation  grating  for  daughter  mask:  ~4  hrs. 

Removal  of  gold  particles:  ~2  hrs. 

X-ray  exposure  for  daughtering  of  mask:  7-30  hrs.,  avg.:  13  hrs. 

Figure  A-1:  Approximate  lab  schedule  for  the  mask  fabrication  process.  To  make 
a  working  daughter,  it  is  necessary  to  run  through  this  process  twice.  Extra  time 
must  be  devoted  for  the  evaporation  of  studs  for  the  mother  mask,  fabrication  of  the 
corrugation  grating  in  case  the  daughter  mask  is  to  be  used  in  intimate  contact  with 
the  substrate,  the  removal  of  gold  particles  by  wet  etching,  and  the  x-ray  exposure 
for  daughtering  of  the  mask. 
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•  Start  a  clean-up  run  on  the  RIE  to  set  up  the  parameters  for  an  oxygen  de-scum 
etch  on  the  Plasma-therm  RIE.  Etch  parameters  right  now  are  He  (10  seem) 
O2  (2.5  seem)  at  6  mTorr  for  12  sec.  Power  settings  are  44  W  at  75  V  DC.  He 
backside  cooling  is  set  to  20  °C  at  7  Torr  and  2.6  seem  He  flow.  This  step  should 
only  be  necessary  if  the  same  RIE  was  used  previously  with  different  gases. 

•  During  the  clean-up  run,  the  plating  bath  can  be  characterized.  This  involves 
measurements  of  pH-value  and  conductivity.  In  some  cases  an  IV-characteriza- 
tion  might  be  advisable.  Brightener  should  be  added  if  needed. 

•  Mix  the  developer  (1:2  MIBKrIPA)  in  a  1  liter  beaker.  The  desired  temperature 
for  the  developer  of  21  °C  can  be  achieved  by  heating  of  the  developer  on  a 
hotplate  or  warming  it  with  the  hands  on  the  outside  of  the  beaker.  The  latter 
method  avoids  an  overshoot  of  the  temperature.  In  case  the  temperature  is 
too  high,  the  developer  can  be  cooled  down  by  spraying  IPA  on  the  outside  of 
the  beaker  and  blowing  with  the  nitrogen  gun.  Monitor  the  temperature  of  the 
developer.  The  temperature  should  drop  to  about  16  °C  upon  mixing.  If  that 
doesn’t  happen  something  is  wrong  with  the  chemicals  used. 

•  Once  the  developer  temperature  has  settled  to  21  °C,  the  mask  is  developed  for 
90  seconds  by  slow  and  constant  agitation  in  the  developer.  The  development 
should  be  kept  the  same  for  all  masks. 

•  When  the  development  time  is  up,  immediately  rinse  the  mask  with  IPA  using 
the  squirt  bottle. 

•  Blow  the  sample  dry  with  the  filtered  N2  gun.  It  is  critical  not  to  blow  on  the 
sample  too  hard.  The  pressure  on  the  N2  gun  can  be  reduced.  Make  sure  to 
blow  dry  starting  in  the  middle  of  the  mask  and  work  your  way  outside. 

•  optically  inspect  the  mask  with  an  optical  microscope  in  clearfield  and  darkfield. 
This  should  be  done  by  a  person  familiar  with  the  pattern. 
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•  Load  the  sample  into  the  RIE  for  de-scum. 

•  Etch  for  approx.  5  seconds  using  the  parameters  from  the  clean  up  run  saved 
in  the  computer. 

•  Keep  the  sample  under  vacuum  in  the  RIE  until  you  are  ready  to  plate. 

•  Electroplate  200  nm  of  gold.  This  should  be  done  by  a  person  familiar  with  the 
pattern.  Plating  thickness  should  be  determined  using  optical  transmission,  but 
other  methods  such  as  alpha-step  can  be  used  as  well. 

•  Check  the  pattern  by  inspection  in  the  optical  microscope  and,  if  needed,  the 
scanning  electron  microscope.  The  inspection  should  be  performed  by  a  person 
familiar  with  the  pattern. 

•  Strip  the  resist  by  rinsing  in  acetone  and  methanol  and  plasma  ashing.  Alter¬ 
natively,  the  resist  can  be  stripped  by  soaking  the  mask  in  hot  NMP  (-^OO^C). 

•  Check  on  an  optical  flat  if  there  are  any  large  particles  on  the  mesa  that  prevents 
a  uniform  gap.  Alternatively,  the  mask  can  be  brought  in  close  proximity  with 
another  mask,  as  if  a  daughter  exposure  was  set  up.  The  process  for  removal 
of  the  particles  is  to  spin  photo  resist  over  the  whole  mask  and  to  expose  the 
suspicious  particles  for  about  2  minutes  in  the  optical  microscope  with  maximum 
white  light  intensity.  Then  use  an  eye  dropper  to  put  drops  of  gold  etch  on  the 
openings  in  the  resist.  The  etch  of  the  largest  particles  should  be  completed  in 
about  20  minutes. 

•  For  masks  to  be  used  in  microgap  exposures,  evaporate  Al-studs.  The  height 
of  the  studs  should  be  about  3  fim.  For  mask  replication,  the  Al-studs  are  not 
always  necessary  because  the  natural  warp  of  the  mesa  sets  an  initial  gap  of 
about  3  /rm. 
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•  Do  a  daughter  exposure.  It  is  typically  advisable  to  do  a  first  exposure  in 
head  1  with  a  gap  determined  by  the  studs.  If  the  effects  of  diffraction  prevent 
successful  replication,  another  exposure  can  be  made  in  head  3  using  the  vacuum 
fixture  and  a  gap  of  about  1  /xm.  The  gap  can  be  held  constant  using  electronic 
feedback  (see  Appendix  I). 

•  Flood  expose  the  daughter  mask  using  deep  UV.  A  TEM  grid  and  a  thin  wire 
should  be  used  to  monitor  development  and  plating.  Typical  exposure  time  is 
30  min. 

•  Process  the  daughter  mask.  The  processing  sequence  is  the  same  as  for  the 
mother  mask,  only  that  the  development  time  is  not  fixed  and  should  be  done 
by  a  person  familiar  with  the  pattern.  Also,  puddle  development  rather  than 
immersion  development  is  used  for  the  daughter  mask  because  of  improved 
particle  control.  The  improved  temperature  control  of  immersion  development 
is  not  necessary  for  daughter  masks  because  of  the  larger  process  latitude  offered 
by  x-ray  exposures  compared  to  e-beam  exposures. 

•  If  the  daughter  mask  is  to  be  used  for  exposures  in  intimate  contact  with  the 
substrate,  a  polyimide  corrugation  grating  is  necessary  for  outgassing  of  the 
resist.  For  these  steps  it  is  a  good  idea  to  process  a  silicon  monitor  in  parallel: 

-  Strip  PMMA  from  exposed  daughter  mask.  One  good  method  of  stripping 
is  to  let  the  mask  soak  in  hot  NMP  (QO^C)  for  15  min.  Be  careful:  NMP 
has  a  fieish  point  of  93  °C.  A  temperature  controlled  hot  plate  should  be 
used. 

-  If  the  plating  base  is  to  be  removed  in  the  alignment  mark  regions:  spin 
photoresist,  e.  g.  Shipley  1813. 

-  Expose  the  alignment  mark  regions  for  appropriate  time.  I  typically  use  the 
substage  condenser  and  expose  through  the  membrane  which  causes  the 
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thick  gold  areas  to  remain  covered  with  resist  (no  loss  of  contrast  during 
gold  etch!).  Typical  exposure  times  I  use  is  10  sec.  with  the  substage 
condenser  aperture  fully  open  and  the  Leitz  Xenon  light  as  the  light  source. 

Develop  with  Shipley  351  (1:5  351  :DI)  developer  for  about  30  sec. 

Plasma  ash  the  sample  in  He/Oa  for  6  seconds  at  about  50  W  power. 
This  step  is  absolutely  critical  to  get  etching  in  all  areas  at  the  same  time 
with  good  uniformity.  I  didn’t  get  any  consistent  results  without  this 
step.  This  step  is  probably  necessary  because  of  the  standing  wave  effect 
in  optical  lithography.  An  intensity  node  can  appear  slightly  above  the 
surface  of  the  sample  and  cause  a  very  fine  resist  layer  to  remain.  It  might 
be  possible  to  cut  down  on  the  exposure  time  though.  The  time  should  not 
be  increased,  however,  because  prolonged  exposure  to  the  plasma  crosslinks 
the  photoresist  which  makes  it  difficult  to  remove. 

Etch  the  gold  in  10:1  dilute  gold  etchant  for  20  seconds. 

Etch  the  NiCr  in  10:1  dilute  NiCr  etchant  for  5  seconds.  This  etch  is  much 
more  dangerous  than  the  gold  etch  because  it  tends  to  etch  laterally  at 
an  incredible  speed  and  undercuts  everything,  causing  Moire  alignment 
marks,  for  example,  to  fall  over  easily. 

Spin  about  300  nm  of  PI,  use  3:2  PI2610:NMP  at  3  krpm  for  60  sec. 

Clean  off  the  polyimide  at  the  edges  of  the  mask  with  a  soaked  Q-tip  for 
good  electrical  contact  with  the  mask  holder.  Remember:  we  use  electro¬ 
static  contact,  so  we  have  to  make  electrical  contact  to  the  mask. 

If  the  polyimide  should  be  stripped  in  the  alignment  mark  area,  bake  this 
layer  at  140  °C  like  the  following  layer.  Spin  photoresist,  develop,  strip, 
then  hardbake  the  polyimide.  If  polyimide  should  not  be  etched  away,  just 
bake  at  130  °C  and  250  "C  for  30  minutes  each. 

Spin  on  5:15  PI:NMP  to  get  roughly  70  nm  of  polyimide. 
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Bake  at  140  "C  for  30  min.  This  temperature  is  critical. 

Immediately  spin  on  photoresist  and  bake  for  30  minutes. 

Expose  with  a  100  /xm-period  grating  flex  mask,  and  if  necessary,  expose 
the  alignment  marks. 

Develop  for  30  seconds  using  5:1  DI  water:Shipley  351  developer.  Use  the 
monitor  before  you  do  your  mask. 

Examine  the  grating  in  the  optical  microscope. 

Strip  the  remaining  resist  using  acetone  and  methanol. 


Bake  at  200  °C  for  1  hour. 
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Appendix  B 


Layout  Transfer  from  MIT  to 
NRL 


After  the  device  layout  is  finished,  and  after  proper  alignment  marks  and  diagnostics 
patterns,  as  described  in  Appendices  F  and  E,  are  dropped  in,  the  KIC  file  needs  to 
be  converted  into  JEOL  readable  format  and  transferred  to  the  e-beam  site,  in  this 
case  NRL. 

•  If  the  KIC  layout  file  was  designed  on  a  machine  other  than  the  NeXT  machines 
transfer  the  file  via  ftp  to  the  NeXT  account.  Every  directory  that  is  used  for 
the  layout  contains  a  file  for  the  definition  of  the  layers,  called  ’.KIC’.  This 
file  should  be  sent  along  with  the  layout  file.  You’ll  need  it  for  the  kictocif 
conversion. 

•  Convert  the  file  from  KIC  format  to  CIF  format  using  the  command  kictocif. 
This  command  should  be  run  on  the  NeXT  machines.  For  some  reason  other 
computers  such  as  MTL  have  problems  handling  the  large  integers  used  for  x- 
ray  masks  (remember:  our  unit  length  A  is  usually  only  10  nm,  which  is  .01 
microns  per  lambda.  This  causes  large  integers  to  appear  in  the  design  of  a 
1x1  cm  die.). 
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Ftp  the  GIF  file  to  Bill  Chu’s  account  on  submicron.ll.mit.edu.  Ask  somebody 
for  the  password  if  you  don’t  know  it.  Put  your  file  in  a  subdirectory  with  your 
name  to  avoid  confusion  with  all  the  .KIC  and  setup. cnv  files  from  other  users. 

Quit  ftp  and  start  X-windows  on  the  machine  you  want  to  use  for  the  conversion. 

Open  the  display  for  submicron.il  by  typing  “xhost  +  submicron. 11” . 

Remote  login  into  submicron.ll.mit.edu,  selecting  X-windows  at  the  first  prompt 
and  then  typing  your  machine  name.  Then  go  to  your  directory. 

Run  the  conversion  program  by  typing  cnv  from  the  shell.  If  you  have  a 
setup. cnv  program  in  your  directory  already,  it  will  ask  you  if  you  want  to  use 
the  existing  one.  If  your  .KIC  file  didn’t  change  from  last  time  and  you  didn’t 
change  part  of  your  layer  to  single  pass  lines,  you  can  use  the  old  setup. cnv. 
Otherwise  you  have  to  edit  the  existing  one. 

Use  the  arrow  keys  to  go  from  one  item  to  the  next.  Use  the  space  bar  to 
change  a  selection.  If  you  edit  your  setup.cnv  file,  here  is  what  you  have  to  set 
for  each  separate  layer:  you  want  to  use  the  5th  lens  rather  than  the  4th  lens. 
The  field  size  should  be  80  /xm  instead  of  100  /xm.  If  you  use  single  pass  lines  in 
one  particular  layer,  line  rank  should  be  different  from  shot  rank.  The  origin  of 
the  field  should  be  set  manually  to  0  /xm/0  ^m.  This  is  the  number  you  want 
to  change  if  you  have  problems  concerning  field  boundaries  going  through  your 
fine  features.  After  you  are  done,  select  “Continue” . 

In  the  next  menu  you  can  change  the  location  of  the  DATA  directory,  where 
the  JEOL  files  are  going  to  be  placed.  In  case  the  user  wants  to  merge  different 
files  that  are  supposed  to  be  written  with  identical  beam  currents  and  different 
doses,  set  Geo  stop  to  “YES”.  If  this  feature  is  not  used  one  can  usually  just 
skip  this  menu  and  select  “Continue” .  The  computer  will  start  crunching  and, 
after  it  is  done,  will  spit  out  huge  files  in  a  subdirectory  called  /DATA.  If  the 


113 


Geo  stop  is  set,  the  data  won’t  be  written  right  into  the  DATA  directory  but 
into  separate  *.geo  files,  which  have  to  be  merged  using  the  merge  program. 

•  In  case  the  above  mentioned  feature  of  merging  different  layers  into  one  layer 
is  used  this  has  to  be  done  at  this  point.  From  the  command  line,  simply 
type  “merge”.  The  program  asks  for  the  layers  to  be  merged.  If  all  layers  are 
supposed  to  be  merged,  type  “*.geo”.  The  program  then  asks  for  the  shot  ranks 
of  each  of  the  layers  to  be  merged.  This  shot  rank  is  assigned  a  different  dose 
at  runtime. 

•  Change  to  this  directory  and  check  where  the  field  boundaries  go  through  your 
layer  by  typing  “xpreview  layerl.jSl  Iayer2.j51”.  Multiple  layers  can  be  pre¬ 
viewed,  in  this  case  two  layers  are  previewed. 

•  There  are  two  very  helpful  programs  on  submicron.il  which  help  estimate  the 
writing  time  for  the  e-beam.  The  first  is  called  stats_51,  a  program  that,  when 
called  as  in  “stats-51  layerl.jSl”  calculates  the  total  area  of  boxes  and  length 
of  lines  to  be  exposed  in  the  e-beam  system.  Those  numbers  should  be  written 
down  for  every  layer.  The  second  program  is  called  calc  and  is  invoked  from 
the  command  line  by  typing  “calc”.  In  this  program  one  can  input  the  area 
of  boxes  and  lengths  of  lines  and  receives  an  approximate  time  for  the  writing. 
The  program  demands  additional  information  as  to  which  lens  is  used  (5th), 
which  dose  (e.g  400  y^C/cm^  and  3  nC/cm),  and  which  current.  The  calculated 
times  should  be  written  down  for  every  layer. 

•  If  the  field  boundaries  cut  through  your  sensitive  features  you  should  go  back 
and  reconvert  your  file  for  a  different  selection  of  the  origin.  It  might  take  a 
couple  of  iterations. 

•  After  you  are  satisfied  with  the  field  boundaries,  run  the  conversion  program 
once  more  after  typing  “script”.  This  will  cause  the  whole  output  of  the  con- 
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version  program  to  be  piped  into  a  file  called  typescript.  After  the  conversion 
program  is  done  exit  the  script  and  move  typescript  to  the  data  directory  under 
the  name  “logfile” ,  or  “script-date” . 

•  Ftp  all  the  data  files,  the  logfile  and  the  setup. cnv  file  to  the  nrl  account  on  the 
NeXT  machines.  The  ftp  connection  should  be  made  binary  (type  “binary”).  If 
you  don’t  know  the  password  of  the  nrl  account,  ask  the  system  administrator. 
Again,  if  there  doesn’t  exist  a  subdirectory  with  your  name,  create  it. 

•  Remove  all  the  files  that  are  not  necessary  on  the  submicron.il  account.  The 
DATA  directory  sometimes  contains  huge  amounts  of  data  and  its  contents 
should  be  removed  if  at  all  possible.  If  you  want  to  keep  anything  in  your 
directory  for  a  while,  like  your  cif  file,  please  compress  it. 

•  Log  out  of  submicron  and  log  into  the  nrl  account  on  the  NeXT  machines.  Go 
into  your  subdirectory. 

•  Run  the  program  that  converts  the  .j51  format  into  an  equivalent  VAX-readable 
format.  NRL  reads  the  data  into  their  machine,  which  happens  to  be  a  VAX. 
The  program  can  be  started  from  your  subdirectory  by  typing  “../convert /layer 
>  junk”.  The  addition  “>  junk”  is  my  creation  because  I  wanted  to  suppress 
the  screen  output  which  dramatically  slows  down  the  conversion  program.  After 
all  my  conversion  is  done  I  delete  the  file.  The  resulting  files  should  have  an 
extension  “.vax”  and  should  have  the  same  size  as  the  corresponding  “.j51”  files. 

•  Please  don’t  forget  to  immediately  compress  all  the  “.j51”  files  since  they  are 
not  needed  any  longer.  The  “.vax”  files  should  be  compressed  or  deleted  once 
the  pattern  is  transferred  to  NRL. 

•  Send  email  to  somebody  at  NRL  so  they  can  pick  up  the  files  and  transfer  them 
to  their  machine. 
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•  print  out  a  copy  of  the  file  ”/NSL/Procedures/E-beam/KIC  >J51.ai”.  This 
table  needs  to  be  filled  out  with  the  layer  specifications  and  desired  beam  cur¬ 
rents  so  folks  at  NRL  get  an  idea  of  what  they  should  do  with  your  files.  Send  or 
fax  this  sheet  to  NRL.  Alternately,  a  document  can  be  written  and  transferred 
to  the  nrl  account  as  a  PostScript  file. 

The  following  is  a  compilation  of  notes  written  down  while  Bill  was  showing  us 
around  the  on  the  NRL  Vax. 

The  following  steps  can  be  done  only  by  persons  who  are  authorized  to  use  the 
account  on  NRL’s  Vax.  The  machine  name  is  mpfvax.nrl.navy.mil  (128.60.15.31),  the 
account  used  for  storing  the  data  is  called  “mitcad”.  This  machine  runs  VMS,  so 
everything  you  have  learned  about  UNIX  is  obsolete  on  this  account.  Translations 
of  most  common  commands:  “Is”— ^ “dir”;  “rm”— >“del”.  VMS  has  a  feature  that 
prevents  overwriting  of  old  files  by  assigning  a  version  number  to  each  file. 

•  Ftp  the  .vax  files  to  the  NRL  machine  (“ftp  mpfvax.nrl.navy.mil”,  user  name 
MITCAD).  Set  ftp  to  binary.  This  process  can  take  hours  because  the  files 
are  frequently  very  large  and  there  is  no  common  compression  scheme  between 
UNIX  and  VMS  to  the  best  of  my  knowledge.  I  have  the  strong  suspicion  that 
the  file  transfer  is  significantly  faster  if  initiated  from  the  remote  machine. 

•  Log  into  mpfvax.nrl.navy.mil  into  the  MITCAD  account. 

•  Verify  that  the  files  you  just  sent  were  received  by  typing  “DIR  *.VAX”. 

•  If  there  is  a  previous  version  of  the  file  to  be  written,  delete  this  file  by  typing 
something  like  “del  NAME.vax:!” 

•  The  .vax  files  are  blocked  off  in  512  byte  blocks.  They  need  to  be  blocked  off 
in  blocks  of  4096  bytes.  This  is  accomplished  by  running  the  “F3”  program. 
First  check  to  see  what  the  largest-number  JEOL51*.DAT  file  is.  To  do  this, 
type  “dir  jeol51*”.  For  example,  if  the  largest  number  is  JEOL51080.DAT;1, 
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then  the  next  data  will  be  converted  into  JEOL51081.DAT;1.  To  initiate  this 
conversion,  type  “RUN  F3”.  This  program  will  ask  you  for  the  input  file.  Type 
“<filename>.VAX”.  Then  it  will  ask  you  for  the  output  file.  If  the  last  unused 
JEOL51  file  is  JEOL51080.DAT,  then  type  “JEOL51081.DAT”.  Repeat  this  for 
each  file.  For  the  next  file,  put  the  output  in  JEOL51082.DAT,  and  so  on. 

•  To  “fix”  the  directory,  run  the  fixing  program  by  typing  “RUN  J51DIRFIX”. 

•  Verify  that  the  new  files  have  been  updated  by  typing  “TYPE  JEOL51DIR- 
.DAT”. 

•  To  figure  out  the  chip-size  of  each  layer,  run  the  program  PATMNG.  Give  in 
“1”  (for  list-full)  and  “2”  (for  Disk  No.  :DUA1:)  and  answer  the  questions. 
This  information  was  given  to  me  by  Kee  Rhee  on  2-16-94. 

After  all  this  is  done  the  operator  at  NRL  should  be  able  to  easily  preview  the 


pattern. 


Appendix  C 
Electroplating 


Power  Supply 

In  order  to  electroplate  x-ray  masks,  it  is  important  to  use  a  well-controlled  and 
reliable  plating  setup.  One  of  the  essential  parts  of  the  setup  is  a  plating  current 
source.  This  appendix  describes  the  design  of  the  circuit. 

There  are  two  principal  ways  to  electroplate  a  metal  onto  a  substrate,  which  are 
referred  to  as  galvanostatic  and  potentiostatic  electroplating.  During  galvanostatic 
plating,  the  current  density  onto  the  substrate  is  held  constant,  while  for  potentio¬ 
static  plating,  a  voltage,  e.g.  the  voltage  between  cathode  and  a  saturated  calomel 
electrode  is  held  constant.  This  power  supply  was  designed  to  feature  both  modes  of 
operation. 

Figure  C-1  shows  a  schematic  of  the  plating  power  supply.  It  was  built  for  battery 
operation,  but  it  can  easily  be  converted  for  use  with  an  outside  power  supply.  The 
plating  bath  appears  on  the  top  left  of  the  diagram,  with  pins  12  and  14  connected  to 
the  cathode  and  anode,  respectively.  The  cathode  is  almost  grounded,  with  less  than 
10  to  ground,  due  to  the  simple  current  divider  into  the  current  amplifier,  which 
measures  the  plating  current  with  lOmV/mA  amplification.  The  anode  is  driven  by 
a  2N3904  transistor  because  the  OP97  used  to  drive  the  transistor  does  not  supply  an 
output  current  above  about  30  mA.  The  OP97  regulates  either  the  saturated  calomel 
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Simplified  Circuit  Diagram  of  Potentio/Galvanostat 

Q)  2N3904J 


Anode_ 


>1 


Bath  3Q^  Potentiostat  ^ 


Cathode 


battery 

indicator 


V. 
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Figure  C-1:  Schematic  circuit  diagram  of  the  plating  power  supply. 
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electrode  (potentiostat),  or  the  voltage  measured  with  the  current  amplifier  to  be 
equal  to  the  reference  voltage  which  is  set  using  a  6.8  Zener  diode  and  a  10-turn 
potentiometer.  The  circuit  on  the  bottom  consists  of  an  integrator,  which  integrates 
the  current  over  time,  J  Idt,  and  a  successive  amplifier.  If  the  ten  turn  potentiometer 
between  pins  7  and  10  is  set  to  be  proportional  to  the  plated  area  and  the  plating 
efficiency  is  unity,  the  voltage  on  pin  8  gives  the  plated  thickness,  in  A/mV.  The 
plating  efficiency  indicates  how  efficiently  the  ionic  current  in  the  solution  transfers 
metal  onto  the  sample.  A  plating  efficiency  of  unity  means  that,  for  a  unit  charge 
injected  into  the  solution  at  the  cathode,  one  gold  atom  is  plated  out  of  the  solution. 

For  trouble  shooting,  please  use  the  schematic  of  Fig.  C-1.  All  elements  used  in 
this  circuit  are  described  in  basic  circuit  books  [97,98]. 

There  are  several  points  that  could  be  improved  on  this  box: 

•  There  seems  to  exist  a  high  current  latch-up  state.  It  might  be  possible  to  avoid 
this  state  if  the  inputs  of  the  OP97  are  separately  buffered. 

•  Change  from  battery  operation  to  an  external  DC  power  supply.  Those  power 
supplies  are  commercially  available,  and  it  shouldn’t  take  more  than  an  hour  to 
change  the  box  accordingly,  once  such  a  power  supply  is  obtained. 

•  It  would  be  nice  if  the  plating  current  and  the  plated  thickness  could  be  read  out 
directly.  It  takes  ordering  of  parts,  in  this  case  two  miniature  digital  voltmeters, 
to  complete  the  job. 

Note:  to  use  the  power  supply,  one  should  turn  the  power  supply  on  only  if  the 
sample  is  disconnected  and  the  current  setting  potentiometer  is  zeroed.  This  will 
avoid  above  latch-up  condition.  Also,  the  circuit  will  latch  if  the  current  drawn  is 
more  than  r\j  45  mA. 
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Figure  C-2:  Qualitative  description  of  the  SCE  voltage  as  a  function  of  time  and 
plating  current  density. 

SCE  voltage  during  plating 

The  voltage  of  the  saturated  Calomel  electrode  should  be  monitored  during  the  plating 
at  all  time  using  a  strip  chart  recorder.  This  way  it  is  easier  to  detect  possible 
problems  if  the  voltage  deviates  from  previous,  similar  runs.  A  qualitative  desciption 
of  the  saturated  calomel  electrode  voltage  as  a  function  of  time  and  plating  current 
is  given  in  Fig.  C-2.  After  turning  on  the  current,  the  voltage  initially  decays  from 
the  starting  voltage,  and  starts  rising  a  again  after  about  two  minutes.  The  rise  in 
the  end  is  probably  due  to  depletion  of  the  solution  close  to  the  cathode. 

Design  Considerations  for  Plating  Fixture 

A  description  of  the  plating  apparatus  used  in  plating  masks  was  already  given  in 
section  2.5.1.  Here  are  some  comments  about  the  design  that  should  be  incorporated 
in  future  fixtures  and  tanks. 
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It  would  be  a  good  idea  to  construct  the  tank  out  of  plexiglass  or  PVC,  materials 
recommended  by  SEL-REX.  Glass  or  pyrex  beakers,  as  used  in  my  experiments,  are 
not  recommended  by  the  manufacturer,  although  it  is  not  on  the  list  of  materials 
to  avoid  either.  It  might  be  a  good  idea  to  avoid  these  materials  because  of  ionic 
impurities  that  might  be  present.  A  new  tank  should  also  provide  space  for  the 
access  of  at  least  two  sensors,  the  saturated  Calomel  electrode  and  a  thermometer.  I 
believe  that  it  is  a  good  idea  to  minimize  the  amount  of  plating  fluid  to  reduce  costs 
in  case  of  a  bath  replacement  as  well  as  to  reduce  temperature  stabilization  times. 

Notes  on  the  bath 

The  following  are  the  manufacturers  recommendations  for  Operation  of  the  bath  in 
an  abbreviated  form: 

•  pH  value  has  to  be  maintained  between  8.5  and  9.0.  This  is,  in  my  experience, 
rarely  a  problem,  but  it  can  be  regulated  by  addition  of  sulfuric  acid  or  sodium 
hydroxide. 

•  dull  deposit  add  2  ml/1  brightener. 

•  if  that  doesn’t  work  add  conducting  salts,  ~7g/l. 

•  if  that  doesn’t  work  — >■  more  serious  problem,  possibly  organic  contamination. 

•  a  sample  should  be  periodically  sent  to  be  analyzed  for  gold  content  and  con¬ 
ducting  salts. 
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Appendix  D 


Fabrication  Sequence 


Following  are  the  lithography  steps  used  to  make  the  devices  described  in  this  report. 
All  steps  are  done  with  Microposit  1813  as  resist,  except  for  the  gate  metallization, 
where  PMMA  is  used. 

•  Solvent  clean  the  sample. 

•  Spin  Microposit  1813  resist  -  I  use  4.3  krpm  for  30  sec  ~  1  /im. 

•  Bake  in  oven  at  90°  C  for  30  minutes. 

•  Perform  alignment  in  OAI  and  expose  for  appropriate  time,  typically  about 
4  sec. 

•  Develop  in  Microposit  351  developer  at  room  temperature. 

•  Perform  mesa  isolation  of  devices.  I  use  1000:20:6  H20:NH40H:H202  for  22 
sec  as  the  etch  and  typically  get  a  65  nm  mesa  (actually,  sometimes  I  get  up  to 
120  nm). 

•  Do  exposure  as  above  using  ohmic  contact  mask  level. 

•  UV  ozone  for  15  s. 
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•  Rinse  in  cold  DI  water  for  1  min. 

•  Rinse  in  Semico  for  1  min.  Semico  is  a  buffered  ammonium  hydroxide  solution 
that  is  supposed  to  strip  any  native  oxide  on  the  GaAs  surface. 

•  Rinse  in  cold  DI  water  for  1  min. 

•  Blow  dry. 

•  E-beam  evaporate  ohmic  contacts:  Ni(2):Au(4):Ge(22):Au(44)Ni(10):Au(30). 

•  Lift  off  and  anneal  sample  at  430°C  for  15  s.  in  RTA. 

•  Solvent  clean  the  sample. 

•  Spin  KTI  PMMA  950  K  4%,  5800  rpm  for  ~250nm.  I  lift  off  up  to  120  nm  with 
this  thickness. 

•  Bake  PMMA  at  180°C  for  1  hour. 

•  Expose  sample  for  a  dose  of  approximately  600  J/cm®,  as  determined  by  the 
x-ray  database  program,  Appendix  G. 

•  Develop  sample  for  5  seconds  in  standard  developer  and  check  development  in 
alpha  step.  The  developer  is  2:3  MIBK:IPA  and  the  development  is  performed 
using  a  squirt  bottle.  A  total  development  time  of  about  30  to  60  seconds  is 
desirable.  Check  development  also  in  the  microscope.  The  gas  relief  grating 
(see  Appendix  A)  should  be  visible  at  this  point,  and  in  the  areas  of  thicker 
polyimide  the  development  should  be  a  little  reduced. 

•  Continue  developing  until  the  development  is  complete  and  add  10  seconds. 
Inspect  again  in  optical  microscope  and  check  that  the  gas  relief  grating  is  not 
visible  any  more. 
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•  UV  ozone  clean  sample  for  30  seconds  to  remove  possible  scum  layer.  Failure 
to  do  so  might  result  in  electrically  unstable  gates. 

•  Before  evaporation  of  the  gate  metal,  the  native  oxide  on  the  surface  of  the 
semiconductor  has  to  be  removed.  The  sample  is  dipped  for  one  minute  into 
DI  water,  then  one  minute  into  Semico  solution,  then  one  minute  again  into 
DI  water.  As  an  alternative  to  Semico  clean,  a  mixture  of  H20:NH40H  can  be 
used. 

•  Use  thermal  evaporator  to  evaporate  80  to  100  nm  of  gate  metal.  One  can 
either  use  Aluminum  or  Gold/Palladium.  Use  of  an  electron  beam  evaporator 
might  result  in  damage  to  the  two  dimensional  electron  gas.  There  are  also  a  lot 
of  successful  device  results  made  with  e-beam  evaporators.  For  electron  beam 
evaporator,  I  usually  use  evaporation  of  10  nm  titanium  and  90  nm  gold. 

•  Perform  liftoff  by  heating  acetone  on  hotplate  for  30  sec.,  putting  sample  in  for 
1  min.  and  then  squirting  sample  with  acetone  from  the  plastic  squirt  bottle. 
Caution,  heating  acetone  is  like  heating  gasoline.  For  difficult  liftoff,  ultrasound 
may  be  used.  Another  solvent  that  can  be  used  for  liftoff  is  NMP,  which  is  also 
a  little  less  dangerous  when  heated.  (Arvind’s  recipe,  I  believe  he  learned  it 
at  IBM:  perform  liftoff  by  soaking  sample  for  30  minutes  in  a  mixture  of  1:1 
acetone:dichloro-methane.)  The  devices  are  now  finished. 

•  Cleaving:  put  the  piece  down  onto  a  paper  towel  and  put  a  glass  slide  on  each 
side.  Put  a  third  glass  slide  across  at  the  position  where  you  want  to  cleave  the 
sample.  It  might  be  necessary  to  use  a  stereo  microscope  if  you  want  to  cleave 
between  devices. 

•  Scribe  the  sample  close  to  the  bottom  edge  using  a  carbon  scribe.  It  is  not 
necessary  to  scribe  along  the  whole  length  of  the  cleaving  line.  GaAs  is  very 
brittle  and  will  easily  cleave  along  the  crystallographic  planes. 
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•  Position  the  sample  on  the  glass  slide  with  the  cleaving  line  along  the  edge  of 
the  glass  slide.  Exert  a  force  on  the  outermost  sample  edges.  The  sample  should 
break  without  the  exertion  of  much  force  if  the  lever  arm  is  large  enough.  The 
smaller  the  sample,  the  more  difficult  it  is  to  cleave. 

•  Take  a  44- pin  header  and  place  a  small  amount  of  silver  paint  on  it.  Be  careful 
not  to  put  too  much  onto  the  header  as  the  silver  paint  might  cover  also  the 
front  side  of  your  sample.  Press  the  sample  down  until  most  of  the  silver  paint 
underneath  the  sample  has  moved  outside. 

•  My  experience  with  bonding  devices  with  the  Kulicke  &  Soffa  bonder  of  the  ICL 
is  that  the  settings  of  three  for  time,  force,  and  power  is  a  good  bet  to  start. 
Much  higher  settings  present  a  physical  danger  to  the  device,  lower  settings 
won’t  complete  the  bond. 

•  Mount  the  sample  onto  the  bonder  chuck.  Make  sure  you  bond  only  the  wires 
which  have  been  grounded  first  to  the  chuck.  The  chuck  and  the  person  doing 
the  bonding  should  both  be  connected  to  ground  via  a  cable.  The  grounding  is 
especially  important  for  devices  made  with  e-beam  lithography  since  the  metal 
tends  to  be  thinner  than  for  devices  made  with  x-ray  lithography. 

•  For  the  first  bond,  position  the  bonder  wedge  appropriately  considering  the 
height  of  your  header  and  sample.  Position  the  mouse  over  the  first  header 
pad  and  press  the  mouse  button.  Adjust  the  search  button  to  place  the  wedge 
slightly  above  the  pad.  Then  go  to  the  sample  pad  and  place  the  wedge  slightly 
above  the  pad  as  well. 

•  For  the  remaining  bonds  you  can  use  the  z-stage  lever  with  your  right  hand  and 
position  the  sample  with  the  mouse. 

•  All  the  cables  in  the  measurement  setup  are  designed  as  twisted  pair  cables 
starting  with  the  pair  1-2.  If  at  all  possible  try  to  bond  source  and  drain  on 
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successive  bonds,  following  this  convention.  This  will  reduce  noise  pickup  due 
to  magnetic  fields  during  measurements. 
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Appendix  E 


Diagnostic  Pattern  for  E-Beam 
Lithography 


The  diagnostic  pattern  was  developed  to  be  dropped  into  existing  patterns  that  are 
written  with  the  electron  beam  on  x-ray  masks.  It  can  also  be  used  by  itself  as  a 
tool  to  evaluate  the  performance  of  the  electron  beam.  The  diagnostic  pattern  was 
developed  with  the  help  of  Scott  Silverman  and  Scott  Hector. 

The  layout  of  the  whole  pattern  is  shown  in  Fig.  E-1.  The  pattern  consists  of 
six  identical  layers  (indicated  as  different  grey  levels)  which  can  be  assigned  separate 
doses  or  currents  that  the  operator  wants  to  get  information  about.  The  different 
layers  are  called  D1-D6  in  the  layout  program.  The  pattern  is  designed  to  fit  within 
one  field  of  80  /xm  x  80  /xm  commonly  used  in  the  e-beam  system  at  NRL.  Under 
normal  circumstances,  it  should  be  possible  to  write  the  whole  pattern  within  about 
two  minutes.  The  individual  patterns  are  described  in  the  following  paragraphs: 


Single  pass  lines  and  star-burst  pattern:  The  single  pass  lines  shown  in  Fig.  E- 
2  and  as  “1”  in  Figure  E-1  were  designed  to  give  information  about  linewidth  of  single 
pass  lines  both  free  standing  and  in  proximity  with  other  single  pass  lines.  The  right 
angle  in  the  lines  insures  that  any  astigmatism  will  be  detected  as  the  difference  in 
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Figure  E!-l:  Layout  of  the  diagnostic  pattern, 
linewidth  for  x-  and  y-direction. 

The  so-called  star-burst  pattern  was  designed  to  write  lines  at  angles  of  30  "  and 
60"  in  addition  to  lines  in  x-  and  y-direction.  The  intention  was  to  show  astigmatism 
by  looking  at  the  superposition  of  the  lines.  In  reality,  however,  the  lines  at  30"  and 
60"  receive  slightly  different  doses  than  the  x-  and  y-directed  lines  because  of  the 
digitization  of  the  beam  into  squares.  The  dose  for  x-  and  y-directed  lines  tend  to  be 
higher.  Since  this  pattern  does  not  take  very  long  to  write,  we  believe  the  writing  of 
the  star-burst  pattern  is  still  worth  the  while. 


Boxes  and  lines:  This  pattern,  shown  in  Fig.  E-3  and  as  “2”  in  Fig.  E-1,  was 
designed  to  print  square  boxes  and  lines  using  areal  doses.  The  nominal  linewidths 


Figure  E-2;  Single  pass  lines  and  star-burst  pattern  in  the  diagnostic  pattern. 


and  side  lengths  of  the  squares  are  20  nm  to  200  nm  in  20  nm  increments,  plus  300 
and  400  nm.  The  square  boxes  are  of  special  significance,  for  example  as  contact  holes 
in  silicon  processing. 

Lines,  grids,  and  inverted  boxes:  This  pattern,  shown  in  Fig.  E-4  and  as  “3”  in 
Fig.  E-1  was  designed  to  print  boxes  of  inverted  polarity  and  single  pass  lines  adjacent 
to  large  areas,  and  to  write  grids  and  gratings  of  varying  periodicity.  The  inverted 
polarity  boxes  have  the  same  sizes  as  the  boxes  in  pattern  “2” ,  but  are  of  the  opposite 
polarity. 

Inverted  Line:  This  pattern,  shown  in  Fig.  E-5  and  as  “4”  in  Fig.  E-1  was  de¬ 
signed  to  check  the  printing  of  an  inverted  line  in  the  presence  of  large  written  areas 
surrounding  the  line.  To  check  for  stigmation,  there  is  one  pattern  in  a:-direction  and 
one  in  y-direction.  The  linewidth  ranges  from  20  nm  to  400  nm  in  20  nm  increments. 
The  first  and  last  segments  are  500nm-long,  all  other  segments  are  200nm-long.  The 
width  of  the  boxes  on  the  sides  of  the  line  is  1  /^m. 
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I  I  I  I  I  I 


Figure  E-3:  Boxes  and  lines  in  the  diagnostic  pattern. 

Effect  of  proximity  effect  on  a  single  line:  This  pattern,  shown  in  Fig.  E-6 
and  as  “5”  in  Fig.  E-1  tests  the  proximity  effect  of  adjacent  boxes  on  a  single  pass 
line.  The  space  between  the  boxes  is  held  constant  at  400  nm,  the  size  of  the  boxes  is 
varied  from  100  nm,  200  nm,  400  nm,  800  nm  up  to  1.6 /rm.  To  check  for  stigmation, 
there  is  one  pattern  in  x-direction  and  one  in  y-direction. 
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Figure  E-6:  Proximity  effect  test  structure  in  the  diagnostic  pattern. 


Appendix  F 


Kic  Layout  Tools 


This  chapter  describes  the  layout  tools  that  were  developed  in  the  course  of  the  report 
research.  The  layout  program  used  for  design  of  devices  is  KIC.  KIC  is  a  relatively 
primitive  editor  which  produces  a  file  which  consists  of  a  short  header,  layer  names 
for  each  new  layer  and  a  list  of  boxes,  wires  and  polygons. 

The  command  for  a  typical  box  is: 

B  390000  5000  605000  390000; 

This  line  describes  a  box  of  width  3900  A  and  height  50  A  and  center  coordinates 

of  x:6050  A  and  y:3900  A. 

The  command  for  a  typical  single  pass  line  is: 

WOO  6400000  30000000  6400000; 

This  line  describes  a  line  of  zero  width  (single  pass)  and  the  x-  and  y-coordinates 
of  the  beginning  and  end  points  of  the  line. 

Moire  Alignment  Mark 

For  high  accuracy  alignment  of  x-ray  masks,  the  alignment  marks  are  of  great  impor¬ 
tance.  This  section  describes  the  code  used  to  generate  alignment  marks  as  described 
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in  ref.  [99].  The  alignment  mark  consists  of  gratings  of  two  different  periodicities 
which  are  perpendicular  to  two  others.  Such  an  alignment  mark  is  shown  in  Fig.  F-1. 
When  aligning  a  mask  to  a  substrate,  a  pair  of  two  horizontal  and  vertical  fringe- 
patterns  each  appear.  The  fringes  of  each  pair  move  in  opposite  directions  when  the 
substrate  is  moved  with  respect  to  the  mask. 

The  program  to  produce  the  alignment  mark  is  listed  below: 

#iiiclude  <stdio.h> 

#include  <math.h> 

/♦this  program  creates  an  alignment  mark  for  fine  alignment  of  x-ray  */ 
/♦masks  .  It  expects  an  input  file  containing  the  periodicity  of  the  ♦/ 
/♦gratings.  The  input  file  should  contain  "periodl,period2,size."  of  ♦/ 
/♦the  alignment  mark.  It  produces  triangular  gratings  with  aspect  ♦/ 
/♦ratio  1:1.  Lambda  is  going  to  be  10  nm  and  the  gratings  are  to  be  a  ♦/ 
/♦triangles  of  base  80  urn  and  40  urn  height.  The  horizontal  grating  ♦/ 

/♦is  going  to  be  of  period  1,  the  vertical  of  period2.  I  am  not  quite  ♦/ 
/♦sure  if  I  have  to  account  for  the  rounding  error  occuring  from  the  ♦/ 
/♦conversion  from  float  to  integer  for  nycl  and  nxc2  ♦/ 

/♦written  by  Martin  Burkhardt  l-94^/ 

mainO 

{ 

float  xcl,wl,ycl,hl;  /♦width,  height,  and  center  coordinates  of  boxes^/ 

float  xc2,w2,yc2,h2; 

int  nycl,nhl,nxcl,nwl; 

int  nyc2,nh2,nxc2,nw2; 

int  i,nol,no2,nwidth; 

int  offset; 

double  periodl,period2, width; 

FILE  ♦fp,  ♦input,  ♦fopenO; 
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fp  =  f open(" align. kic" , "w") ; 
input  =  fopen("input.dat","r") : 

/*the  following  input  is  expected  from  the  input  file  "input .dat"*/ 

fscanf  (input, "y,lf  7.1f  %lf " ,&periodl,&period2,&width) ; 
nol  =  width/ (2*periodl) ; 
no2  =  width/ (2*period2) ; 

periodl  =  periodl  *  10000;  /*convert  microns  into  lambda*/ 

period2  =  period2  *  10000; 

width  =  width  *  10000; 

nwidth  =  width;  /^integer  conversion*/ 

offset  =  nwidth  /  2; 

fprintf (fp,"(S3raibol  align. kic) ;\n") ; 
fprintf (fp, "9  align.kic;\n") ; 
fprintf (fp,"DS  0  1  l;\n"); 
fprintf (fp,"L  HRA;\n"); 

/*central  lines*/ 

xcl  =  width/4;  /*x-coordinate  of  center  1*/ 
xc2  =  0;  /*x-coordinate  of  center  2*/ 
nxcl  =  xcl; 
nxc2  =  xc2; 

ycl  =  0;  /*y-coordinate  of  center  1*/ 
yc2  =  width/4;  /*y-coordinate  of  center  2*/ 

nycl  =  ycl;  /*might  have  to  add  half  a  unit  to  account  for  rounding*/ 
nyc2  =  yc2; 
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hi  =  periodl  /  2;  /*height  of  box  1*/ 
h2  =  width  /  2;  /^height  of  box  2*/ 
nhl  =  hi; 
nh2  =  h2; 

wl  =  width  /  2;  /*width  of  box  1*/ 
w2  =  period2  /  2;  /*width  of  box  2*/ 
nwl  =  wl ; 

Iiw2  =  w2; 

fprintf  (fp,"B  '/.d  7.d  7,d  7.d;\n"  ,nwl  ,nhl .nxcl+off  set ,nycl+of f  set) ; 
fprintf  (fp,"B  7.d  7.d  7.d  7.d:\n" , nwl  ,nw2, -nxcl+off set .nycl+off set)  ; 
fprintf  (fp,"B  7.d  7.d  7.d  7.d;\n"  ,nw2,nh2,nxc2+off  set,nyc2+offset)  ; 
fprintf  (fp."B  7.d  7.d  7.d  7.d;\n"  ,nhl  ,nh2,nxc2+off set,-nyc2+off set)  ; 

/*all  other  lines*/ 


f or (i=l ; i<=nol ; i++){ 

xcl  =  width/4  +  i  *  periodl  /  2;  /*x-coordinate  of  center  1*/ 
ycl  =  i  *  periodl;  /*y-coordinate  of  center  1*/ 
nxcl  =  xcl; 

nycl  =  ycl;  /*might  have  to  add  half  a  unit  to  account  for  rounding*/ 

hi  =  periodl  /  2;  /*height  of  box  1*/ 

wl  =  width  /  2  -  i  *  periodl;  /*width  of  box  1*/ 

nhl  =  hi; 

nwl  =  wl; 

fprintf  (fp,"B  7.d  7.d  7.d  7.d;\n" , nwl , nhl, nxcl+off  set , nycl+off  set)  ; 
fprintf  (fp,"B  7.d  7.d  7.d  7.d  ;\n", nwl, nhl. nxcl+off  set, -nycl+off  set)  ; 
fprintf  (fp,"B  7.d  7.d  7.d  7.d;  \n", nhl, nwl, nycl+off  set, -nxcl+off  set) ; 
fprintf  (fp,"B  7.d  7.d  7od  7.d;\n",nhl,nwl,-nycl+offset,-nxcl+offset) ; 

> 

f or(i=l ; i<=no2; i++){ 
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yc2  =  width/4  +  i  *  period2  /  2;  /*x-coordinate  of  center  2*/ 
xc2  =  i  *  period2;  /*x-coordinate  of  center  2*/ 
nyc2  =  yc2; 
nxc2  =  xc2; 

w2  =  period2  /  2;  /^height  of  box  2*/ 

h2  =  width  /  2  -  i  *  period2;  /^height  of  box  2*/ 

nw2  =  w2; 

nh2  =  h2; 

fprintf  (fp,"B  7.d  '/.d  7.d  7.d;\n" .nw2,nh2,nxc2+offset,nyc2+offset) ; 
fprintf  (fp,"B  7.d  7,d  7.d  7.d;\n",nw2,nh2,-nxc2+offset,nyc2+offset) ; 
fprintf  (fp,  "B  7.d  7.d  7.d  7.d;\n"  ,nh2,nw2,-nyc2+offset,nxc2+offset)  ; 
fprintf  (fp,"B  7.d  7.d  7.d  7.d;\n"  ,nh2,nw2,-nyc2+offset,-nxc2+offset) ; 

} 

fprintf (f p , "DF ; \n" ) ; 
fprintf (fp, "E\n") ; 

} 

} 


Appendix  G 


X-Ray  Lithography  Run-sheet 


This  appendix  describes  the  X-ray  run-sheet  developed  for  the  electron  bombardment 
x-ray  sources  in  the  Nanostructures  Laboratory  at  MIT.  The  run-sheet  uses  a  database 
file  that  runs  off  a  NeXT  computer  in  the  lab.  An  example  of  the  run-sheet  is  shown 
in  Fig.  G-1.  The  advantage  of  running  the  run-sheet  off  a  network  computer  is  that  all 
the  parameters  for  the  run  can  be  accessed  from  remote  locations.  Another  advantage 
is  that  the  parameters  given  for  the  run  can  be  converted  into  doses  with  the  use  of 
corresponding  formulas,  and  development  times  can  be  estimated.  In  this  database, 
the  dose  incident  and  dissipated  in  the  resist  is  calculated,  the  development  time 
for  the  samples  is  recorded,  and  values  for  diffraction  coefficient  a  and  penumbra 
coefficient  /3  are  calculated  for  a  given  linewidth,  gap,  and  source  size.  Furthermore, 
it  is  an  easy  task  to  summarize  the  data  for  several  months  of  operation  and  to 
statistically  evaluate  the  parameters. 

The  database  was  created  with  the  help  of  Mark  Schattenburg  and  Scott  Hector. 

Dose  incident  in  the  resist 

The  incident  power  Pi,  of  the  x-rays  on  the  nitride  membrane  is  given  by 
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Figure  G-1;  Sample  of  a  run-sheet  for  an  exposure  in  helium  environment. 
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which  states  that  the  incident  x-ray  power  is  the  product  of  the  electrical  power 
Pe  times  the  brightness  of  the  source  divided  by  the  square  of  the  distance.  The 
brightness  is  given  in  units  of  W/SrW.  The  values  for  the  brightness  are  taken  from 
Henke  [100]  for  anode  voltages  of  8  kV  for  copper  and  6  kV  for  graphite,  which  are 
the  voltages  of  maximum  brightness,  and  given  in  tab.  G.l. 


Line  Source 

Brightness  [ph  /  SrW  s] 

Brightness[W /SrW] 

Cu-L 

4.8x10'" 

7.15xl0"6 

C-K  1 

1.83x10" 

8.14x10-6 

Table  G.l:  Brightness  of  X-ray  Sources. 


To  calculate  the  areal  dose,  the  power  is  multiplied  by  the  total  time  of  the  ex¬ 
posure.  For  the  volume  dose,  one  has  to  calculate  the  absorption  in  the  PMMA. 
The  mass  absorption  coefficient  in  PMMA  for  Cu-L  x-rays  is  m=3324.9  cm^/g,  the 
density  of  PMMA  is  r  =  1.2  g/cm^.  This  results  in  a  coeflScient  rm  =  3990  cm“T  We 
calculate  the  areal  dose  at  the  resist-substrate  interface  and  the  volume  dose  is  the 
difference  between  the  two  areal  doses  divided  by  the  resist  thickness.  Here  we  assume 
that  for  a  Taylor  expansion  of  the  exponential  intensity  decay,  =  l+ax+Oix/2+... 
we  can  neglect  all  but  the  constant  and  linear  term. 


Calculation  of  attenuation  in  Au,  Cr,  SiN,  Polyimide,  and  Helium  Atmo¬ 
sphere 

The  intensity  of  the  x-rays  is  proportional  to  7  oc  where  a  is  the  absorption 
coefficient  and  t  is  the  thickness  of  the  resist.  We  can  write  for  the  output  power  Po, 
after  traveling  a  distance  t: 

P  —  Pe"*^* 
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where  Pi  is  the  input  power.  Attenuations  are  usually  given  for  different  materials 
in  dB  per  micron  for  a  specific  wavelength.  ^  The  attenuation  coefficients  for  the 
Cu-L  x-rays,  used  in  our  lab,  are  given  in  tab.  G.2.  For  completeness,  the  attenuation 


Material 

A[dB///m] 

(j[l/m] 

Silicon  Nitride 

3.4517 

7.948x10^ 

Polyimide 

1.9897 

4.581x10^ 

Au 

49.694 

1.144x10^ 

Cr 

27.248  1 

6.274x10*^ 

Table  G.2:  Attenuation  coefficients  for  Cu-L  line. 


coefficients  for  the  C-K  x-rays  are  given  in  tab.  G.3. 


Material 

A[dB//im] 

a[l/m] 

Silicon  Nitride 

30.002 

6.908x10^ 

Polyimide 

1.7384 

4.003x10^ 

Au 

119.46  ^ 

2.751x10^ 

Cr 

37.020 

8.524x10® 

Table  G.3;  Attenuation  coefficients  for  C-K  line. 

For  x-ray  exposures  in  a  helium  ambient,  the  attenuation  in  helium  needs  to  be 
calculated.  The  transmission,  T,  is  given  by 


^The  conversion  can  be  done  using  the  definition  for  dB 

A[dB]  = -10  X  log(^)  , 

that  is  10  times  the  logarithm  of  output  to  input  power.  Since  we  have  dB  per  micron  given,  we  can 
equate  above  equations: 

a  =  -10®ln(10-^)  - 
m 
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where  p  is  the  density  of  the  gas.  Since  the  chemical  composition  of  air  is 


A^79O2liye(^_l)x50  5 


we  can  calculate  for  its  density,  p, 


p=(.18  + 


%02 

21 


)  X  10  ^p/cm^  . 


The  magnitude  of  p  can  be  determined  from  the  formula; 


HaSaWa 


where  S  is  the  percentage  of  each  gas  in  room  air,  and  Wa  is  the  atomic  weight  of 
each  gas.  The  values  for  p  are  tabulated  (in  cmVg)  for  Cu-L  x-rays  (A  =13.3  A),  and 
are  given  in  tab.  G.4.  For  this  particular  case,  we  can  calculate  the  magnitude  of  p. 


Element 

N 

3.9191x10^ 

0 

5.4575x10^ 

He  1 

7.1734x10^ 

Table  G.4:  Attenuation  coefficients  for  Cu-L  line  in  gases. 

79  X  p^Vj  X  14  -t"  21  X  poj  X  16  -h  50  1)  ^  Phc  ^  4 

79  X  14  -H  21  X  16  4-  50(^  -  1)  x  4 

It  should  be  noted,  however,  that  the  doses  calculated  this  way  are  not  completely 
consistent  among  the  three  different  x-ray  sources  in  our  lab.  For  head  3,  which 
features  a  vacuum  window  and  He-atmosphere,  the  doses  are  typically  lower  for  the 
same  development  time  than  doses  in  head  1.  This  draws  the  absolute  dose  informa¬ 
tion  calculated  in  this  database  into  question.  The  need  to  develop  an  absolute  dose 
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measurement  system  must  be  emphasized.  It  is  still  possible  to  compare  exposures 
done  on  the  same  x-ray  source. 

Two  additional  parameters,  a  and  /?,  were  introduced  into  the  database  [18,101, 
20].  a  is  a  measure  of  diffraction,  as  defined  in  eq.  (2.1); 


G  is  the  maximum  allowable  gap,  W  is  the  minimum  linewidth,  and  A  is  the  x-ray 
wavelength.  /3  is  a  measure  of  penumbra  5,  defined  by  the  formula 


5  =  , 


which  states  that  fl  is  the  ratio  of  the  penumbra  to  the  minimum  linewidth.  The 
penumbra  can  be  calculated  if  the  source  extent,  the  source-to-sample  distance,  and 
the  gap  are  known. 


Appendix  H 


Electrical  Measurements  Box 


In  order  to  perform  low-temperature  measurements  at  low  noise  levels,  it  is  necessary 
to  design  an  interface  between  the  probe  and  the  measurement  equipment.  The  basic 
design  of  this  interface  was  done  with  the  help  of  Arvind  Kumar  and  Nathan  Belk. 
A  good  reference  for  low  noise  measurement  systems  was  written  by  Ott  [102], 


Board  layout:  The  layout  of  the  board  is  shown  in  Fig.  H-1.  The  interface  box 
consists  of  four  gate  filters,  an  active  voltage  divider  and  two  current  amplifiers.  To 
keep  the  layout  clear,  most  of  the  circuitry  was  concentrated  on  the  backside  of  the 
board.  Only  the  operational  amplifiers  (opamps),  trimpots,  connectors,  and  dip- 
switches  were  put  onto  the  top  of  the  board  to  facilitate  easy  servicing  and  switching. 
All  wires  that  need  to  be  connected  to  the  board  are  brought  onto  the  board  through 
connectors  so  that  it  is  relatively  easy  to  take  out  and  service  the  board.  All  circuits 
are  assembled  on  one  circuit  board  and  tied  to  a  single  ground  point. 

The  four  low-pass  filters  for  the  gates  have  a  cut-off  frequency  of  3  Hz  and  have 
color  coded  wires  and  usually  use  OP27  opamps.  The  colors  used  for  the  wires  of 
the  corresponding  sub-circuits  are  indicated  on  Fig.  H-1.  A  circuit  diagram  of  the 
gate  filters  is  shown  in  Fig.  H-2.  The  active  voltage  divider,  whose  circuit  diagram 
is  shown  in  Fig.  H-3,  reduces  the  voltage  by  a  factor  of  10,000  and  filters  the  voltage 
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4  gate  filters  (3  Hz) 
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l-amp  1 
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Trimpots 

ground  _  I  l-amo  1 _ I 

I  l-amo  2  I 
I  V-divider  I 


1-8 

9-16 

17-22 

dip  switches:  open->connect,  close->ground 

23-30 

31-38 

39-44 

l-amp  2 
blue 


active  voltage 
divider  10000:1 


green 


2  9V  transistor  batteries 
lifetime  ~  30  hrs. 


Dip  switch  assignment: 

1- 1  e6  V/A  for  l-amp  1  when  closed 

2- 1  e8  V/A  for  l-amp  1  when  closed 

3-  NC 

4-  series-C  for  l-amp  1  when  open 

5- 1  e6  V/A  for  l-amp  2  when  closed 

6- 1  e8  V/A  for  l-amp  2  when  closed 

7- NC 

8-  series-C  for  l-amp  2  when  open  - 


Figure  H-1:  Layout  of  interface  box. 


Dip-switches 
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V1-V2 
1  +  sRC 


Figure  H-2:  Circuit  diagram  of  gate  filter. 


at  the  same  time  with  a  160  Hz  low-pass  filter.  The  voltage  divider  uses  OP27  as 
well;  any  offset  voltage  can  be  eliminated  using  a  trimpot.  The  current  amplifiers, 
shown  in  Fig.  H-4,  have  a  gain  of  either  10®  or  10®V/A,  which  can  be  selected  using 
the  dip  switches  in  the  right  top  corner  of  the  board.  Another  dip  switch  connects 
a  10  /LiF  capacitor  in  series  with  the  input  of  the  current  amplifier.  The  capacitor 
is  an  AC-short  for  high  device  impedances,  which  is  the  case  for  transistors  below 
threshold.  The  board  is  powered  with  two  9  V  transistor  batteries,  which  have  a 
hold  time  of  about  36  hours  under  a  load  of  all  six  opamps  as  shown  in  figure  H-5. 
The  44  leads  of  the  two  cables  are  connected  to  the  dip  switches  on  the  bottom  of 
the  board.  A  closed  dip  switch  connects  the  corresponding  lead  to  ground,  an  open 
dip  switch  floats  the  lead  that  can  now  be  set  to  any  potential.  On  the  top  right 
of  the  box  are  some  connectors,  which  can  be  used  to  calibrate  the  voltage  divider 
and  current  amplifiers.  The  two  middle  contacts  are  simply  shorted  so  that  one  can 
directly  measure  the  applied  voltage  by  connecting  the  green  wire  (output  of  voltage 
divider)  to  one  of  the  black  wires.  In  between  the  right  two  contacts  there  is  a  1  kf2 
resistor  and  between  the  left  two  contacts  a  1  Mf)  resistor,  which  can  be  used  to 
calibrate  the  current  amps  by  connecting  the  green  wire  (voltage  divider)  and  the 
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Figure  H-3:  Circuit  diagram  of  active  voltage  divider. 


Figure  H-4:  Circuit  diagram  of  current  amplifier. 


151 


Figure  H-5:  Supply  voltage  versus  time. 


blue  wire  (current  amp)  and  reading  out  the  current  value. 
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Appendix  I 


Gap  Control  Circuit 


Circuit  description 

The  circuit  that  enables  one  to  control  the  gap  between  the  mother  and  daughter 
masks  during  mask  replication  is  shown  in  Fig.  I-l.  It  is  a  straightforward  regulat¬ 
ing  circuit  using  operational  amplifiers  to  keep  the  capacitance  between  the  masks 
constant.  For  simplicity  and  availability  OP27s  were  used  as  operational  amplifiers 
(opamps),  but  more  primitive  opamps  such  as  the  OP07  could  be  used  as  well.  A 
Boonton  72B  capacitance  meter  is  used  to  measure  the  capacitance  and  its  analog 
output  is  supplied  to  the  circuit  on  pin  W.  This  voltage,  weighted  by  the  feedback 
resistors,  is  compared  to  the  reference  voltage  which  is  set  by  the  operator,  using  a 
variable  voltage  divider  and  a  Zener  diode,  at  the  central  operational  amplifier.  This 
opamp  drives  an  npn-transistor  which  regulates  the  MKS  248A  valve  to  keep  the 
capacitance  constant.  The  transistor  can  be  switched  off  manually  or  automatically 
in  case  a  short  between  the  masks  occurs.  The  short  is  detected  using  the  opamp  in 
the  top  of  the  diagram  which  acts  as  a  switch  in  a  conventional  Wheatstone  bridge 
configuration.  The  opamp  on  the  bottom  measures  the  current  going  through  the 
valve  for  diagnostic  purposes. 

Not  shown  in  this  figure  is  a  separate  capacitor  which  can  be  switched  in  parallel 
with  the  masks.  This  has  the  purpose  of  extending  the  measurement  range  of  the 
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Board  Pinout: 

A.V+ 

L.  LED  short  indicator 

M.  Mask  Contact 

N.  Mask  Contact 

P.  Switch  for  close  valve 

R.  Switch  for  close  valve 

S.  Valve  control  return 

T.  Current  output  -10mV/mA 

U.  Valve  control 

V. Pot 

W. Capacitance  Meter 

X. Pot 

Y.  ground 

Z. V- 


There  is  a  capacitance  in  series  with  the  mask  to  extend  the  range  of  the  C-meter. 
This  capacitance  is  about  9.44  nF,  nominal  value  is  10  nF. 

Valve  specs:  current  controlled,  completely  open  at  less  thanlOO  mA, 
should  be  protected  to  less  than  140  mA.  There  is  a  hysteresis,  and  the 
threshold  is  around  30  to  50  mA,  depending  on  the  setting  on  the  bottom. 

It  takes  about  1 0  to  1 5  mA  from  turn  on  to  completely  open.  The  internal 
resistance  of  the  coil  is  100  Ohms,  so  that  the  current  is  automatically 
limited  to  1 37  mA  in  this  case  because  of  V+  =  1 5  V.  1=1 5V/1 09Ohms  =  1 37  mA! 

Figure  I-l;  Diagram  of  gap-control  circuit. 
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capacitance  meter  above  3.8  nF.  The  series  capacitance  has  a  nominal  value  of  10  nF, 
and  the  capacitance  can  be  extracted  from  the  displayed  capacitance  using  C  = 
l/{l/Cdisp  —  1/10  nF).  This  way,  the  mask  capacitance  can  be  extended  to  about 
6.1  nF. 

Operation 

First,  the  daughter  mask  should  be  checked  for  particles  and  placed,  membrane  side 
up,  onto  the  o-ring  on  the  metal  base  plate,  after  being  connected  to  one  of  the 
contact  leads.  Then  the  mother  mask  should  be  checked  for  particles,  connected  to 
the  other  contact  lead  (polarity  doesn’t  matter)  and  placed  membrane  side  down 
onto  the  daughter  mask.  The  capacitance  meter  should  be  connected  to  a  strip  chart 
recorder  and  monitored  during  the  entire  run.  The  gap  can  be  monitored  easily  with 
green  light,  as  described  in  2.3.3.  After  about  a  minute,  the  masks  should  be  at  a 
gap  determined  by  the  mesa  warp  or  by  aluminum  studs,  typically  about  3-5  yum. 
To  verify  that  the  leads  make  contact  to  both  masks,  the  capacitance  should  read 
approximately  2  nF.  In  case  there  is  a  short,  the  LED  on  the  box  will  light  and  the 
capacitance  meter  will  be  off  scale.  In  this  situation,  the  mask  should  first  be  rotated 
by  a  small  amount  to  see  if  the  masks  are  still  shorted  together.  If  this  does  not  work, 
either  particles  have  to  be  removed,  or  a  different  daughter  mask  has  to  be  used.  The 
potentiometer  in  the  front  of  the  box  should  now  be  dialed  until  the  valve  starts 
to  open  up,  which  can  be  determined  by  observing  the  interference  fringes  in  green 
light.  Alternatively,  the  valve  current  can  be  monitored  using  the  second  channel 
of  the  strip  chart  recorder.  The  capacitance  can  be  increased  up  to  3.8  nF  without 
the  series  capacitor  mentioned  in  the  previous  paragraph.  The  gap  can  easily  be  set 
to  about  1  /urn.  Below  a  gap  of  1  yum  there  is  a  possibility  that  the  masks  will  come 
into  intimate  contact,  which  is  not  a  desirable  state.  During  exposure,  resist  on  the 
daughtermask  outgases  and  creates  bubbles  and  large  gaps  associated  with  it.  Contact 
exposures  typically  don’t  turn  out  successful  unless  there  is  some  provision  for  the 
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gas  to  escape.  Furthermore,  if  the  masks  are  in  intimate  contact,  they  typically  stay 
that  way  until  forced  apart,  endangering  both  the  mother  and  the  daughter  mask.  A 
state  in  which  the  masks  are  in  intimate  contact  appears  to  be  energetically  favorable 
over  a  state  in  which  the  masks  are  at  a  small  gap.  The  reasons  for  this  effect  have 
not  been  fully  explored.  To  get  the  masks  out  of  intimate  contact,  I  typically  use  thin 
wires  to  exert  a  small  lever  on  opposite  sides  of  the  masks  to  force  them  apart. 


Appendix  J 


Calculation  for  Double  Dot 
System 


This  appendix  contains  a  calculation  of  the  charging  energy  using  the  method  used  by 
Ruzin  et  al.[89].  In  order  to  keep  calculations  consistent  with  section  3.5,  the  indices 
are  adjusted  to  that  convention.  The  two  dots  and  the  charges  in  the  dot  are  shown 
schematically  in  Fig.  J-1.  The  charging  energy  U  depends  on  the  charges  Qi  on  the 
capacitors  Ci  in  the  following  way; 

2  2  2  2 

U {qag,  Qbg,  QaU  Qbl,  Qab)  =  ~  Ki^ag  +  Qbg)  ■  (J-1) 

Cag  and  Cbg  are  the  capacitances  between  the  gate  and  dot  a  and  b,  Ci  is  the  dot  to 
lead  capacitance,  and  Cab  is  the  dot  to  dot  capacitance.  Because  the  charge  in  the 
dots  is  discrete,  charge  balance  must  be  fulfilled^, 


Qal  +  Qab  -  Qag  =  -NaC 
Qbl  —  Qbg  —  Qab  =  —NbC  . 


(J.2) 


^  Ruzin  [89]  uses  qu  +  qtg  —  qab  =  —Nae,  which  is  incorrect  with  the  chosen  sign  convention. 
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Figure  J-1:  Schematic  of  the  double  dot  device,  depicting  the  important  capacitances 
in  the  double  dot  system. 


For  every  state,  the  charging  energy  should  be  at  a  minimum,  so  that  we  get 

dU/dqi  =  0,  i  =  al,bl,ab  .  (J-3) 


This  leaves  us  with  five  equations  for  five  unknowns  and  we  have  to  solve  the  following 
matrix  equation 


1  0 
0  1 
0  0 
0  0 


-1 

0 

C,  ^  Cas 


0  0  7^ 

C'ag 


0 

-1 

0 

J.  +  -L. 

Cl  Cbg 

_ 

Cbg 


Cab 


^ag 

_ 

Cbg 

4.  -L^ 

I  '  r'. 


Cbg  / 


^  Qaa  ^ 
Qbg 
Qal 
Qbl 

y  Qab  J 


eNa 
eNb 

if  i^ag 

y  -  slIiL 

Cbg 

eN„  _  eNh 
\  Cag  Cbg  ) 


■  (J.4) 


The  linear  algebra  of  this  particular  problem  turns  out  to  be  rather  complicated  and 
is  best  calculated  with  a  software  package  such  as  Mathematica.  For  Ci  —  Cab,  the 


solution  of  this  matrix  equation  is  that  of  equation  (8)-(10)  in  Ref.  [89]^.  Those 
formulas  were  given  in  eq.  3.63. 

We  are  interested  what  effect  the  coupling  capacitance  has  on  the  charging  energy. 
It  is  clear  that  the  coupling  to  the  leads  has  some  influence  on  the  energy  scales 
within  the  system,  as  it  has  for  the  single  dot  in  the  form  of  a  lever  arm.  To  get  some 
qualitative  information  about  how  the  system  behaves,  it  is  interesting  to  consider 
the  case  Q  0.  This  should  result  in  an  expression  identical  to  eq.  3.67.  This  would 
be  equivalent  to  solving  the  linear  system 


Solving  eq.  J.4  with  Q  0  yields  the  same  results  as  solving  eq.  J.5,  as  expected. 
The  solution  is 


U{Na,Nb)  = 


Ca,Nl  +  C,gNl  +  2Ca,NaN,  +  +  CggN^ 

2  CabCag  +  CabCbg  +  CagCbg 


e{Na  +  Nb)Vg  . 


For  tunneling  to  take  place,  we  must  fulfill  the  conditions 


U{N,,Nb)  =  U{Na-l,Nb) 

UiNa,Nt)  =  U{Na,Nb-l)  .  (J.7) 


This  leads  to  the  conditions: 

%  _  {Cab  +  Cbg){Na  -  \)  +  CgbNb 

e  CabCag  +  CabCbg  +  CagCbg 

_  {Cab  +  Cag){Nb  -  I)  +  CabNg 
C  CgbCgg  "1“  CgbCbg  CggCbg 

^We  again  mention  that  there  is  a  mistake  in  formula  (9)  in  I.  M.  Ruzin’s  paper.  The  index  Ci 
should  read  C-.i,  meaning  “not  i”. 


(J.8) 


160  APPENDIX  J.  CALCULATION  FOR  DOUBLE  DOT  SYSTEM 

These  conditions  have  the  right  limits,  so  that  in  the  independent  dot  limit,  Cab  ^ 

Cagt  Cbg 

Vi  =  -  i)  .  (J.9) 

Dbg  ■  2, 

and  in  the  single  dot  limit  Cab  ^  Cag,  Cbg 

=  r  Ir  +  ■  (•'■1“) 

This  means  that  for  vanishing  values  of  the  coupling  capacitance  we  expect  conduc¬ 
tance  peaks  only  if  both  conditions  J.9  are  fulfilled,  yielding  random  peaks. 
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